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The  CREAMS-WT  and  BASIN  water  quality  models  are  designed  to 
reflect  the  hydrologic,  nutrient  (nitrogen  and  phosphorus),  and 
management  characteristics  typical  of  flat,  sandy,  high-water-table 
South  Florida  flatwoods  watersheds.  Characteristics  to  be  considered 
include  1)  the  storage  based  hydrologic  system;  2)  low  chemical 
reactivity  with  the  sand  soils;  3)  nutrient  uptake  in  wetlands; 
4)  management  practices  appropriate  for  pasture  and  dairy  land  use. 

CREAMS-WT,  a  modified  version  of  CREAMS,  simulates  the  dynamic 
water  table,  limits  deep  seepage,  and  uses  a  modified  curve  number 
method  for  predicting  runoff.  Predicted  annual  water  balance  com- 
ponents compare  favorably  with  data  from  three  watersheds  in  South 
Florida  and  are  dramatically  better  than  CREAMS  predictions.  A 
20-year  simulation  shows  that  the  model  is  stable  and  performs 
predictably  over  the  range  of  natural  rainfall  variability. 


The  nutrient  component  of  CREAMS  was  changed  to  represent  the 
lower  phosphorus  buffering  capacity  of  flatwoods  soils.  Predicted 
annual  nutrient  loads  compare  well  with  data  from  three  watersheds  that 
cover  the  range  of  expected  areal  nutrient  loading  rates. 

The  BASIN  model  predicts  the  delivery  of  field/cell  nutrient  loads 
to  sub-basin  and  basin  outlets.  A  first-order  rate  equation  is  used  to 
attenuate  nutrient  loads  in  streamflow  and  flow  through  wetlands.  Rate 
coefficients  were  determined  from  nutrient  uptake  studies  reported  in 
the  literature.  BASIN  also  computes  nutrient  loads  from  cattle  in 
streams  and  wetlands,  simulates  effects  of  fencing  and  impoundments, 

and  provides  background  loads  for  non-agricultural  land  use. 

2 
The  464  km  Taylor  Creek--Nubbin  Slough  (TCNS)  basin  was  simulated 

without  calibration  using  5  ha  grid  cells  as  "fields".  Parameters  were 

selected  from  the  literature  and  from  knowledge  of  the  area.  Predicted 

long  term  average  annual  nutrient  yields  compare  closely  with  mean 

observed  values  at  the  basin  outlet  and  at  gaged  sub-basin  outlets. 

The  CREAMS-WT  and  BASIN  models  can  be  used  for  evaluating  the 

relative  effectiveness  of  alternative  management  practices  at  both  the 

field  and  basin  level  in  the  TCNS  basin,  and  should  be  applicable  to 

other  watersheds  in  Coastal  Plain  flatwoods  regions. 
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CHAPTER  I 
INTRODUCTION 


Protecting  and  managing  our  water  resources  is  an  issue  of 
national  interest.  Of  the  leading  water  resource  problems  facing  the 
nation,  nonpoint  source  stream  pollution  has  been  identified  as  one  of 
continuing  concern  (DeCoursey,  1985).  The  State  of  Florida  is  experi- 
encing increasing  stress  on  its  v/ater  resources  due  to  the  growing 
population  and  the  accompanying  land  development.  In  the  1960s, 
development  in  the  Miami  area  resulted  in  dairy  farmers  moving  from 
that  vicinity  to  the  open  range  lands  north  of  Lake  Okeechobee,  locat- 
ing in  the  Taylor  Creek  and  Nubbin  Slough  watersheds.  These  large 
dairies  grow  very  little  of  their  own  feed,  choosing  rather  to  buy  it 
from  outside  sources.  The  result  is  that  large  amounts  of  nitrogen  and 
phosphorus  are  imported  to  these  watersheds.  Not  surprisingly,  some  of 
this  nitrogen  and  phosphorus  reaches  the  watershed  outlet.  In  calcu- 
lating the  nutrient  budget  of  Lake  Okeechobee,  Federico  et  al .  (1981) 
found  that  on  an  annual  basis,  the  Taylor  Creek--Nubbin  Slough  (TCNS) 
basin  contributes  approximately  30%  of  the  phosphorus  and  5%  of  the 
nitrogen  to  the  lake  while  accounting  for  only  4%  of  the  water  inflow. 

Concern  over  the  eutrophication  status  of  Lake  Okeechobee  has 
focused  attention  on  reducing  the  nutrient  load  flowing  into  the  lake. 
Accordingly,  both  state  and  federal  funds  have  been  allocated  to  the 
TCNS  basin  with  the  goal  of  reducing  nutrient  loading  to  the  lake.  The 
state  funded  Taylor  Creek  Headwaters  Program  (TCHP)  and  the  federally 


funded  Rural  Clean  Waters  Program  (RCWP)  have  both  provided  cost 
sharing  to  land  owners  for  implementing  management  practices  which  have 
been  identified  as  having  potential  for  reducing  the  nonpoint  source 
nutrient  load  coming  from  this  basin. 

The  "Best  Management  Practice"  Concept 

Effort  has  gone  into  identifying,  often  conceptually,  the  manage- 
ment practices  which  might  be  expected  to  reduce  nonpoint  sources  of 
pollution  when  compared  with  the  standard  management.  These  "best 
management  practices"  (BMPs)  focus  on  ways  to  minimize  the  negative 
water  quality  impact  of  a  particular  land  use  while  keeping  it  profit- 
able for  the  land  user.  While  reduction  of  point  source  pollution 
loads  has  focused  on  collection  and  treatment  immediately  before 
discharge,  the  emphasis  for  reducing  nonpoint  source  loads  has  been  on 
reducing  the  pollution  potential  at  the  source  through  the  use  of  BMPs. 

With  the  emphasis  on  utilizing  BMPs  for  controlling  nonpoint 
pollution,  an  obvious  need  has  been  to  have  some  measure  of  the  effec- 
tiveness of  different  BMPs.  Proper  management  of  any  type  of  system 
depends  on  having  reasonable  estimates  of  the  expected  outcomes  of  the 
different  alternatives  being  considered.  So  too  with  water  resources, 
the  land  user  and  the  governing  water  quality  agency  must  resolve  the 
often  conflicting  management  objectives  of  maximizing  return  and 
minimizing  pollution.  A  logical  decision  can  be  reached  only  with 
adequate  data.  Estimates  of  the  effectiveness  of  BMPs  can  be  used  to 
a)  select  the  most  appropriate  BMP  for  a  particular  site;  b)  estimate 
the  overall  impact  of  BMPs  applied  in  a  watershed;  c)  rank  the 
cost-effectiveness  of  BMP  alternatives;  d)  determine  an  "optimum"  BMP 
program  based  on  a  specified  objective. 


The  objective(s)  used  in  determining  what  is  "most  appropriate"  or 
for  an  optimization  scheme  will  depend  on  the  individual  situation  but 
may  include  meeting  some  desired  water  quality  goal,  or  achieving  the 
best  possible  reduction  in  pollutant  load.  The  possible  financial 
objectives  are  numerous.  Some  examples  would  be  maximizing  the  utility 
of  cost-sharing  monies,  minimizing  cost  (initial  investment  or  mainte- 
nance costs  or  both)  to  the  farmer,  or  maximizing  return  while  meeting 
some  minimum  water  quality  standard. 

Evaluating  BMP  Effectiveness 

Two  approaches  have  been  used  for  evaluating  the  effectiveness  of 
BMPs:  watershed  data  collection  and  mathematical  modeling.  While 
measuring  the  actual  effects  may  be  the  most  desirable  way  to  evaluate 
the  response  of  a  system,  there  are  significant  drawbacks  to  using  data 
collection  for  general  BMP  evaluation.  Early  data  collection  programs 
comparing  paired  watersheds  found,  in  many  cases,  that  differences 
between  the  plots  masked  the  treatment  effects  being  evaluated 
(DeCoursey,  1985).  There  is  tremendous  variability  in  the  watershed 
characteristics  that  affect  the  rainfall-runoff  process  and  chemical 
transformations.  A  large  number  of  sampling  sites  (watersheds)  would 
be  required  to  get  a  statistically  significant  result,  and  this  is 
prohibitively  expensive.  Time  is  a  second  limiting  factor.  Due  to  the 
variability  in  rainfall  and  other  meteorological  data,  it  is  desirable 
to  have  a  long  term  record.  Not  only  is  it  expensive  to  maintain  a 
sampling  program,  but  we  cannot  wait  years  to  get  these  data  before 
taking  action  on  the  water  resource  issues  that  face  us  today. 

The  second  approach  to  BMP  evaluation  is  the  use  of  mathematical 
models.  Physical  process  models  attempt  to  describe  the  physical, 


chemical  and  biological  processes  that  define  the  natural  system  being 
modeled.  Models  developed  around  our  understanding  of  physical 
processes  enable  us  to  examine  in  greater  detail  the  various  inter- 
relationships between  those  processes.  The  other  main  class  of  models, 
empirical  models,  are  less  applicable  for  use  in  BMP  evaluation. 
Empirical  models  are  mathematical  formulations  based  on  data  analysis 
and  do  not  necessarily  represent  the  actual  processes  controlling  the 
system.  They  should  not  be  used  beyond  the  range  of  data  used  in  their 
development,  and  in  addition,  empirical  models  often  use  fitted 
parameters  thus  are  not  well  suited  for  simulating  BMPs. 

Mathematical  Models 

Although  long  term  watershed  data  are  needed  for  model  development 
and  testing,  it  is  now  recognized  that  the  primary  means  of  evaluating 
the  effectiveness  of  BMPs  is  through  the  use  of  detailed  water  quality 
models  designed  for  that  purpose  (Beasley  et  al . ,  1982;  DeCoursey, 
1985).  The  impact  of  a  BMP  on  offsite  water  quality  is  determined 
through  modeling  by  simulating  the  pre-BMP  and  post-BMP  conditions, 
and  comparing  the  results.  BMPs  are  simulated  by  making  changes  in  the 
model  parameter  values  that  are  representative  of  the  actual  changes 
that  would  be  expected  were  the  BMP  actually  implemented. 

However,  models  also  have  their  limitations.  Extensive  amounts  of 
data  are  often  required,  some  of  which  will  not  be  readily  available. 
In  some  cases,  more  is  known  about  the  governing  processes  than  can  be 
implemented  in  a  model  due  to  computational  (cost)  and  data  limita- 
tions. We  are  also  limited  in  our  understanding  of  the  details  of  the 
processes  that  govern  the  natural  system.  There  is  a  tremendous  amount 
of  variability  in  the  response  of  any  natural  system,  and  while  current 


physically-based  models  may  be  able  to  predict  mean  response  with  some 
certainty,  much  of  the  variability  is  not  addressed  as  the  causal 
processes  are  not  fully  understood  (McDowell  et  al.,  1980). 

In  practice,  BMPs  are  implemented  on  a  field  by  field  basis;  thus, 
the  most  accurate  assessment  of  the  effectiveness  of  a  BMP  can  be  ob- 
tained with  field-scale  simulation.  However,  it  may  also  be  desirable 
to  evaluate  the  effects  of  a  BMP  implementation  program  over  a  water- 
shed or  basin  to  provide  information  for  evaluating  broader  planning, 
funding  and  water  quality  objectives. 

Research  Problem  and  Objective 

Management  practices  with  potential  for  reducing  the  nonpoint 
nutrient  load  discharged  from  the  TCMS  basin  have  been  identified 
(Ritter  and  Allen,  1982).  Under  the  TCHP  and  RCWP  programs,  the  BMPs 
deemed  most  effective  have  been  implemented  (RCWP,  1982;  1984),  and 
already,  the  question  of  the  effectiveness  of  this  basin-wide  BMP 
program  is  being  debated.  The  water  quality  data  collection  effort 
faces  the  problems  discussed  above  in  attempting  to  categorize  the 
effectiveness  of  the  BMP  program.  A  water  quality  model  could  provide 
some  indication  of  the  long  term  water  quality  improvement  that  might 
be  expected  from  various  BMP  programs. 

Bottcher  and  Baldwin  (1981)  developed  a  BMP  application  model 
which  uses  a  digitized  database  of  the  TCNS  basin  and  matches  BMPs  with 
appropriate  land  use.  Selection  of  BMPs  is  based  on  land  use,  cattle 
density,  and  distance  from  waterways.  With  the  inclusion  of  cost  data, 
the  model  can  be  used  to  determine  the  total  amounts  of  BMPs  applied 
and  the  associated  costs  for  different  basin-wide  BMP  scenarios.  For  a 
more  complete  evaluation  of  BMP  alternatives,  Bottcher  and  Baldwin 


(1981)  recognized  the  need  for  a  water  quality  model  which  could  pro- 
vide estimates  of  the  water  quality  improvement  associated  with  the  BMP 
scenarios. 

Numerous  water  quality  models  have  been  designed  for  general 
application.  However,  their  use  in  evaluating  BMP  effectiveness  in  the 
TCNS  basin  (and  in  similar  South  Florida  flatwoods  watersheds)  is  sig- 
nificantly limited  by  several  somewhat  unique  features  of  this  region. 

1)  The  hydrologic  response  is  dominated  by  the  flat  terrain,  sandy 
soils,  and  a  high  water  table. 

2)  Phosphorus  transport  and  transformations  in  the  sandy  soils  do 
not  follow  conventional  concepts. 

3)  Wetlands  comprise  a  significant  portion  of  the  area  and  often 
receive  (by  design  and  default)  nutrient-ladened  runoff. 

4)  This  region  is  predominantly  pasture  with  very  little  tilled 
land,  and  appropriate  BMPs  may  not  be  readily  reflected  in  the 
parameters  of  the  existing  models  which  tend  to  be  oriented 
towards  row  crop  land  use. 

The  objective  of  this  research  is  to  develop  a  water  quality  model 
that  can  be  used  to  evaluate  the  effectiveness  of  BMPs  in  South  Florida 
flatwoods  watersheds.  The  model  should  encompass  both  field  and  basin 
scale  simulation  and  must  address  the  unique  hydrologic,  chemical,  and 
management  conditions  characteristic  of  these  watersheds. 


CHAPTER  II 
LITERATURE  REVIEW 


Description  of  the  Taylor  Creek  —  Nubbin  Slough  Basin 
General  Characteristics 

The  Taylor  Creek  —  Nubbin  Slough  (TCNS)  basin  lies  in  the  Southern 
Florida  Flatwoods  physiographic  region  (Figure  2-1).  The  region  is 
very  flat  with  sandy  soils.  The  term  "flatwoods"  is  still  used  even 
though  in  most  areas  the  pine  forests  have  been  cleared  for  range  and 
pasture.  The  TCNS  basin  is  very  similar  to  the  larger  adjoining  Lower 
Kissimmee  River  (LKR)  basin  which  has  similar  physiographic  and  land 
use  characteristics. 

The  soils  of  this  area  are  generally  classified  as  fine  sand, 
often  with  less  than  1%  clay  content,  and  typically  have  very  high 
hydraulic  conductivities  (>16  cm/hr) .  Internal  drainage  may  be  very 
rapid  to  slow,  so  a  dual  hydrologic  classification  is  used,  A/D  or  B/D. 
with  the  appropriate  classification  dependent  on  what  drainage  improve- 
ments have  been  made.  Slopes  range  from  0  to  2  percent  but  are 
typically  less  than  0.5  percent.  As  a  result  of  the  low  slopes,  there 
is  little  or  no  erosion. 

Water  tables  vary  with  topography  and  soil  classification  but 
generally  range  from  the  surface  to  2  meters  deep.  The  water  table  in 
the  phreatic  aquifer  is  underlaid  by  the  impermeable  Hawthorne 
Formation,  and  beneath  that  is  the  artesian  Floridan  aquifer.  Thus, 
deep  seepage  is  generally  assumed  to  be  negligible  (Stewart,  1980; 
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Figure  2-1.  Florida  Flatwoods  Dhysiographic  regions  and  site  location. 


Yates  et  al . ,  1982).  During  the  wet  summer  months,  the  water  table  is 
near  the  surface  and  is  the  controlling  factor  in  the  hydrologic 
response  of  the  area. 

Wetlands  are  an  important  component  of  the  basins,  impacting  both 
hydrology  and  water  quality.  There  are  many  isolated  ponded  areas, 
some  seasonal  and  many  year-round.  Due  to  the  lack  of  sufficient 
gradient,  these  areas  often  remain  wet  even  when  ditched.  The  primary 
drainage  ways  are  poorly  defined  channels  through  marshes  and  wooded 
wetlands.  Total  wetland  area  in  the  TCMS  basin  accounted  for  approx- 
imately 8%  of  the  total  area  in  1980. 

Land  use  in  the  area  has  moved  towards  more  intensive  management, 
with  dairy  farming  and  ranching  now  dominating.  For  the  Upper  Taylor 
Creek  watershed,  Allen  et  al .  (1982b)  report  1959  land  use  as  24% 
improved  pasture,  44%  unimproved  pasture,  17%  forest  and  range,  and  15% 
miscellaneous  (including  roads,  buildings,  and  ponded  wetlands).  By 
1980  it  was  3%  citrus,  69%  improved  pasture,  16%  forest  and  range,  4% 
urban,  and  8%  miscellaneous.  The  Nubbin  Slough  basin  has  similar  land 
use,  while  the  larger  Lower  Kissimmee  River  basin  is  less  developed, 
with  higher  percentages  of  unimproved  pasture  and  range  land. 

The  basins  are  ditched  in  varying  degrees  depending  on  the  land 
use.  Ditching  densities  increase  from  unimproved  pasture,  to  improved 
pasture,  with  intensively  managed  pasture  along  with  citrus  and  row 
crops  having  the  greatest  ditching  densities  (18  km/sq  km  reported  by 
Huber  et  al . ,  1976).  Ditches  often  help  blur  watershed  boundaries  as 
they  inter-connect  over  wide  areas.  Watershed  boundaries  are  generally 
poorly  defined,  and  direction  of  flow  in  a  ditch  may  be  dependent  on 
prevailing  wind,  type  of  storm,  and  the  depth  of  water  in  the  ditch. 
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Several  factors  uniquely  impact  the  water  quality  of  this  region. 
There  is  evidence  that  phosphorus  will  leach  through  sandy  flatwoods 
soils  (Neller  et  al . ,  1951;  Blue,  1970;  Chaiwanakupt  and  Robertson, 
1976).  This  is  a  rare  phenomenon  for  phosphorus  in  soils.  Very 
reactive  in  solution,  phosphorus  normally  is  readily  adsorbed  by  soil 
minerals  with  negligible  concentrations  remaining  in  the  soil  solution 
(Buckman  and  Brady,  1969;  Larsen,  1967).  Because  of  the  very  low  clay 
content,  the  adsorptive  capacity  of  these  soils  is  greatly  reduced 
resulting  in  significant  phosphorus  concentrations  in  the  soil  water 
and  in  runoff. 

Wetlands  play  an  important  role  in  water  quality.  They  are  most 
commonly  considered  as  sinks  for  nutrients  although  under  certain  con- 
ditions they  may  be  a  source  of  nutrients.  A  number  of  studies  have 
been  done  in  the  LKR  basin  (Davis,  1982;  Federico  et  al . ,  1978; 
Goldstein,  1986),  but  pertinent  data  can  also  be  obtained  from  studies 
along  the  Eastern  Coastal  Plain  (Simpson  et  al . ,  1983;  Craig  and 
Kuenzler,  1983)  and  Louisiana  (Conner  and  Day,  1982). 

Hydrologic  Processes 

Runoff  from  these  flat,  sandy  soils  is  not  produced  by  the  common 
concept  of  an  infiltration  limiting  'rainfall  excess1.  Runoff  is  more 
likely  a  result  of  rainfall  on  saturated  areas  and  subsurface  flow  to 
ditches  and  drainage  ways  both  during  and  following  a  storm.  This 
concept  of  the  rainfall/runoff  process  is  very  similar  to  the  'variable 
source  area'  concept  (Hewlett  and  Nutter,  1970;  Hewlett  and  Troendle, 
1975;  Betson  and  Marius,  1969).  These  researchers  were  studying  forest 
hydrology  where  infiltration  rates  were  also  rarely,  if  ever,  a  limiting 
factor. 
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The  very  flat  terrain  also  challenges  some  of  the  concepts  of 
classical  hydrology.  The  standard  terms  of  overland  flow,  interflow, 
and  base  flow  are  difficult  to  apply  to  flatwoods  watersheds.  Accepted 
interpretations  associate  these  terms  with  particular  flow  processes 
and  specific  time  intervals  of  discharge  hydrographs.  Although  the 
same  flow  processes  occur  on  flatwoods  watersheds,  they  are  often 
difficult  to  separate  and  not  associated  with  the  same  hydrograph  time 
interpretations.  Direct  runoff  is  typically  defined  as  including  over- 
land flow  and  interflow.  Base  flow  is  defined  as  being  nearly  constant 
and  originating  from  groundwater.  More  appropriate  terms  for  describ- 
ing runoff  from  flatwoods  watersheds  are  rapid,  intermediate,  and  slow 
flow.  These  terms  do  not  attempt  to  explain  process  but  refer  only  to 
flow  rate  (Speir  et  al.,  1969). 

A  water  balance  for  this  area  shows  that  the  rainfall  and  soil 
moisture  storage  leave  the  system  primarily  by  evapotranspiration  and 
secondarily  by  runoff  with  percolation  or  deep  seepage  being  negligible 
(Stewart,  1980;  Yates  et  al.,  1982).  Thus,  the  level  of  the  water 
table  (i.e.  change  in  soil  moisture  storage)  directly  reflects  input 
from  rainfall,  and  losses  by  evapotranspiration  and  runoff. 

Data  Sources 
There  are  numerous  sources  of  data  describing  the  hydrology  and 
water  quality  of  Southern  Florida  Flatwoods.  The  Upper  Taylor  Creek 
watershed  has  been  studied  in  greatest  detail,  with  initial  data 
collection  beginning  in  1955  (Speir  et  al.,  1969).   In  the  past  20 
years,  concern  over  the  nutrient  status  of  Lake  Okeechobee  has  spawned 
a  spate  of  special  studies  specifically  structured  to  survey  nutrient 
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dynamics  in  the  Lower  Kissimmee  River  flatwoods.  Additional  data  can 
be  gathered  from  studies  of  other  watersheds  in  the  flatwoods  of 
Florida  (Riekerk  et  al.,  1979)  and  the  Eastern  Coastal  Plain  (Craig  and 
Kuenzler,  1983). 

Taylor  Creek—Nubbin  Slough 

Initial  research  interests  in  the  Upper  Taylor  Creek  watershed 
were  to  determine  rainfall-runoff-evapotranspiration  relationships  for 
natural  watersheds  in  Central  Florida  (Knisel  et  al.,  1985).  This  was 
broadened  to  include  analysis  of  the  effects  of  channel  improvements 
and  water  control  structures  on  storm  runoff,  water  yield  and  ground 
water  in  the  1960 ' s .  Beginning  in  1972,  water  quality  sampling  was 
begun  at  15  sites.  Additional  flow  gaging  and  water  quality  sampling 
sites  have  been  added  in  the  intervening  years.  Three  agencies  have 
cooperated  in  collecting  data:  the  USDA-ARS,  USGS,  and  the  South 
Florida  Water  Management  District  (SFWMD). 

Most  of  the  data  collected  in  this  basin  are  from  the  Upper  Taylor 
Creek  watershed  (also  identified  as  watershed  W-2) ,  and  from  sub- 
watersheds  within  this  watershed,  Northwest  Taylor  Creek  (W-3) , 
Williamson  Ditch  (W-5)  and  Otter  Creek  (W-13)  (See  Figure  2-2).  In 
June  1973,  Upper  Taylor  Creek  was  diverted  into  the  Nubbin  Slough  canal 
to  be  discharged  into  Lake  Okeechobee  through  structure  S-191.  The 
drainage  area  contributing  to  the  flow  at  S-191  has  since  been  referred 
to  collectively  as  the  Taylor  Creek--Nubbin  Slough  basin. 

Comprehensive  reports  on  the  hydrology  and  hydrogeology  of  the 
Upper  Taylor  Creek  watershed  (W-2  and  W-3)  have  been  published  by  the 
USDA-ARS  (Speir  et  al . ,  1969;  Knisel  et  al . ,  1985).  Additional  reports 
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Figure  2-2.  Taylor  Creek— Nubbin  Slough  basin  and  sub-watersheds. 
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have  focused  on  specific  aspects  of  the  data:  hydrology  (Allen  et  al., 
1976;  Allen  et  al . ,  1982a),  effects  of  channel  modifications  (Yates  et 
al.,  1982),  and  water  quality  (Stewart  et  al . ,  1978;  Allen  et  al . , 
1982b). 

With  the  initiation  of  the  TCHP  and  RCWP  cost  sharing  programs, 
the  focus  has  turned  to  water  quality,  and  additional  sampling  was 
started  in  the  lower  portion  of  the  basin.  Water  quality  sampling  and 
analysis  has  been  done  by  the  SFWMD,  with  samples  routinely  collected 
on  a  biweekly  basis.  Samples  have  been  collected  at  all  flow  gaging 
stations  and  at  numerous  other  locations  in  the  watershed  in  an  attempt 
to  categorize  the  nutrient  concentrations  that  result  from  different 
land  use  practices.  Ritter  and  Allen  (1982)  provide  a  detailed  report 
and  analysis  of  the  water  quality  data  from  Upper  Taylor  Creek.  Data 
collected  since  publication  of  that  report  (and  including  data  from 
Nubbin  Slough  and  at  S-191)  are  summarized  in  the  RCWP  Annual  Progress 
Reports  (RCWP,  1982;  1983;  1984). 

There  are  two  problems  that  plague  collection  and  analysis  of  flow 
data  from  these  watersheds.  First,  due  to  the  very  low  slopes  in  some 
channels  accurate  flow  measurement  is  difficult.  Particular  uncer- 
tainty exists  in  the  streamflow  measurement  at  the  W-2  outlet,  and  flow 
conditions  may  also  affect  the  chemical  concentrations  in  the  channel 
(Knisel  et  al . ,  1985).  Thus,  calculated  nutrient  loads  from  this 
watershed  have  been  noted  as  being  "estimated"  (RCWP,  1984).  Knisel  et 
al.  (1985)  consider  "observed"  concentrations  and  calculated  loads  to 
be  in  the  right  order  of  magnitude,  but  caution  that  uncertainty  in  the 
"observed"  values  must  be  kept  in  mind. 
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A  second  problem  in  hydrologic  analysis  of  these  watersheds  is 
what  appears  to  be  the  dynamic  nature  of  the  watershed  boundaries. 
This  is  due  to  the  low  slopes  and  the  widespread  use  of  ditches  for 
land  drainage.  The  areas  of  all  the  gaged  watersheds  in  the  TCNS  basin 
(with  few  exceptions),  have  been  changed  during  the  course  of  the  data 
collection  program  as  documented  by  Knisel  et  al.  (1985)  and  in  the 
RCWP  reports.  While  most  of  the  major  changes  in  ditching  and  drainage 
patterns  have  surely  been  identified,  the  very  fact  that  this  occurs  is 
indicative  of  the  uncertainty  that  exists  in  defining  the  watershed 
area. 

Upland  Detention/Retention  Demonstration  Project 

Data  have  been  collected  from  several  small  watersheds  in  the  LKR 
and  TCNS  basins  under  the  Upland  Detention/Retention  Demonstration 
Project  (UD/RDP)  implemented  by  the  SFWMD  (Goldstein,  1986).  The 
objectives  of  this  project  were  to  obtain  information  on  the  water 
quality  problems  associated  with  typical  agricultural  practices  in 
these  basins,  and  to  evaluate  the  potential  of  using  wetland  areas  for 
pollution  abatement.  Data  from  this  project  have  been  valuable  in 
helping  to  characterize  hydrologic  processes  on  the  field  scale,  and  in 
identifying  the  water  quality  response  of  individual  land  uses. 

Some  characteristics  of  these  watersheds  are  listed  in  Table  2-1, 
with  their  locations  shown  in  Figure  2-3.  This  data  collection  program 
was  maintained  for  approximately  three  years  (Sept.  1979  through  1982) 
with  measurement  of  precipitation,  runoff,  water  table  elevation,  and 
biweekly  samples  taken  for  water  quality  analysis.  Unfortunately, 
1980-1981  was  one  of  the  driest  periods  on  record.  Not  only  were  there 
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Table  2-1.  Characteristics  of  the  Upland  Detention/Retention 
Demonstration  Project  watersheds  (Capece,  1984). 


Parameter  Armstrong  Peavine    SEZ    Ash  W.   Ash  E. 

Slough   Pasture   Dairy  Pasture  Pasture 


Drainage  area  (ha) 
Channel  slope  (%) 
Overland  slope  (%)   ~ 
Drainage  dens,  (km/km  ) 
Ponds  and  marsh  (%) 
Length/width  ratio 


fewer  data,  but  how  well  any  of  the  data  represent  "typical"  conditions 
must  be  questioned. 

Analysis  of  the  runoff  data  from  several  of  these  sites  has  raised 
questions  as  to  their  integrity.  In  some  cases  runoff  exceeds  rainfall 
even  over  a  period  of  several  months.  In  addition  to  the  typical  prob- 
lems of  defining  the  watershed  boundary  in  this  flat  region,  there  were 
problems  with  maintenance  of  the  sites  (Capece,  1984).  The  estimated 
range  of  expected  errors  in  flow  measurement  were  5  to  15  percent,  with 
the  flumes  designed  for  optimum  measurement  of  low  to  normal  flow 
events  (Goldstein,  1986).  Thus,  as  runoff  volumes  increase,  measure- 
ment errors  will  also  increase  considering  the  backwater  conditions 
that  tend  to  occur.  This  may  partially  explain  the  unreasonably  high 
runoffs  "measured"  in  1982. 

Data  from  the  three  smaller  sites  were  used  in  the  development  and 
testing  of  the  field  scale  model  and  will  be  described  briefly  here. 
Additional  description  of  the  sites  can  be  found  in  the  report  by 
Goldstein  (1986)  and  from  the  hydrologic  analysis  of  Capece  (1984). 
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Ash  Slough  (Bass  West  and  Bass  East).  These  two  pastures  on 
Myakka  fine  sand  are  on  either  side  of  the  Ash  Slough  wetland  (Figure 
2-4).  Both  pastures  are  used  for  grazing  beef  cattle.  The  East  pas- 
ture drainage  boundary  is  poorly  defined,  with  the  contributing  area 
estimated  to  be  between  8  and  20  ha.  The  West  pasture  has  a  deep  pe- 
rimeter ditch  with  an  internal  network  of  shallow  (2  ft)  ditches. 
Cattle  movements  resulted  in  numerous  breaches  in  the  low  levee  around 
the  pasture  and  no  estimate  of  possible  flow  errors  could  be  made.  In 
addition,  there  was  no  adjustment  for  possible  subsurface  flow  into  the 
perimeter  ditch  from  "outside"  the  watershed. 

SEZ  dairy.  The  SEZ  dairy  watershed  is  a  291  ha  area  surrounded  by 
a  deep  ditch  (Figure  2-5).  The  soils  are  predominantly  fine  sands 
(Immokalee,  Myakka),  with  approximately  7%  of  the  watershed  being  wet- 
lands. As  a  dairy,  part  of  the  area  (56  ha)  is  holding/staging  pasture 
for  the  dairy  herd,  with  217  ha  used  for  heifer  and  dry  cow  pasture  and 
hayland.  The  remaining  18  ha  is  the  farmstead  area. 

Lower  Kissimmee  River  Basin 

There  have  been  several  studies  in  the  Lower  Kissimmee  River  basin 
that  have  related  primarily  to  water  quality.  Federico  (1982)  charac- 
terized the  water  quality  of  the  Lower  Kissimmee  River  canal  (C-38)  and 
its  tributaries  with  monthly  sampling  over  a  five  year  period  (1973-78). 
Federico  et  al.  (1978)  report  on  environmental  studies  in  the  Chandler 
Slough  marsh.  A  study  in  the  Boney  Marsh  was  designed  to  provide  data 
on  the  potential  use  of  wetland  areas  for  removing  nutrients  from 
runoff  (Davis,  1981).  Pertinent  information  from  these  studies  will  be 
included  in  following  sections. 
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Figure  2-3.  Location  of  the  Upland  Detention/Retention 
Demonstration  Project  watersheds. 
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Phosphorus  Transformations  in  Sandy  Soils 
Phosphorus  Retention  in  Soil 

It  has  long  been  recognized  that  when  soluble  phosphate  is  added 
to  the  soil  it  rapidly  disappears  from  the  soil  solution.  In  his 
review  article,  Wild  (1949)  noted  that  this  phenomenon  had  been  demon- 
strated 100  years  earlier.  Volumes  of  reports  and  related  phosphorus 
(P)  research  have  followed  in  the  past  130  years.  It  has  been  sug- 
gested that  the  subject  of  P  retention  in  soils  has  prompted  more 
research  than  any  other  aspect  of  soil-fertilizer-plant  interactions 
(Sample  et  al.,  1980). 

Along  with  the  hundreds  of  studies  demonstrating  rapid  P  retention 
by  the  soil,  there  have  also  been  a  few  that  did  not  observe  this 
phenomenon  but  found  that  P  would  actually  pass  through  the  soil 
profile  in  significant  amounts.  Phosphorus  leaching  has  been  reported 
in  sandy  soils  in  Florida  (Neller,  1946;  Blue,  1970;  and  others)  and  in 
Australia  (Hingston,  1959;  Russell,  1960;  Ozanne  et  al . ,  1961)  and  is 
somewhat  unusual  in  soil/phosphorus  interactions.  This  response  of 
some  sandy  soils  can  lead  to  problems  with  water  quality  as  a  result  of 
excessive  P  in  ground  and  surface  waters,  concerns  that  do  not  apply  to 
most  soils. 

Soluble  phosphorus  added  to  the  soil  as  fertilizer  is  often  ren- 
dered unavailable  even  under  the  most  ideal  field  conditions.  In 
general,  the  overall  phosphorus  problem  for  agriculture  is  1)  the 
small  total  amount  present  in  soils,  2)  the  unavailability  of  this 
native  P,  and  3)  the  rapid  sorption  of  added  soluble  P  (Buckman  and 
Brady,  1969).  Thus,  the  fixation  of  P  in  soils  is  most  commonly  viewed 


22 

as  a  problem  because  it  decreases  the  efficiency  of  fertilization, 
rapidly  rendering  the  applied  P  unavailable  to  plants  (Sanchez  and 
Uehara,  1980).  This  is,  in  general,  directly  opposite  the  concerns  of 
P  management  in  the  TCNS  basin. 

The  availability  of  inorganic  P  is  primarily  controlled  by:  1) 
soil  pH,  2)  soluble  Fe,  Al ,  and  Mn,  3)  minerals  containing  Fe,  Al ,  and 
Mn,  4)  available  calcium  and  Ca  containing  minerals,  5)  amount  and 
decomposition  of  organic  matter,  and  6)  microorganism  activity  (Buckman 
and  Brady,  1969).  The  first  four  are  related,  largely  through  their 
dependence  on  pH.  Under  low  pH  conditions,  P  is  fixed  by  Fe  and  Al 
compounds,  and  in  high  pH  calcareous  soils,  Ca  is  the  primary  fixing 
agent.  In  a  more  recent  review,  Sanchez  and  Uehara  (1980)  define 
factors  affecting  the  amounts  of  P  fixed  in  terms  of  soil  charac- 
teristics: clay  mineralogy,  clay  content,  x-ray  amorphous  colloid 
content,  exchangeable  Al,  and  soil  organic  matter.  They  state  that 
these  factors  are  not  strictly  additive  in  their  effects  but  interact 
with  each  other. 

Modeling  Soil/Phosphorus  Reactions 

The  phosphorus  present  in  the  soil  has  been  categorized  as  being 
in  three  states:  phosphate  in  the  soil  solution,  adsorbed  phosphate 
which  is  in  some  state  of  equilibrium  with  the  solution,  and  nonlabile 
P  which  is  more  firmly  held  (Larsen,  1967).  A  variety  of  reactions 
have  been  proposed  as  those  maintaining  the  state  of  equilibrium 
between  the  soil  P  pools,  with  precipitation  and  adsorption  receiving 
the  most  attention.  The  chemistry  of  these  soil  phosphorus  reactions 
is  described  by  Larsen  (1967),  Olsen  and  Khasawneh  (1980)  and  others. 
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Barrow  (1980)  presents  evidence  that  on  the  addition  of  P  to  a 
soil,  there  is  a  two  step  reaction.  The  first  step,  the  adsorption  of 
P,  is  rapid,  with  over  90%  of  the  P  commonly  fixed  within  the  first 
hour  of  contact  (Sanchez  and  Uehara,  1980).  The  second  step  is  a  slow 
one  in  which  some  of  the  rapidly  adsorbed  P  is  converted  to  more  firmly 
held  forms.  These  reactions  may  last  for  days  or  months  and  are 
suggested  to  be  the  diffusion  of  P  into  inner  surfaces  of  hydroxides  or 
a  rearrangement  of  the  crystal  structure  (Barrow,  1974).  Sanchez  and 
Uehara  (1980)  also  note  that  precipitates  formed  with  Al  and  Fe,  though 
highly  insoluble,  will  gradually  release  phosphate  ions  into  the  soil 
solution  over  several  years. 

Hemwall  (1957),  after  summarizing  precipitation  and  adsorption 
studies  concluded  that  regardless  of  the  type  of  the  reaction,  the 
compounds  formed  and  the  mechanism  of  reaction  seem  to  be  essentially 
the  same.  Now  there  is  widespread  acceptance  of  the  concept  of  a 
continuum  of  reaction  mechanisms  in  a  media  as  diverse  as  soil.  There 
is  little  concern  for  distinguishing  between  precipitation  and  adsorp- 
tion reactions,  and  both  phenomena  are  generally  considered  together  as 
"sorption"  (Berkheiser  et  al.,  1980). 

The  phosphorus  isotherm  defines  the  equilibrium  relationship 
between  the  P  content  of  the  soil  and  the  concentration  in  the  soil 
solution.  Phosphorus  sorption  in  soils  has  most  frequently  been  char- 
acterized using  the  Freundlich  or  the  Langmuir  adsorption  isotherms. 
However,  many  additional  equations  have  been  proposed  (Barrow,  1978; 
Berkheiser  et  al.,  1980),  some  being  modifications  of  the  Langmuir  or 
Freundlich  equations.  There  has  been  a  great  deal  of  controversy  in 
discussing  the  particular  advantages  of  these  two  models  and  in 
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arguing  for  one  as  being  'best'.  However,  theoretical  distinctions 
between  the  two  have  been  narrowed,  and  Kurtz  (1981)  concludes  that  it 
makes  little  difference  which  equation  is  used  as  both  have  empirically 
fitted  parameters. 

The  adsorption/desorption  process  is  not  completely  reversible, 
with  desorption  occurring  at  a  slower  rate  than  adsorption  (Selim  et 
al.,  1975;  Enfield  and  Ellis,  1983).  The  Langmuir  and  Freundlich 
isotherms  assume  complete  reversibility  thus  do  not  adequately  describe 
the  adsorption/desorption  process.  The  rate  of  desorption  of  P  has 
been  found  to  decrease  proportional  to  the  cube  root  of  elapsed  time 
(Barrow,  1980).  Research  on  desorption  is  largely  from  the  viewpoint 
of  plant  availability  of  residual  P  and  is  and  not  concerned  with  the 
source  of  the  P.  Thus,  it  is  not  clear  if  this  release  of  P  to  the 
soil  solution  comes  from  adsorbed  P  or  from  more  tightly  held  P.  No 
distinct  discontinuity  has  been  observed  in  the  rate  of  release;  thus 
no  clear  differentiation  on  the  source  of  the  P  can  be  made  (Barrow, 
1980). 

With  the  growing  need  for  predicting  nonpoint  sources  of  phos- 
phorus, algorithms  describing  P  desorption  are  being  developed  (Chien 
and  Clayton,  1980;  Sharpley  et  al.,  1981;  Sharpley  et  al.,  1985). 
Sharpley  (1983)  evaluated  several  models,  relating  physical  and 
chemical  properties  of  60  soils  to  the  models  parameters.  Sharpley 
concludes  that  these  equations  may  have  practical  application  in 
predicting  soluble  P  release  in  modeling  water  quality,  particularly 
the  model  by  Sharpley  et  al .  (1981)  whose  parameters  were  found  to  be 
proportional  to  a  soil's  percent  clay/organic  C  content  ratio,  thus 
readily  determinable. 
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Phosphorus  Relationships  in  Sandy  Soils 

The  factors  conducive  to  leaching  are  essentially  the  absence  of 
the  soil  constituents  that  are  responsible  for  P  fixation.  Thus 
leaching  of  P  has  been  observed  primarily  in  sandy  soils  with  low  clay 
content  (typically  less  than  5%)  and  is  primarily  correlated  to  low 
content  of  extractable  Al  and  Fe  (Ballard  and  Fiskell,  1974;  Humphreys 
and  Pritchett,  1971;  Yuan  and  Lucas,  1982;  Chaiwanakupt  and  Robertson, 
1976).  Barrow  (1980)  also  attributes  the  low  water  retention  capacity 
of  sands  as  a  contributory  factor  to  leaching  of  P,  presumably  because 
water  may  move  through  the  profile  faster,  and  there  is  less  reaction 
time  for  P  in  solution. 

Despite  the  ready  leaching  of  P  from  flatwoods  soils  as  described 
above,  it  has  also  been  shown  that  P  can  accumulate  in  the  surface  soil 
in  significant  amounts  provided  a  higher  pH  is  maintained  with  regular 
liming.  Blue  (1970)  found  that  with  an  average  soil  pH  of  6.0,  70%  of 
the  P  added  to  a  Leon  fine  sand  over  a  period  of  18  years  was  retained 
in  the  surface  soil  of  the  white  clover-perennial  grass  pasture.  In 
similar  pastures,  phosphorus  contents  of  156  and  327  ppm  have  been 
found  (Rodulfo  and  Blue,  1970).  Spencer  (1957),  working  with  Lakeland 
fine  sand  planted  to  citrus,  found  that  over  a  15  year  period  where  the 
pH  had  been  maintained  at  5.6  or  above,  26.4%  of  the  added  P  had 
accumulated  in  the  surface  6  inches,  but  where  the  pH  had  been  allowed 
to  drop  to  4.0,  no  P  was  retained.  This  build-up  of  P  has  also  been 
observed  over  a  five  year  period  in  the  top  three  inches  of  an 
Immokalee  fine  sand  with  pH  greater  than  5.5  (Neller  et  al . ,  1951). 


26 

Even  though  the  relationship  between  liming  and  P  retention  in 
sandy  soils  has  been  observed  repeatedly,  the  mechanism  involved  in 
this  P  retention  has  not  been  clearly  identified.  Under  the  conditions 
of  high  Ca  content,  pH  approaching  neutrality  (6.0)  and  low  Fe  and  Al 
in  the  profile,  it  is  expected  that  the  P  would  be  held  as  Ca  compounds 
(Blue,  1970).  However,  this  has  not  been  observed.  Yuan  et  al.  (1960) 
analyzed  for  water-soluble,  aluminum,  iron  and  calcium  phosphates 
fractions  in  three  Florida  soils.  In  a  Leon  fine  sand  (pH  5.4)  to 
which  100  ppm  P  was  added,  90%  of  the  P  was  found  associated  with  Al 
and  7%  with  Fe,  with  the  water-soluble  and  Ca  forms  negligible. 

Thus,  it  is  uncertain  how  much  P  may  accumulate  in  a  soil  under 
these  conditions.  In  a  normal  P  fixing  soil,  it  has  been  observed  that 
a  decrease  in  P  sorption  was  related  to  previous  P  applications  (Smyth 
and  Sanchez,  1980),  which  is  consistent  with  the  concept  of  adsorption 
sites  becoming  saturated.  For  sandy  soils  with  lime  added,  such 
assumptions  are  less  sure  because  of  the  uncertainty  of  the  mechanism 
of  P  retention.  The  soil  P  content  may  increase  with  continued  liming, 
but  this  P  is  poorly  buffered  and  may  be  released  readily. 

Release  of  P  held  in  limed  sandy  soils  is  most  often  associated 
with  a  decrease  in  the  pH  when  liming  is  halted  (Neller  et  al.,  1951; 
Spencer,  1957;  Forbes  et  al.,  1974).  Rodulfo  and  Blue  (1970)  have  also 
shown  that  a  Leon  fine  sand  with  high  P  content  (327  ppm)  had  suffi- 
cient P  available  to  plants  so  that  no  supplemental  fertilizer  P  was 
needed.  The  logical  conclusion  is  that  when  the  soil-P  pools  are  in  a 
relatively  stable  state,  changing  one  aspect  of  the  balance  (such  as 
decreasing  the  pH)  will  affect  the  others. 
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Nutrient  Uptake  in  Wetland  Systems 
Among  the  management  practices  suggested  for  reducing  nonpoint 
source  pollution  loading  to  receiving  waters  is  the  possibility  of 
using  wetlands  as  a  water  purifier.  Using  wetlands  as  nutrient  sinks 
has  been  proposed  as  a  BMP  for  South  Florida  range  and  dairy  land 
(Allen  et  al.,  1982b).  Wetlands  have  been  used  for  municipal  waste- 
water polishing  for  several  years  (Boyt  et  al . ,  1977;  Ti 1  ton  and 
Kadlec,  1979).  However,  this  use  of  wetlands  has  not  gone  without 
debate.  The  use  of  wetlands  for  filtering  of  agricultural  runoff  must 
be  re-evaluated  in  light  of  recent  reports.  Can  wetlands  effectively 
filter  nutrients  (in  particular,  phosphorus)  from  runoff?  What  type  of 
wetland  systems  might  be  appropriate  for  such  use,  and  what  uptake 
rates  should  be  considered? 

To  understand  nutrient  movement  and  storage  in  wetlands,  Valiela 
and  Teal  (1978)  suggest  study  on  two  scales:  first,  the  whole  eco- 
system, looking  at  inputs,  outputs,  and  internal  processes,  and  second, 
analyzing  in  depth  the  component  processes  of  the  system.  Component 
processes  that  have  been  studied  are  nitrogen  and  phosphorus  transfor- 
mations, interactions  with  sediment,  and  nutrient  uptake  and  release  by 
macrophytes.  System  studies  have  been  done  with  natural  wetlands  and 
with  artificial  wetland  reservoirs. 

Component  Studies 

Nitrogen  transformations  in  the  water  column  and  sediment 

Soluble  nitrogen  in  floodwater  occurring  as  ammonia  and/or  nitrate 

can  undergo  several  transformations.  Ammonia  in  floodwater  can  undergo 
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volatilization,  assimilation  by  plants,  nitrification,  and  can  be 
adsorbed  and  retained  in  the  soil  (Chen  and  Patrick,  1981;  Reddy  and 
Graetz,  1981).  In  an  aerobic  water  column,  ammonia  is  readily  oxidized 
to  nitrate  completing  one  side  of  the  nitrification/denitrification 
process  (Patrick  and  Reddy,  1976).  In  a  water  column  with  high  pH  (8.5 
up),  ammonia  volatilization  can  take  place  and  nitrification  is  slowed 
(Reddy,  1983). 

Nitrate  transformations  include  immobilization  by  incorporation 
into  organic  forms,  plant  assimilation,  reduction  of  nitrate  to  ammon- 
ia, and  denitrification.  Denitrif ication  has  been  shown  to  be  the 
primary  pathway  of  nitrate  removal  from  soil  samples  incubated  with 
floodwater  (Bartlett  et  al . ,  1979;  Reddy  et  al . ,  1980;  Reddy  and 
Graetz,  1981;  Sompongse,  1982).  This  is  true  even  at  lower  tempera- 
tures where  there  is  an  increase  in  nitrate  reduction  to  ammonia 
(Bartlett  et  al.,  1979;  Reddy  et  al.,  1980),  and  immobilization  to 
organic  N  (Reddy  et  al.,  1980).  Denitrification  is  effective  for  both 
mineral  and  organic  soils  (Krottje,  1980;  Reddy  and  Graetz,  1981),  and 
in  laboratory  studies  has  accounted  for  over  95%  of  the  nitrate  lost 
(Reddy  et  al.,  1980;  Sompongse,  1982). 

The  rate  of  denitrification  is  a  function  of  floodwater  depth, 
nitrate  concentration,  presence  of  organic  carbon  as  an  energy  source, 
pH,  and  temperature  (Reddy  et  al.,  1980;  Krottje,  1980).  Since  deni- 
trification of  wastewater  occurs  when  nitrate  moves  down  into  the 
anaerobic  soil  profile,  the  rate  of  diffusion  of  nitrate  into  the  soil 
is  often  the  limiting  factor  (Engler  et  al.,  1976;  Reddy  and  Graetz, 
1981).  Adequate  residence  time  of  the  wastewater  is  crucial  to  nitrate 
removal,  especially  if  ammonia  is  also  present  in  the  wastewater. 
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Rapid  nitrification  of  the  ammonia  to  nitrate  may  initially  increase 
the  concentrations  of  nitrate  over  several  days  (Reddy  and  Graetz, 
1981).  Reddy  and  Graetz  (1981)  found  a  12  to  24  day  residence  time 
necessary  for  effective  nitrate  removal  using  agricultural  drainage 
water  incubated  on  undisturbed  cores  from  wetland  soils.  However, 
actual  nitrate  removal  rates  and  the  required  residence  time  will  vary 
greatly,  depending  in  part  on  the  factors  mentioned  above. 
Phosphorus  transformations  in  the  water  column  and  sediment 

Various  processes  contribute  to  the  equilibrium  concentration  of 
soluble  phosphorus.  Phosphorus  is  removed  from  wetland  water  through 
plant  uptake  and  by  adsorption  and  precipitation.  Desorption  can 
release  sediment  bound  P  into  solution.  There  is  no  process  that 
removes  P  from  the  system  comparable  to  the  nitrogen  processes  of 
ammonia  volatilization  and  denitrification,  a  fact  having  important 
consequences  in  considering  the  ultimate  fate  of  P  added  to  a  wetland. 

There  are  many  factors  that  control  adsorption-desorption  (and 
precipitation).  Three  that  recur  in  the  literature  are  redox  poten- 
tial, pH,  and  availability  of  reactive  compounds  of  calcium,  aluminum 
and  iron.  Redox  potential  affects  the  release  of  P  to  the  water.  The 
general  characterization  is  that  P  is  released  into  solution  under 
anaerobic  conditions  (low  redox  potential)  (Ponnamperuma,  1972; 
Berkheiser  et  al.,  1980;  Sloey  et  al.,  1978).  In  reality,  the  inter- 
action of  other  factors  is  very  important.  Reddy  and  Graetz  (1981) 
found  that  an  organic  soil  with  an  aerobic  water  column  increased 
soluble  P  concentration,  while  the  anaerobic  reservoir  sorbed 
phosphorus.  Sompongse  (1982)  found  that  P  was  removed  from  nutrient 
enriched  water  by  sorption  under  both  aerobic  and  anaerobic  conditions. 
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The  effect  of  pH  on  phosphate  adsorption  appears  to  show  no  uni- 
fied trend.  However,  it  is  recognized  as  an  important  factor:  directly 
through  effects  on  chemical  reactions,  and  indirectly  through  effects 
on  plant  decomposition  (Berkheiser  et  al.,  1980;  Sloey  et  al . ,  1978). 

A  third  dominating  factor  in  P  sorption  is  the  presence  of  Ca,  Fe, 
and  Al  in  the  water  and  sediment.  The  availability  of  Ca  compounds  is 
strongly  related  to  precipitation  of  P  (Sloey  et  al.,  1978;  Reddy  and 
Graetz,  1981),  and  Al  and  Fe  oxides  in  the  sediment  are  the  primary 
sorption  agents  for  P  (Berkheiser  et  al.,  1980;  Sompongse,  1982). 
Richardson  (1985)  found  that  the  phosphorus  adsorption  potential  in  a 
variety  of  wetland  soils  could  be  predicted  solely  by  the  extractable 
Al  content  of  the  soil.  The  generally  low  Al  and  Fe  content  of  organic 
soils  means  they  have  a  limited  capacity  for  removing  P  from  the  over- 
lying water.  Thus,  while  there  may  be  good  uptake  initially,  removal 
rates  decline  quickly  as  the  system  becomes  saturated,  and  P  export  may 
eventually  occur  (Kadlec  and  Hammer,  1982;  Richardson,  1985). 

The  variability  in  wetland  response  is  demonstrated  in  the  uptake 
of  P  in  wetlands  with  organic  soils.  Several  studies  have  found 
organic  soils  to  be  a  source  of  soluble  P  under  anaerobic  conditions 
(Peverly,  1982;  Reddy  and  Graetz,  1981;  Reddy  et  al.,  1982b).  Reddy 
and  Graetz  (1981)  characterized  flooded  organic  soils  as  having  poor 
adsorptive  capacity.  They  found  equilibrium  phosphorus  concentration 
(EPC)--the  concentration  at  which  P  is  neither  adsorbed  or  desorbed— to 
be  2.25  mg/L  for  a  flooded  organic  soil,  and  0.05  mg/L  for  a  calcareous 
marl  soil  used  to  line  a  reservoir.  Other  researchers  have  noted  good 
potential  for  adsorption  in  organic  soils  (Tilton  and  Kadlec,  1979; 
Sloey  et  al.,  1978).  It  appears  difficult  to  characterize  in  simple 
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terms  the  sorption  of  P  in  a  wetland  system.  The  exact  physico- 
chemical  nature  of  the  water  and  sediment  control  the  response,  and 
changes  in  one  or  more  factors  could  dramatically  affect  that  response. 
Nutrient  uptake/release  by  aquatic  macrophytes 

Following  the  pioneering  work  of  Seidel  (1976)  in  Germany  in  the 
mid-1950s,  there  has  been  growing  interest  in  the  ability  of  aquatic 
macrophytes  to  purify  wastewater.  Some  of  the  species  that  have  been 
studied  are  waterhyacinth,  bulrush,  cattail,  reeds,  rush,  bamboo, 
pennywort  and  elodea  (Lakshman,  1979;  Wolverton  et  a!.,  1983;  Reddy, 
1983;  Ogwada,  1983).  Macrophytes  affect  the  nutrient  dynamics  of  the 
wetland  system  through  nutrient  assimilation,  nutrient  release  by 
decomposition,  altering  the  physico-chemical  environment  of  the  water 
column  and  through  sediment/water  nutrient  cycling. 

Macrophytes  can  alter  the  physico-chemical  environment  of  the 
water  affecting  nutrient  transformations.  Denitrification,  ammonia 
volatilization  and  precipitation  of  ortho-phosphate  can  be  increased 
significantly  as  a  result  of  plant  presence  creating  anaerobic 
conditions  (Reddy,  1983). 

Prentki  et  al . ,  (1978)  studied  nutrient  movement  in  a  Wisconsin 
marsh  that  did  not  have  external  nutrient  inputs.  Macrophytes  rooted 
in  the  sediment  withdrew  nutrients  during  the  growing  season,  returning 
those  nutrients  to  the  water  and  sediment  surface  upon  decomposition. 
In  this  habitat,  nutrient  cycling  by  aquatic  plants  was  found  to  be  a 
significant  source  of  nutrients  in  the  water. 

However,  most  of  the  interest  in  studying  macrophytes  is  for 
nutrient  uptake.  Nutrient  removal  by  macrophytes  depends  on  the 
climate,  type  of  plant,  and  nutrient  availability.  The  ability  of 
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macrophytes  to  remove  nutrients  from  municipal  and  agricultural  waste- 
water has  been  well  documented  in  micro-environments,  reservoirs,  and 
natural  wetlands  (Boyd,  1970;  Spangler  et  al.,  1976;  Lakshman,  1979; 
Reddy,  1983;  DeBusk  et  al.,  1983). 

The  waterhyacinth  has  attracted  most  attention  because  of  its 
prodigious  growth  rate  and  accompanying  large  uptake  of  nutrients.  The 
waterhyacinth  takes  up  nutrients  in  proportion  to  the  concentrations  in 
the  water.  Tissue  concentrations  of  37  g  N/kg  and  9.4  g  P/kg  have  been 
observed  for  plants  grown  in  sewage  effluent.  Average  tissue  concen- 
trations for  plants  taken  from  19  eutrophic  water  bodies  in  Florida  had 
1.61  +  5.0  g  N/kg  and  3.1  +  1.8  g  P/kg  (Reddy  and  Sutton,  1984). 

Nutrient  uptake  by  waterhyacinth  under  ideal  conditions  in  a 
nutrient  rich  solution  has  been  documented  up  to  5350  kg  N/ha/yr  and 
1260  kg  P/ha/yr  (Reddy  and  Tucker,  1983).  Removal  rates  for  natural 
stands  have  been  estimated  as  480  kg  N/ha/yr  and  93  kg  P/ha/yr  (Reddy 
and  Sutton,  1984).  In  organic  soil  drainage  water  at  Zellwood, 
Florida,  730  kg  N/ha/yr  and  159  kg  P/ha/yr  were  removed  (Reddy  et  al . , 
1982a).  OeBusk  et  al .  (1983)  measured  waterhyacinth  nutrient  uptake 
from  secondary  sewage  effluent  (9.96  mg/1  total  N,  and  4.68  mg/1  total 
P)  between  May  and  August.  This  period  corresponds  to  50%  of  the 
annual  biomass  yield  in  Florida  (Reddy  and  Sutton,  1984)  giving 
projected  annual  uptakes  of  648  kg  N/ha/yr  and  207  kg  P/ha/yr. 

Plant  assimilation  of  nutrients  only  occurs  during  the  growing 
season,  and  plant  decay  releases  many  of  these  nutrients  back  to  the 
water.  Decay  rates  and  nutrient  fate  have  been  studied  under  a  variety 
of  conditions  for  different  plants,  and  frequently,  a  rapid  release  of 
a  major  fraction  of  N  and  P  has  been  observed  (Simpson  et  al.,  1978; 
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Hill,  1979;  Ogwada,  1983).  Regeneration  of  N  and  P  varies  greatly 
between  different  species  and  is  dependent  on  temperature,  pH  and 
water  oxygen  concentration  (Jewell,  1971;  Davis  and  van  der  Valk,  1978; 
Reddy  and  Sacco,  1981;  Sompongse,  1982).  Depending  on  these  factors 
(but  primarily  on  the  species),  a  small  fraction  of  plant  material  may 
resist  decomposition  and  retain  N  and  P  (Jewell,  1971;  DeBusk  et  al . , 
1983). 

Rates  of  detritus  accumulation  have  not  been  well  documented. 
Estimates  of  annual  accumulation  range  from  3.5  -  70  kg/ha  for  N,  and 
0.05  -  2.4  kg/ha  P  (Richardson,  1985).  After  17  weeks  of  decomposition 
of  waterhyacinth,  DeBusk  et  al.  (1983)  observed  a  fairly  constant  rate 
of  detritus  accumulation.  Assuming  a  six  month  growing/senescence 
season,  those  rates  would  project  to  accumulation  of  25  kg/ha  N  and  1.8 
kg/ha  P  per  year.  These  values  are  clearly  very  low  in  terms  of  a 
annual  budget  and  are  only  slightly  higher  than  the  average  loadings 
from  rainfall  (19.5  kg/ha  N  and  1.2  kg/ha  P)  in  the  LKR  and  TCNS  basins 
(Goldstein,  1986). 

Many  researchers  have  concluded  that  plant  nutrient  uptake  is  pri- 
marily a  temporary  storage  factor  in  the  nutrient  dynamics  of  wetlands 
(Prentki  et  al.,  1978;  Simpson  et  al.,  1983;  Hill,  1979).  However, 
this  time  lag  can  be  very  beneficial  in  reducing  summer  algal  blooms  in 
lakes  receiving  wetland  water  (Lee  et  al.,  1975).  Harvesting  is  the 
best  approach  for  permanent  removal  of  plant  nutrients  and  is  important 
for  maintaining  optimum  growth  and  nutrient  uptake  in  artificial 
wetland  systems  (Boyd,  1970;  DeBusk  et  al . ,  1983;  Ogwada,  1983). 

In  considering  the  role  of  nutrient  removal  by  macrophytes,  only 
two  conditions  result  in  permanent  removal:  physical  removal  of  plant 
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material  through  harvesting  and  detritus  accumulation.  In  considering 
the  annual  nutrient  budget,  the  only  process  with  a  significant  impact 
is  harvesting. 

System  Studies 
Natural  wetlands 

The  fate  of  nutrients  in  natural  wetlands  on  an  ecosystem  basis, 
looking  at  inputs,  outputs  and  sometimes  internal  processes,  has  been 
the  subject  of  numerous  reports  in  the  past  ten  years.  Investigators 
have  examined  the  use  of  wetlands  as  a  final  treatment  of  municipal 
wastewater  (Fetter  et  al.,  1978;  Boyt  et  al.,  1977;  Sloey  et  al., 
1978;  Tilton  and  Kadlec,  1979).  Agricultural  nonpoint  source  runoff 
flowing  through  wetlands  has  been  studied  in  Minnesota  (Hill,  1979), 
New  York  (Peverly,  1982),  Delaware  (Simpson  et  al.,  1983),  Louisiana 
(Conner  and  Day,  1982),  and  Florida  (Kangas,  1977;  Davis,  1982; 
Goldstein,  1986). 

High  nutrient  removal  efficiencies  characterize  the  majority  of 
municipal  wastewater  applications.  Tilton  and  Kadlec  (1979)  tabulate 
nine  such  studies  with  most  showing  nutrient  reductions  of  above  90%, 
and  only  one  case  each  of  reductions  of  ortho-P  and  ammonia  being  below 
10%.  Sloey  et  al.  (1978)  present  a  similar  summary  but  also  include  a 
few  cases  showing  less  promising  results.  In  general,  results  of 
municipal  wastewater  "treatment"  using  wetlands  has  been  positive  as 
indicated  in  the  literature. 

Studies  on  the  fate  of  nutrients  from  agricultural  runoff  in 
wetlands  have  been  more  variable  in  outcome.  Peverly  (1982)  studied  a 
3000  ha  muckland  in  New  York  over  two  years.  The  upper  region  is 
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drained  and  cultivated,  and  the  drainage  water  and  runoff  drains 
through  the  lower  region  which  is  maintained  as  a  natural  wetland.  The 
wetland  acted  as  an  N  and  P  sink  only  during  the  second  year,  when  flow 
was  about  half  that  of  the  first  year.  Otherwise,  N,  P,  K,  C,  and  Ca 
were  exported  from  the  wetland. 

Conner  and  Day  (1982)  found  that  a  Louisiana  swamp  was  effective 
in  the  removal  of  nitrate-N  from  agricultural  runoff,  but  that  on  the 
average,  ortho-P,  organic-N  and  organic-P  were  added  to  the  water. 
During  the  winter  the  oxygen  concentration  in  the  water  increased,  and 
the  increase  in  redox  potential  led  to  phosphorus  removal  by  the 
sediment. 

There  have  been  several  investigations  in  the  Kissimmee  River 
basin  in  southern  Florida.  Davis  (1982)  reports  that  the  142  ha  Boney 
Marsh  functioned  as  a  year-round  phosphorus  sink.  Phosphorus  storage 
in  both  living  and  dead  plant  material  continued  to  increase  over  the 
three  year  study.  An  average  P  retention  of  4200  kg/ha/yr  was  calcu- 
lated, a  figure  Davis  (1982)  says  compares  closely  with  the  uptake  of 
Chandler  Slough  as  calculated  by  Federico  et  al . ,  (1978).  Since  the  P 
input  to  Chandler  Slough  was  six  times  the  Boney  Marsh  input,  Davis 
(1982)  suggests  this  uptake  rate  may  be  the  maximum  assimilative 
capacity  for  this  region. 

Davis  (1982)  noted  that  nutrient  uptake  took  place  over  a  short 
distance  after  which  the  concentrations  were  reduced  to  background 
levels.  This  nutrient  uptake  front  has  moved  downstream,  with  a  perma- 
nent change  in  plant  community  noted  behind  the  front  (S.M.  Davis  per- 
sonal communication,  1985).  A  similar  progression  of  a  nutrient  removal 
front  has  been  documented  by  Kadlec  and  Hammer  (1982).  The  implication 
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is  that  any  significant  long  term  storage  must  be  reflected  by  a  change 
in  the  flora  to  a  community  which  contains  higher  nutrient  content  per 
unit  area.  The  desirability  of  such  a  change  in  the  plant  community  is 
a  factor  that  must  be  considered  in  any  management  decision. 

Goldstein  (1986)  calculated  nutrient  budgets  for  two  wetland  sys- 
tems, Ash  Slough  and  Armstrong  Slough,  and  concluded  that  these  systems 
appeared  to  be  in  equilibrium  with  nutrient  loading  from  the  surround- 
ing pastures.  Reduction  in  exported  nutrient  load  was  attributed  to 
the  reduction  in  outflow  volume  over  inflow,  and  one  year  with  high 
outflows  had  a  net  export  of  total  N  and  total  P. 

Allen  et  al .  (1982b)  report  on  a  12  ha  wetland  used  as  a  second 
stage  impoundment  for  dairy  barn  wash  water  which  first  enters  an 
anaerobic  lagoon.  Outflow  from  the  second  impoundment  had  high  total  P 
(9.55  mg/1)  and  ortho-P  (5.12  mg/1)  concentrations.  They  concluded 
that  the  receiving  area  was  too  small  and  the  high  water  and  nutrient 
inflow  overtaxed  the  wetland. 

For  all  natural  wetland  systems,  the  indications  are  that  if 
regular  inputs  are  continued  over  a  long  term,  a  steady  state  response 
will  be  reached.  ("Steady  state"  is  used  here  as  being  descriptive  of 
the  average  long  term  response.  Variability  due  to  meteorologic  inputs 
will  still  certainly  be  evident.)  This  apparent  equilibrium  has  been 
observed  in  the  wetlands  that  have  been  receiving  runoff  from  agri- 
cultural areas  over  many  years  (Peverly,  1982;  Conner  and  Day,  1982; 
Goldstein,  1986).  Wetlands  where  significant  uptake  has  been  observed 
after  initiating  wastewater  application  have  with  time  demonstrated  an 
apparent  nutrient  saturation  where  active  uptake  no  longer  occurs 
(Kadlec  and  Hammer,  1982;  Davis,  1982). 
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Artificial  wetlands 

Flooded  reservoirs  stocked  with  aquatic  plants  have  been  used  for 
waste  treatment  in  Europe  for  many  years  (Seidel,  1976;  DeJong,  1976), 
and  have  been  evaluated  for  purification  of  municipal  (DeBusk  et  al . , 
1983)  and  agricultural  water  (Reddy  et  al.,  1982a, b).  These  artificial 
wetlands  utilize  the  processes  of  natural  wetlands  without  endangering 
or  disturbing  those  resources,  and  in  addition,  enable  better  control 
and  management  of  the  system. 

Research  has  focused  on  quantifying  the  response  of  the  system, 
identifying  the  important  nutrient  pathways,  and  determining  optimum 
management  techniques.  Again,  results  of  studies  reflect  findings  for 
a  particular  set  of  conditions,  and  caution  must  be  used  in  general 
application  of  the  results.  Many  of  the  studies  looking  at  nutrient 
uptake  by  macrophytes  as  discussed  previously,  were  done  using  arti- 
ficial wetlands,  and  are  directly  applicable  to  this  topic. 

Permanent  removal  via  plant  uptake  in  these  systems  is  dependent 
on  frequent  harvesting  of  the  plants.  This  removes  nutrients  and  also 
enables  optimum  growth  to  continue.  Plant  uptake  has  been  found  to  be 
the  most  significant  pathway  of  P  removal  (DeBusk  et  al.,  1983;  Reddy 
et  al.,  1982a),  although  Reddy  (1983)  demonstrated  that  under  the  right 
physico-chemical  conditions,  P  could  also  be  removed  by  sorption  to  the 
sediment.  Denitrification  can  equal  or  exceed  plant  uptake  in  nitrogen 
removal  (Reddy,  1983;  DeBusk  et  al.,  1983). 

Artificial  wetland  reservoirs  have  been  found  to  effectively 
remove  N  and  P  from  wastewater.  Regular  harvesting  is  necessary,  and 
proper  management  is  crucial.  The  system  must  be  maintained  at  desired 
operating  conditions  to  be  effective  (Reddy,  1983). 
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Summary 

Two  methods  of  presenting  experimental  results  have  been  encoun- 
tered which  need  some  comment  in  that  they  tend  to  obscure  the  "bottom 
line."  A  nutrient  budget,  with  total  average  annual  input  and  output 
loadings,  is  the  only  way  to  determine  if  a  wetland  is  acting  as  a  sink 
or  a  source  (Richardson  et  al.,  1978).  Data  presented  as  concentra- 
tions without  flow  information  say  nothing  about  net  nutrient  flux,  as 
lower  output  concentration  may  simply  be  a  result  of  dilution. 

A  second  problem  arises  in  presenting  only  the  percentage  of 
reduction.  This  is  less  than  adequate  since  it  says  nothing  about  out- 
fall concentrations  and  total  load  exported,  which  are  more  important 
criteria  (Richardson  et  al.,  1978).  The  results  are  doubly  obscure 
when  only  the  percent  reduction  in  concentration  is  reported— an  all 
too  frequent  occurrence. 

In  contrast  to  the  enthusiasm  shown  by  some  early  proponents,  more 
caution  is  now  being  expressed  regarding  the  use  of  natural  wetlands  as 
wastewater  filters.  Caution  is  due  for  the  following  reasons. 

1.  Wetland  ecosystems  are  vastly  different  depending  on  type  of 
wetland,  climate,  soil,  plant  population,  hydrologic  regime. 
Tilton  and  Kadlec  (1979)  appropriately  warn  that  their  results 
apply  only  to  the  specific  wetland  being  studied.  Each  site 
must  be  considered  separately  with  careful  testing  and  eval- 
uation (Whigham,  1982;  Stearns,  1978;  Sloey  et  al.,  1978). 
Loading  rate  is  a  key.  A  small  load  in  a  very  large  area  may 
be  a  fraction  of  the  total  natural  nutrient  assimilation  and 
ultimately  be  negligible  even  on  a  long  term  basis. 
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2.  Most  research  has  been  short-term  studies,  and  long-range 
effects  are  largely  unknown  (Stearns,  1978;  Sloey  et  al . , 
1978).  To  determine  the  extended  impact  on  all  flora  and  fauna 
will  require  very  detailed  monitoring,  and  some  current 
information  suggests  that  the  risk  of  damage  may  not  be  worth 
the  gain  (Stearns,  1978). 
In  a  system  view  of  a  natural  wetland,  inflowing  nutrients  can  be 
removed  from  the  system  through  volatilization  and  plant  removal,  can 
be  stored  in  the  sediment,  organic  matter,  and  plants,  or  can  be 
exported  from  the  system.  Wetlands  can  unequivocally  function  as  a 
nitrogen  sink  as  a  result  of  denitrif ication.  The  fate  of  phosphorus  is 
less  certain.  If  it  is  stored,  for  what  length  of  time  can  the  wetland 
continue  to  assimilate  phosphorus?  Because  of  the  dynamics  of 
adsorption/desorption  and  the  natural  fluctuations  in  the  hydrologic 
regime  of  a  wetland,  phosphorus  export  will  most  likely  occur,  whether 
annually,  intermittently  or  continuously,  and  will  then  have  a  greater 
storage  supply  from  which  to  draw. 

Thus  it  is  difficult  to  try  to  characterize  a  given  wetland  as  a 
source  or  a  sink.  Because  of  the  higher  natural  nutrient  levels, 
wetlands  often  have  higher  background  nutrient  output.  Nutrient  output 
is  very  sensitive  to  the  hydrology  of  the  system,  to  draining  (Lee  et 
al.,  1975;  Klopatek,  1978)  and  variations  in  rainfall.  A  compounding 
problem  is  that  the  agricultural  runoff  of  greatest  concern  is  the 
large  storms,  for  which  the  treatment  capacity  of  a  wetland  is  greatly 
reduced  because  of  the  decreased  residence  time.  Management  of 
artificial  systems  is  possible,  for  control  of  inputs,  maintenance  of 
desirable  conditions,  harvesting,  and  control  of  the  hydrologic  regime. 
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The  use  of  artificial  wetlands  has  the  potential  of  capitalizing  on  the 
benefits  of  marshes  without  the  prospect  of  damaging  natural  wetlands. 

Mathematical  Models  for  Hydrology  and  Water  Quality 
The  Curve  Number  Method  for  Predicting  Runoff 

The  SCS  runoff  equation  was  developed  to  provide  a  simple  means  of 
estimating  runoff,  while  being  able  to  reflect  differences  due  to  soils, 
landform,  land  use,  and  management  practices.  Due  to  its  simplicity 
(it  is  a  single  parameter  model)  and  its  in-house  use  by  the  USDA-SCS, 
the  method  has  become  widely  used  for  estimating  runoff.  The  history 
of  the  development  of  the  equation  has  recently  been  presented 
(Rallison,  1980). 

The  SCS  runoff  equation  is  written  as 


(P  -  0.2S)2 

Q  -  [2-1] 

(P  +  0.8S) 


where  Q  is  runoff,  P  is  rainfall,  and  S  is  a  watershed  storage  param- 
eter, with  all  having  units  of  depth  (ins).  The  equation  was  intended 
to  be  used  for  individual  storms  but  is  often  used  with  daily  rainfall 
amounts  since  those  are  readily  available.  The  single  parameter,  S,  is 
usually  determined  by  way  of  a  "curve  number,"  a  transformation 
introduced  to  linearize  the  value  of  the  parameter.  The  curve  number, 
CN,  is  related  to  S  by  the  relationship: 


1000 

CN  =  — -  [2-2] 

S+10 
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The  SCS  has  published  tables  of  values  of  curve  numbers  based  on 
different  land  use,  management  practices  and  other  conditions. 

While  enjoying  widespread  use,  some  researchers  have  questioned  the 
reliability  of  the  curve  number  method.  It  has  been  found  that  the 
method  does  not  respond  to  rainfall  intensity  (Hawkins,  1978)  or  storm 
duration  (Rallison,  1980),  and  that  the  relationship  between  initial 
abstraction  and  S  used  in  the  derivation  of  the  equation  may  be 
dependent  on  location  (Hjelmfelt,  1980)  and  storm  characteristics 
(Rallison,  1980). 

Capece  (1984)  evaluated  several  methods  of  predicting  storm  runoff 
volume  to  determine  their  applicability  for  use  on  flatwoods  water- 
sheds. Data  for  his  analysis  were  from  the  Upland  Detention/Retention 
Demonstration  Project  (Goldstein,  1986).  The  runoff  prediction  methods 
analyzed  were  all  forms  of  the  SCS  runoff  equation  (equation  2-1).  The 
difference  between  the  various  methods  lies  in  the  approach  used  to 
determine  S. 

The  three  most  accurate  methods  related  S  to  an  estimate  of  the 
actual  physical  storage  available  in  the  profile,  using  curves  that 
relate  available  profile  storage  as  a  function  of  the  water  table 
depth.  The  "ARS"  curve  (Figure  2-6)  was  developed  from  Taylor  Creek 
watershed  data  (Speir  et  al.,  1969).  The  South  Florida  Water 
Management  District  (SFWMD,  1983)  developed  a  curve  for  storage  on 
natural  sites,  along  with  an  alternate  curve  for  developed  sites 
(Figure  2-7).  The  ARS  curve  predicts  less  available  storage  than  does 
the  SFWMD  curve,  with  the  differences  becoming  very  significant  for 
water  table  depths  greater  than  two  feet.  Capece  (1984)  found  that 
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Figure  2-6.  Water  storage  characteristics  of  sandy  soils  from  the 
Upper  Taylor  Creek  watershed  (Speir  et  al.,  1969). 
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Figure  2-7.  Available  soil  profile  storage  as  a  function  of  depth 
to  the  water  table  for  watersheds  in  South  Florida 
(SFWMD,  1983). 
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directly  predicting  available  storage  using  the  ARS  curve  gave  the  most 
accurate  runoff  prediction. 

Water  Quality  Models 

Numerous  water  quality  models  have  been  developed  in  the  past  ten 
years.  Many  of  these  have  been  designed  for  specific  situations,  but 
there  are  several  models  which  have  been  designed  for  general  use,  and 
have  been  adequately  documented,  distributed  and  supported.  There  are 
several  compilations  of  hydrologic/water  quality  models  which  include 
both  specialized  and  general  application  models  (Huber  and  Heaney, 
1982;  Renard  et  al.,  1982). 

There  are  several  models  which  are  applicable  for  use  in  eval- 
uating the  impact  of  various  management  practices  on  offsite  water 
quality.  Models  developed  for  this  purpose  have  generally  been 
physically  based,  distributed  parameter  models.  These  models  have 
parameters  that  are  based  on  physical  characteristics  of  the  system  so 
parameter  values  can  be  obtained  by  measurement  rather  than  through 
calibration  (Beasley  et  al.,  1982).  Several  of  the  currently  available 
models,  and  some  just  recently  developed  will  be  described  briefly. 

A  field  scale  model  for  Chemicals,  Runoff,  and  Erosion  from 
Agricultural  Management  Systems  (CREAMS),  developed  by  a  group  of  USDA 
scientists,  is  designed  specifically  for  evaluating  the  relative 
effects  of  different  management  practices  (Knisel,  1980).  Conceptually 
it  is  a  distributed  parameter,  physically  based  model,  and  requires  no 
calibration.  Hydrology,  sediment,  nitrogen  and  phosphorus  transfor- 
mations and  pesticides  are  modeled  continuously  for  up  to  twenty  years 
using  a  daily  time  step.  A  "field"  is  defined  as  an  area  of  uniform 
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hydrologic  and  chemical  characteristics  under  a  single  land  use. 
CREAMS  has  been  widely  distributed  and  used  (Bengtson  and  Carter,  1983; 
Foster  and  Ferreira,  1981;  DelVecchio  and  Knisel,  1982;  Crowder  et  al . , 
1985). 

Several  models  have  been  developed  which  rely  heavily  on  the 
CREAMS  algorithms.  The  Simulator  for  Water  Resources  in  Rural  Basins 
(SWRRB)  is  a  basin  scale  model  which  builds  on  the  daily  rainfall 
hydrology  component  of  CREAMS  (Williams  et  al . ,  1985).  The  model  adds: 
a  return  flow  component,  spatial  linking  of  watersheds  within  a  basin, 
reservoir  storage,  simple  flood  routing,  and  a  weather  simulation 
model.  In  addition  to  water,  SWRRB  simulates  sediment  yield,  but  no 
nutrients  are  considered. 

The  Erosion-Productivity  Impact  Calculator  (EPIC)  simulates 
erosion  and  its  effects  on  the  long  term  productivity  of  soils 
(Williams  et  al.,  1983,  1984).  This  is  a  field  scale  model,  based  on 
the  CREAMS  hydrology  component.  It  has  a  much  more  complex  nutrient 
model  than  does  CREAMS.  Where  CREAMS  has  a  very  simple  phosphorus 
model,  EPIC  contains  a  comprehensive  soil  phosphorus  cycle,  with 
phosphorus  pools  and  reactions  calculated  separately  in  up  to  ten  soil 
layers  (Jones  et  al.,  1984a, b;  Sharpley  et  al.,  1984).  However,  this 
model  has  not  been  released  for  public  use. 

The  AGNPS  model  (Young  et  al.,  1985)  uses  CREAMS  runoff,  sediment 
and  nutrient  algorithms,  with  some  simplifications.  There  are  two 
major  changes  in  model  design:  AGNPS  is  an  event  model,  and  it 
simulates  yield  from  multiple  source  areas.  A  grid  is  used  with  each 
cell  assumed  to  have  uniform  conditions.  Flow  and  sediment  are  routed 
in  overland  and  channel  flow  to  the  watershed  outlet  for  a  storm  event. 
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The  ANSWERS  model  (Area!,  Nonpoint  Source  Watershed  Environmental 
Response  Simulation)  predicts  runoff,  along  with  erosion,  sediment 
transport  and  deposition,  from  field  or  basin  scale  watersheds  on  a 
storm  event  basis  (Beasley  and  Huggins,  1982).  The  watershed  is 
overlayed  with  a  square  grid  (typically  1-4  hectares  per  cell)  with 
physically-based  parameter  values  specified  for  each  cell.  The  ability 
to  define  this  level  of  spatial  variability  in  the  watershed  results  in 
a  more  accurate  representation  of  the  watershed  and  enables  evaluation 
of  changes  in  land  use  and  management  practices  at  the  cell  level.  In- 
cluded is  the  ability  to  simulate  various  structural  practices  such  as 
diversions  and  field  detention  basins.  More  recent  versions  of  the 
model  include  simple  nutrient  yield  predictions,  and  a  version  includ- 
ing a  detailed  phosphorus  model  has  been  planned  (Beasley  et  al.,  1982). 

Several  water  quality  models  have  been  developed  around  the  hydro- 
logic  framework  of  the  Stanford  Watershed  Model  (SWM).  The  history  of 
the  development  of  this  "family"  of  models  is  detailed  by  Barnwell  and 
Johanson  (1981).  Since  an  adaptation  of  the  SWM  serves  as  the  hydro- 
logic  base,  these  models  all  follow  the  lumped  parameter  concept  and 
must  be  calibrated  for  each  site. 

The  first  use  of  the  SWM  for  water  quality  was  the  Pesticide 
Transport  and  Runoff  (PTR)  model  (Crawford  and  Donigian,  1973).  This 
was  expanded  into  the  Agricultural  Runoff  Management  (ARM)  model 
(Donigian  and  Crawford,  1976a),  which  was  further  improved  by  Donigian 
et  al.  (1977).  ARM  predicts  runoff,  sediment  loss,  nutrient  uptake  by 
plants,  and  soil/water  interactions  of  nutrients  and  pesticides.  The 
Non-point  Source  (NPS)  model  (Donigian  and  Crawford,  1976b)  was  devel- 
oped as  a  simpler  version  of  ARM  and  included  algorithms  to  make  it 
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applicable  to  urban  modeling.  Both  ARM  and  NPS  are  continuous  simu- 
lation models;  neither  have  channel  routing  routines,  thus  provide  only 
edge-of-field  loadings. 

The  Hydrologic  Simulation  Program— Fortran  (HSPF)  (Johanson  et 
al.,  1980)  is  a  culmination  of  these  earlier  modeling  efforts,  combin- 
ing components  of  ARM  and  NPS  in  a  model  framework  that  includes  a  data 
manager  (Time  Series  Management  System),  channel  routing  and  instream 
chemical  processes.  Improved  instream  processes  and  statistical  anal- 
ysis for  risk  assessment  have  been  included  in  recent  versions  of  the 
model  (Johanson  et  al.,  1984;  Donigian  et  al.,  1983).  Although  it  is  a 
lumped  parameter  model  in  design,  HSPF  can  account  for  some  watershed 
variability  through  its  ability  to  route  flow  from  different  sub- 
watersheds  delineated  within  the  watershed  of  interest.  HSPF  has  been 
proposed  for  use  in  agricultural  nonpoint  source  modeling  (Woodruff  et 
al.,  1982;  Donigian  et  al.,  1983),  and  BMP  evaluation  (Bicknell  et  al., 
1984),  however,  as  a  lumped  parameter  model,  calibration  is  required, 
and  the  availability  of  suitable  data  will  often  be  a  limiting  factor. 

Of  the  current  water  quality  models,  none  adequately  assess  the 
effectiveness  of  BMPs  on  both  the  field  and  watershed  level.  CREAMS  is 
specifically  designed  for  evaluating  alternative  practices,  but  is 
strictly  a  field  scale  model.  HSPF  predicts  runoff  and  water  quality 
constituents  at  the  basin  outlet,  but  provides  no  information  at  inter- 
ior points  in  the  basin.  Networking  field  sized  areas  together  would 
be  possible,  but  this  approach  is  neither  practical  nor  is  it  consis- 
tent with  the  model's  design  and  intended  application.  ANSWERS  and  the 
AGNPS  model  do  provide  the  desired  modeling  framework,  predicting  both 
field  (cell)  and  basin  scale  results.  A  major  limitation  of  both  is 
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that  they  are  event  models,  and  for  water  quality  modeling,  long  term 
simulation  is  desirable  (DelVecchio  and  Knisel,  1982;  Frere  et  al., 
1980;  DeCoursey,  1985) . 


CHAPTER  III 
THE  CREAMS-WT  FIELD  SCALE  MODEL:  HYDROLOGY  COMPONENT 


Development 
Description  of  the  CREAMS  Model 

Background.  CREAMS  was  developed  as  a  tool  for  evaluating  the 
relative  effectiveness  of  alternative  management  practices  in  reducing 
nonpoint  pollution  (Knisel,  1980).  Objectives  in  the  design  of  CREAMS 
were  that  the  model  should  1)  simulate  major  physical  processes  that 
control  water  balance,  erosion,  sediment  yield,  and  movement  of  plant 
nutrients  and  pesticides,  2)  use  physically  based  parameters  that  can 
reflect  changes  in  management  systems,  3)  operate  on  a  daily  time  step 
for  computational  efficiency,  and  4)  be  field  scale,  since  this  is  the 
common  base  for  selection  of  alternative  management  practices 
(DelVecchio  and  Knisel,  1982).  Model  documentation  includes  an 
extensive  users  manual  (Knisel,  1980)  and  there  are  several  good 
summaries  of  the  model  in  the  literature  (Knisel  et  al . ,  1983;  Knisel 
and  Foster,  1980) . 

CREAMS  is  composed  of  hydrology,  erosion  and  chemistry  sub-models. 
The  three  sub-models  are  run  independently  (in  sequence),  with  "pass" 
files  used  to  link  the  models  providing  the  necessary  data  for  the 
next  sub-model.  The  model  was  designed  this  way  so  that  the  sub-models 
could  be  used  separately  if  desired  and  so  the  model  could  be  used  on 
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smaller  computers.  The  obvious  disadvantage  is  that  there  is  no  feed- 
back between  the  modules.  This  structure  also  made  the  model  easier  to 
modify,  as  each  module  can  be  handled  separately. 

Hydrology  algorithms.  The  hydrology  sub-model  uses  rainfall, 
temperature,  average  monthly  radiation,  soils  and  land  use  parameters 
to  predict  runoff,  evapotranspiration,  percolation  and  soil  water 
content.  Using  daily  rainfall  records,  a  modified  SCS  curve  number 
method  is  used  to  predict  runoff.  A  second  hydrology  option  is 
available  which  uses  hourly  or  breakpoint  rainfall  with  an  infiltration 
equation.  The  daily  rainfall  option  was  selected  as  the  base  for 
CREAMS-WT  because  1)  breakpoint  rainfall  data  are  generally  not 
available,  2)  an  infiltration  model  does  not  conceptually  represent  the 
rainfall /runoff  process  on  flatwoods  watersheds,  and  3)  the  analysis  by 
Capece  (1984)  indicated  that  the  curve  number  method  has  good  potential 
for  predicting  runoff  from  flatwoods  watersheds. 

The  soil  profile  is  divided  into  seven  layers  with  accounting  of 
the  soil  moisture  within  each  layer.  Percolation  is  computed  by  moving 
soil  water  down  through  the  soil  layers,  to  be  lost  out  the  bottom  as 
deep  seepage.  Evapotranspiration  (ET)  is  calculated  using  an  adapta- 
tion of  the  Penman  equation.  Leaf-area-index,  radiation  and  available 
soil  moisture  determine  actual  ET  from  potential  ET.  A  simple  plant 
root  growth  model  is  used  to  partition  the  ET  extraction  of  soil  water 
among  the  different  soil  layers. 

The  daily  rainfall  option  of  CREAMS  uses  the  SCS  runoff  equation 
(equation  2-1)  for  predicting  runoff.  The  storage  parameter,  S,  is 
calculated  as, 
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(ULE  -  SM) 

Smx [3-1] 

ULE 


where  ULE  is  the  upper  limit  of  soil  water  storage  in  the  profile,  SM 
is  the  soil  water  content,  and  Smx  is  the  maximum  value  of  S.  This 
maximum  value  is  determined  from  the  SCS  equation  relating  curve  number 
and  S  as, 
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Smx  =  10  [3-2] 

CNj 


where  CNy  is  the  curve  number  for  antecedent  moisture  condition  I.  The 
curve  number  (condition  II)  is  an  input  parameter  and  is  converted  by 
the  model  to  the  equivalent  CNj  value.  The  relationship  between  CNjj 
and  Smx  is  described  in  greater  detail  in  Appendix  A. 

As  calculated  in  equation  3-1,  all  layers  of  the  profile  are 
equally  important  in  determining  the  soil  moisture,  hence  storage. 
Runoff  is  normally  expected  to  be  more  sensitive  to  the  moisture 
content  in  the  surface  layers  of  the  soil.  To  better  represent  this, 
CREAMS  uses  a  weighted  form  of  equation  3-1.  The  surface  layer  is 
weighted  heaviest  with  weighting  factors  decreasing  with  depth. 

Limitations  for  flatwoods  watersheds.  The  modified  SCS  Curve 
Number  method  used  by  CREAMS  for  predicting  runoff  volumes  was  one  of 
the  methods  evaluated  by  Capece  (1984).  This  algorithm  was  compared 
with  several  other  predictive  methods  and  found  to  be  one  of  the  best 
ones  tested.  However,  good  estimates  of  runoff  volume  were  found  to  be 
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dependent  on  knowing  the  available  storage  in  the  soil.  The  CREAMS 
model  uses  an  internal  moisture  accounting  procedure  to  determine  S, 
whereas  Capece  (1984)  used  the  ARS  storage  curve. 

For  use  in  flatwoods  watersheds,  a  significant  problem  with  the 
CREAMS  model  is  its  inability  to  account  for  a  fluctuating  water  table 
in  the  soil  profile.  There  are  two  associated  effects.  First,  the 
high  saturated  conductivity  of  the  sandy  soil  results  in  the  model 
predicting  a  significant  amount  of  deep  seepage,  while  in  reality,  deep 
percolation  in  these  basins  is  negligible.  Second,  available  storage 
in  the  soil  is  overestimated  when  the  water  table  is  in  the  profile, 
resulting  in  underprediction  of  runoff.  Thus,  the  effect  on  runoff 
prediction  is  especially  significant  since  runoff  occurs  only  when  the 
water  table  is  near  the  surface,  and  this  is  the  condition  where  the 
model's  estimate  of  available  storage  is  least  accurate. 

Bengtson  and  Carter  (1983)  applied  CREAMS  to  a  Mississippi  water- 
shed that  is  also  very  flat  and  has  a  seasonally  high  water  table. 
They  found  that  CREAMS  underestimated  runoff  during  wet  periods  and 
overestimated  runoff  during  dry  periods.  A  similar  result  is  reported 
by  Pathak  et  al.  (1982),  who  simulated  native  rangeland  runoff  in 
Oklahoma. 

Concepts  of  CREAMS-WT 

For  application  to  flat,  high  water  table  watersheds,  CREAMS  must 
be  modified  to  limit  the  amount  of  deep  seepage  that  is  predicted. 
Without  modification,  the  CREAMS  estimate  of  the  available  storage  in 
the  profile  is  not  accurate.  This  error  is  then  reflected  in  the 
runoff  volume  predicted  by  the  curve  number  method. 
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The  CREAMS-WT  model  corrects  these  limitations.  A  high  water 
table  is  modeled  with  limited  or  no  deep  seepage  and  a  modified  method 
is  used  to  calculate  the  available  storage.  The  result  is  more 
accurate  estimation  of  the  actual  storage  available,  thus  better 
prediction  of  runoff  volumes. 

The  following  objectives  guided  the  development  of  the  CREAMS-WT 
model . 

1)  The  model  should  represent  the  physical  processes  occurring  in 
flatwoods  watersheds:  a  high  water  table,  and  limited  deep 
seepage. 

2)  The  method  of  estimating  available  storage,  a  crucial 
parameter  in  the  curve  number  method,  must  be  improved. 

3)  The  general  modeling  objectives  for  CREAMS  should  be  followed: 

-  use  physically-based  parameters,  and 

-  be  computationally  efficient. 

4)  The  current  operation  of  CREAMS  should  not  be  disrupted.  This 
modification  can  then  be  incorporated  as  a  third  hydrology 
option  if  desired.  It  will  also  ease  incorporation  of  this 
modification  in  future  versions  of  CREAMS. 

Description  of  Algorithms 
Water  table  simulation 

CREAMS  considers  only  the  defined  root  zone  depth  (input  parameter 
RD)  in  its  soil  water  accounting.  Additional  routines  were  thus  needed 
in  two  areas.  First,  to  follow  the  water  table  in  the  root  zone,  and 
second,  to  follow  the  water  table  when  it  falls  below  the  root  zone 
into  the  "lower  zone". 
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Water  content  above  field  capacity  is  "routed"  by  CREAMS,  layer  by 
layer,  down  through  the  root  zone.  For  sandy  soil,  the  high  saturated 
hydraulic  conductivity  can  easily  result  in  no  effective  detention  of 
downward  water  movement.  Thus,  on  a  daily  time  step  and  given  adequate 
rainfall,  the  water  content  of  the  profile  is  simply  restored  to  field 
capacity  and  any  excess  becomes  deep  seepage.  The  way  this  deep  seep- 
age is  handled  in  CREAMS-WT  depends  on  which  zone  the  water  table  is  in. 

Water  table  in  the  root  zone.  CREAMS-WT  retains  all  the  CREAMS 
routines  for  calculating  the  moisture  balance  in  the  root  zone.  Added 
routines  take  the  "deep  seepage"  predicted  by  the  CREAMS,  and  beginning 
at  the  bottom  and  working  upwards,  fill  to  saturation  successive  layers 
in  the  root  zone.  An  estimate  of  the  position  of  the  water  table  is 
then  made  by  interpolating  to  find  the  depth  of  saturation  in  the  layer 
just  above  the  topmost  layer  that  is  fully  saturated. 

A  new  parameter,  DSP,  specifies  a  rate  (in  inches/day)  of  subsur- 
face water  movement  out  of  the  modeled  watershed  system.  This  may  be 
used  to  simulate  a  limited  degree  of  deep  percolation,  or  lateral 
seepage.  Deep  percolation  occurs  at  the  rate  specified  by  DSP  while 
the  water  table  is  in  the  root  zone.  When  the  water  table  falls  below 
the  root  zone,  no  deep  percolation  occurs.  One  possibility  is  to  use 
this  component  to  model  baseflow  on  larger  watersheds.  In  this  case, 
the  predicted  "percolation"  should  be  added  to  the  volume  of  runoff 
predicted. 

When  the  water  content  in  the  root  zone  is  at  or  below  field 
capacity  the  water  table  is  assumed  to  be  at  the  boundary  of  the  two 
zones.  On  the  following  day  (next  time  step),  provided  there  is  no 


54 

rainfall  to  raise  the  water  content  in  the  root  zone  above  field 
capacity,  the  water  table  is  considered  to  be  in  the  lower  zone. 

Water  table  in  the  lower  zone.  Water  movement  into  and  out  of  the 
lower  zone  is  tracked  with  the  aid  of  a  water  table  recession  curve. 
Speir  et  al.  (1969)  used  well  data  from  the  Taylor  Creek  watershed  to 
develop  a  recession  curve  showing  water  table  depth  versus  number  of 
rainless  days.  With  an  estimate  of  the  drainable  porosity  of  the  lower 
zone,  a  storage  versus  time  relationship  can  be  derived  from  the 
recession  curve. 

The  function  used  to  represent  the  recession  curve  is, 


D  =  Tk  [3-3] 


where  D  is  depth  of  the  water  table  from  the  surface,  in  feet,  T  is 
time  in  days,  and  k  is  a  coefficient.  A  value  of  k  =  0.33  closely 
represents  the  curve  presented  by  Speir  et  al.  (1969).  The  recession 
curve  begins  at  the  surface,  but  the  only  portion  of  the  curve  used  in 
the  model  is  that  which  lies  below  the  root  zone  depth  (RD). 

Introducing  a  parameter  for  the  drainable  porosity,  PORS  (in/ft), 
depth  in  the  lower  zone  can  be  converted  to  storage,  LZS  (in  inches), 

LZS  =  PORS  *  D  [3-4] 

which  combined  with  equation  [3-3]  becomes, 


LZS  =  PORS  *  Tk  [3-5] 
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This  storage  versus  time  relationship  is  used  to  track  the  water  table 
and  provide  a  water  balance  in  the  lower  zone.  In  the  following 
discussion,  the  subscripts  '1-1'  and  'i'  refer,  respectively,  to  the 
previous  and  current  time  step  (day). 

When  the  water  table  is  identified  as  being  in  the  lower  zone, 
rainfall  first  fills  the  root  zone  to  field  capacity  before  any  excess 
moves  into  the  lower  zone  (variable  SEP).  For  days  without  water  re- 
charging the  lower  zone,  the  time  location  on  the  curve  is  increased  by 
one  day,  with  one  exception.  For  days  with  rainfall  greater  than  0.1 
inch,  no  change  is  made  in  the  water  table  since  ET  is  assumed  to  be 
satisfied  from  the  root  zone.  For  days  with  seepage  into  the  lower 
zone,  the  position  on  the  curve  is  moved  up  representative  of  that 
volume  of  seepage.  These  relationships  can  be  expressed  in  the 
following  way. 

When  SEP  =  0  and  precipitation  <  0.1  inch, 

T.   =  Ti_1  +  1  [3-6] 


DPET  =  PORS  *  (  T.k  -  T.  ,k  )  [3-7] 


where  DPET  is  the  volume  of  water  lost  from  the  lower  zone  in  that  time 
step.  These  withdrawals  from  the  lower  zone  are  summed  and  reported  as 
"Lower  Zone  ET". 

When  SEP  >  0,  the  new  storage  is  found  as, 


LZSi  =  PORS  *  Ti_1  -  SEP  [3-8] 
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and  the  corresponding  reference  time  on  the  recession  curve  is 

Ti  =  (  LZS^PORS  )**(l/k)  [3-9] 

The  total  volume  of  seepage  into  the  lower  zone  is  reported  as  "Lower 
Zone  Recharge".  For  any  period  of  time  where  the  water  table  starts  in 
the  root  zone,  moves  to  the  lower  zone,  and  ends  up  back  in  the  root 
zone,  the  Lower  Zone  Recharge  will  equal  the  Lower  Zone  ET. 

When  the  volume  of  seepage  (SEP)  from  the  root  zone  is  greater 
than  the  available  storage  in  the  lower  zone  (between  the  current  water 
table  and  the  bottom  of  the  root  zone),  the  excess  is  transferred  back 
to  the  root  zone.  The  soil  water  is  then  modeled  with  the  root  zone 
component  in  the  following  time  step. 
Runoff  prediction 

The  CREAMS  curve  number  method  requires  an  estimate  of  soil 
moisture  (SM  in  equation  3-1)  which  is  translated  to  an  effective  curve 
number.  CREAMS  uses  a  weighting  of  the  actual  storages  in  the  seven 
layers  of  the  root  zone.  The  top  layer  is  weighted  most  heavily  and 
the  weighting  decreases  moving  down  the  profile.  This  is  an  attempt  to 
more  closely  represent  the  effects  of  soil  moisture  on  infiltration 
but  is  not  desirable  for  the  storage  based  system  being  modeled.  This 
weighting  was  removed  so  that  for  CREAMS-WT  the  actual  available 
storage  in  the  root  zone,  rather  than  a  weighted  value,  is  used  in 
predicting  runoff. 

Runoff  is  the  first  of  the  various  components  of  the  daily  water 
balance  that  is  computed.  The  difference  between  precipitation  and 
runoff  is  then  assumed  to  infiltrate.  To  ensure  that  this  volume  does 
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not  exceed  the  available  storage  in  the  profile,  the  "strict  storage" 
prediction  of  runoff  as  precipitation  minus  available  storage  is  cal- 
culated, and  the  larger  of  these  two  runoff  predictions  is  used.  The 
curve  number  prediction  will  generally  be  higher  but  comparing  against 
the  strict  storage  prediction  ensures  that  the  available  storage  in  the 
profile  will  not  be  exceeded. 

Calibration  and  Testing 

Calibration.  Few  parameters  were  adjusted  during  the  calibration 
of  the  model.  A  minor  change  was  made  in  the  exponent  of  the  water 
table  recession  curve  ('k'  in  equation  3-3,  equivalent  to  model 
parameter  RXP).  This  was  changed  from  0.33  to  a  value  of  0.36  which 
better  matched  recession  curves  from  Bass  West  well  data.  This  also 
matches  closely  with  a  value  of  0.355  determined  by  Capece  (1984). 

The  most  important  calibration  was  of  the  annual  predicted  ET.  As 
will  be  discussed  later,  the  ET  model  component  is  very  sensitive  to 
parameter  changes.  This  is  especially  critical  since  ET  is  a  large 
component  in  the  annual  water  balance.  Leaf-area-index  (LAI)  is  the 
primary  parameter  affecting  ET  that  can  be  reasonably  varied.  (For 
instance,  it  would  not  be  reasonable  to  vary  temperature  from  recorded 
values.)  While  LAI  is  a  physically  based  parameter,  it  is  difficult  to 
determine  appropriate  values.  Data  from  Jones  et  al.  (1984c)  were  used 
to  determine  initial  estimates,  and  these  values  were  then  adjusted  so 
the  model  output  matched  annual  ET  estimates  better.  Additional  com- 
ments on  estimating  LAI  appear  in  the  following  sections  on  sensitivity 
analysis  and  parameter  estimation. 
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With  additional  calibration  it  would  be  possible  to  match  the 
observed  watershed  response  more  closely.  While  this  would  improve  the 
perception  of  the  model's  performance,  additional  calibration  was  not 
done  for  the  following  reasons.  An  objective  for  the  model  is  that  it 
be  applicable  for  predicting  the  hydrologic  and  water  quality  response 
of  ungaged  watersheds.  Thus,  the  ability  to  achieve  reasonable  simula- 
tions without  extensive  calibration  demonstrates  the  model's  capacity 
for  use  in  ungaged,  uncalibrated  application.  In  addition,  uncertain- 
ties in  parameter  values  as  well  as  in  the  observed  watershed  response 
imply  some  expected  variability  between  predicted  and  observed  values. 

Bass  West  simulation.  Data  from  the  Bass  West  site  of  the  Upland 
Detention/Retention  Demonstration  Project  (Goldstein,  1986)  were  used 
during  model  development,  for  initial  evaluation  of  the  model  and  for 
calibration  of  some  parameters.  Comparison  with  runoff  volumes  was  not 
possible  for  this  site,  as  precipitation  and  water  table  elevation  were 
the  only  data  judged  to  be  suitable  for  use  in  model  evaluation. 

Parameters  for  CREAMS-WT  are  physically  based  and  were  estimated 
from  various  sources:  Goldstein  (1986),  the  CREAMS  manual,  Soil  Survey 
reports,  Capece  (1984)  and  Speir  et  al.  (1969).  The  simulation  of  the 
Bass  West  pasture  was  run  from  June  1979  to  December  1982,  the  length 
of  the  rainfall  record,  although  water  table  observations  did  not  begin 
until  1980.  This  allowed  a  period  for  model  equilibration  which  mini- 
mized the  effects  of  the  initial  parameter  values  chosen.  The  para- 
meter file  is  listed  in  Appendix  B. 

Accurate  runoff  prediction  is  the  goal  of  the  model.  However,  the 
importance  of  the  soil  moisture  content  in  the  runoff  model  would  sug- 
gest that  the  accuracy  of  the  predicted  soil  water  content  is  crucial 
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to  the  accurate  representation  of  the  hydrologic  system  and  to  runoff 
prediction.  The  predicted  water  table  was  thus  used  as  one  measure  of 
model  evaluation.  Several  comments  on  the  observed  water  table  data 
are  appropriate.  First,  recording  the  water  table  elevation  is  a 
simple  and  reliable  procedure,  thus  it  provides  good  data  for  use  in 
model  evaluation.  However,  one  well  in  a  20  to  3600  acre  watershed  is 
clearly  limited  in  how  well  it  represents  the  average  storage  over  the 
entire  watershed.  This  must  be  considered  in  comparing  water  tables 
since  the  CREAMS-WT  predicted  water  table  is  a  value  representative  of 
the  watershed  and  thus  may  at  times  reasonably  differ  from  a  water 
table  elevation  measured  at  one  particular  location. 

The  simulated  and  observed  water  tables  are  compared  in  Figure 
3-1.  The  general  trends  in  the  water  table  are  followed  fairly  well. 
The  high  frequency  'spikes'  in  the  observed  water  table  are  not 
followed  by  the  model  due  to  two  factors:  1)  the  model  operates  on  a 
daily  time  step,  and  2)  the  runoff  volume  predicted  by  the  curve  number 
method  is  subtracted  from  the  water  balance  the  same  day,  whereas  in 
reality,  the  runoff  hydrograph  will  extend  over  several  days.  One 
consistent  discrepancy  is  that  the  simulated  water  table  is  too  high  in 
the  winter  months.  This  is  due  to  inaccuracies  in  estimating  ET, 
resulting  from  the  lack  of  specific  site  radiation,  temperature,  and 
leaf  area  index  data.  The  low  water  table  predicted  in  1983  may  in 
part  be  due  to  uncertainties  in  estimating  ET  parameters,  and  may  also 
indicate  that  values  for  the  parameters  PORSUB  and  FUL  are  too  high. 
However,  the  primary  cause  is  attributed  to  rainfall  events  missing  in 
the  recorded  data. 


60 


WATER  TABLE  DEPTH,  feet 


61 

Evaluation  of  the  CREAMS-WT  Runoff  Equation 

One  aspect  of  model  evaluation  is  to  determine  if  the  algorithms 
respond  in  a  way  that  is  reasonable,  or  expected,  of  the  system  being 
modeled.  The  following  discussion  will  explore  the  characteristics  of 
the  rainfall /runoff  relationship  used  in  CREAMS-WT.  As  mentioned  in 
the  literature  review  section,  the  appropriateness  of  the  curve  number 
method  has  come  into  question  for  a  wide  range  of  reasons.  Since  the 
method  was  developed  from  data  on  soils  and  watersheds  that  are  more 
typical  of  the  rest  of  the  country,  it  is  even  more  important  to 
address  the  issue  of  the  theoretical  and  conceptual  appropriateness  of 
the  method  for  simulating  flatwoods  watersheds.  (Note  that  Capece 
[1984]  addresses  the  pragmatic  question  of  "does  it  work"?) 

Many  of  the  current  hydrologic  models  are  based  on  the  concept  of 
limited  infiltration  producing  runoff.  The  curve  number  method  of 
estimating  runoff  is  based  on  a  storage  concept  rather  than  on  infil- 
tration. It  does  not  consider  rainfall  intensity  and  duration  —  crucial 
factors  in  the  generation  of  runoff  in  most  of  the  country.  However, 
for  modeling  a  physical  system  that  is  storage  based,  as  are  the 
flatwoods  watersheds,  the  curve  number  method  is  a  good  model. 

The  SCS  runoff  equation  (equation  2-1)  can  be  arranged  in  the 
form, 


Q    [  1  -  0.2S/P  ]2 

-  =  - —  -  [3-10] 

P    [  1  +  0.8S/P  ] 
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By  defining  available  storage  in  the  profile  as 

Sav  =  ULE  -  SM  [3-11] 

equation  3-1  can  be  written 

S  =  Sav  *  Smx/ULE  [3-12] 

Substituting  this  in  equation  3-10  gives  the  effective  relationship  which 
is  used  to  calculate  runoff  in  CREAMS-WT, 


q         [1-0.2   (Smx/ULE)    *   (Sav/P)    ]2 

_     = [3-13] 

P         [1+0.8   (Smx/ULE)   *   (Sav/P)    ] 


For  ease  in  graphical  presentation,  consider  unit  precipitation 
which  simplifies  the  expression  to, 


[1-0.2   (Smx/ULE)    *  Sav  ]2 

Q     = - [3-14] 

[1+0.8   (Smx/ULE)    *  Sav   ] 


A  graph  of  Q  versus  Sav  is  shown  in  Figure  3-2  for  several  values  of 
Smx/ULE.  The  factor  Smx/ULE  is  fixed  for  a  given  simulation  and 
represents  the  ratio  of  the  maximum  profile  storage  as  determined  by 
the  curve  number  to  the  maximum  storage  as  specified  by  the  input 
parameter  UL.  When  this  ratio  is  one,  equation  3-13  reduces  to  a 
formula  analogous  to  the  SCS  equation  (3-11).  Also  shown  in  Figure  3-2 
is  the  straight  line  represented  by  the  strict  storage  relationship, 
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Q  =  P  -  Sav  [3-15] 

which  for  unit  precipitation  becomes, 

Q  =  1  -  Sav  [3-16] 

While  the  flatwoods  watersheds  respond  as  a  storage  based  system, 
the  strict  storage  relationship  is  less  than  ideal.  Lower  areas  of  a 
watershed  (such  as  sloughs)  begin  contributing  to  runoff  before  the 
entire  watershed  is  saturated.  The  curve  number  method  (the  Smx/ULE=1 
curve  in  Figure  3-2)  is  more  representative  of  this  process.  Runoff 
begins  before  available  storage  is  completely  filled,  and  predicted 
runoff  is  always  greater  than  the  strict  storage  relationship  would 
indicate.  For  the  original  curve  number  method  (Smx/ULE  =  1),  runoff 
is  predicted  out  to  the  point  where  storage  is  five  times  the  amount  of 
precipitation. 

The  CREAMS-WT  runoff  algorithm  has  an  additional  variable,  the 
Smx/ULE  ratio.  With  a  higher  curve  number  (inversely  related  to  Smx) 
this  ratio  is  less  than  unity  and  increased  runoff  is  predicted  for  a 
given  value  of  Sav.  When  the  ratio  is  greater  than  unity,  decreased 
runoff  is  predicted.  When  the  Smx/ULE  ratio  is  greater  than  one,  the 
curve  crosses  the  straight  line  representing  the  strict  storage 
relationship.  As  discussed  previously  in  the  "Runoff  Prediction" 
section,  the  model  forces  the  runoff  prediction  to  satisfy  the  strict 
storage  relationship.  Thus  when  the  Smx/ULE  ratio  is  greater  than  one, 
the  runoff  prediction  follows  the  curve  (moving  from  right  to  left) 
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1.8 


AVAILABLE  STORAGE  (Sav),  ins 


Figure  3-2.  Runoff  predicted  by:  the  CREAMS-WT  runoff  equation 

(graphed  for  several  values  of  Smx/ULE),  and  the  linear, 
strict  storage  relationship  of  Q=P-Sav.  These  functions 
are  graphed  considering  unit  precipitation  (P=l). 
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until  it  intersects  the  straight  line,  then  it  follows  the  straight 
line  relationship  to  zero  storage. 

While  the  runoff  prediction  must  conform,  at  minimum,  to  the 
strict  storage  relationship,  the  length  of  the  "tail"  to  the  right, 
which  represents  increased  runoff,  can  be  controlled  through  the 
Smx/ULE  ratio.  This  feature  may  be  useful  in  representing  morphologic 
variability  between  watersheds  as  well  as  the  effects  of  some  manage- 
ment practices.  Thus,  the  implementation  of  the  curve  number  method  as 
used  in  CREAMS-WT  appears  to  be  particularly  appropriate  for  predicting 
runoff  from  flat,  high  water  table  watersheds.  In  a  simple  manner,  it 
represents  the  variable  source  area  concept  and  is  descriptive  of  the 
rainfall/runoff  process  on  flatwoods  watersheds.  Further  analysis  is 
needed  to  determine  the  relationship  between  the  various  parameters  and 
watershed  characteristics.  Unfortunately,  analysis  is  limited  at  this 
point  by  a  lack  of  applicable  data. 

Verification 

Two  sources  of  data  were  used  in  verification  of  the  hydrology 
model.  Small  watershed  data  from  the  UD/RDP  were  used  for  specific 
comparison  of  model  predictions  with  observed  data,  and  long  term  water 
balance  studies  on  the  Taylor  Creek  watershed  were  used  to  test  how 
reasonable  and  representative  the  model  predictions  are  over  a  long 
simulation. 
Armstrong  Slough 

The  Armstrong  Slough  site  was  simulated  without  calibration.  Both 
CREAMS-WT  and  the  original  CREAMS  models  were  run  using  this  same 
parameter  set  (included  in  Appendix  B).  Parameters  were  taken  largely 
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from  the  Bass  West  pasture  simulation.  A  small  amount  of  percolation 
was  modeled,  in  part  as  a  demonstration  of  this  aspect  of  the  model. 
Since  there  is  little  specific  data  related  to  this  pathway,  its 
inclusion  here  is  merely  conjecture.  In  a  larger  watershed  there  may 
be  some  losses  due  to  deep  seepage  or  lateral  outflow,  or  this  could  be 
used  to  represent  baseflow. 

The  simulated  and  observed  water  tables  are  shown  in  Figure  3-3. 
The  inconsistencies  in  the  water  table  prediction  are  due  in  part  to 
inaccuracies  in  simulating  ET  as  mentioned  previously.  However,  it 
also  appears  that  most  of  the  discrepancy  may  be  due  to  errors  in  the 
rainfall  data.  In  June  and  July  of  1980,  there  are  two  significant 
rises  in  the  observed  water  table  which  are  not  echoed  in  the  simulated 
water  table.  This  is  typical  of  the  type  of  discrepancy  that  would 
occur  when  a  rainfall  event  on  the  watershed  does  not  appear  in  the 
recorded  rainfall  data.  Errors  in  the  simulated  water  table  are 
perpetuated  until  the  water  table  is  "reset"  by  reaching  the  surface, 
which  in  this  case  is  more  than  a  year  later.  It  is  important  to 
recognize  these  limitations  in  the  observed  data  when  making 
comparisons  with  model  results. 

The  components  of  the  annual  water  balance  as  observed  and  simula- 
ted by  the  two  models  are  compared  in  Table  3-1.  The  CREAMS  simulation 
was  included  to  demonstrate  the  significance  of  the  simulated  water 
table.  The  excessive  percolation  predicted  by  the  unmodified  CREAMS 
version  clearly  biases  the  other  components  of  the  water  balance. 
Specifically,  runoff  and  ET  are  both  underpredicted. 

In  evaluating  the  performance  of  CREAMS-WT,  the  following  points 
must  be  kept  in  mind.  The  error  in  the  simulated  water  table  due  to 
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missing  rainfall  data  resulted  in  a  low  water  table  through  most  of 
1981.  This  would  produce  high  amounts  of  "lower  zone  ET"  thus 
inflating  the  ET  estimate,  and  would  also  under-predict  runoff  because 
of  the  excess  storage  in  the  profile.  Thus,  the  effects  of  this  error 
carry  through  to  the  water  balance.  Water  extracted  from  profile 
storage  in  1980  inflates  ET  for  that  year.  This  storage  is  refilled  in 
1981,  with  both  ET  and  runoff  underpredicted  because  of  the  very  low 
water  table.  For  1982,  the  accuracy  of  the  observed  runoff  data  must 
be  questioned,  as  "observed"  ET  is  unreasonably  low  for  a  year  with 
above  average  rainfall. 

Considering  the  above  comments,  CREAMS-WT  models  the  site  very 
well.  Runoff  and  other  components  of  the  water  balance  are  all  pre- 
dicted reasonably  well.  And  again,  this  simulation  was  not  calibrated. 

SEZ  Dairy 

The  SEZ  Dairy  watershed  was  simulated  without  calibration.  The 
water  table  and  annual  water  balance  are  compared  with  observed  data  in 
Figure  3-4,  and  Table  3-2.  The  accuracy  of  the  observed  1982  runoff 
must  again  be  questioned  based  on  the  unreasonably  low  ET  values.  The 
above  comments  on  the  accuracy  of  the  water  table  simulation  also  apply 
here. 

Long  Term  Simulation 

Two  to  three  years  of  data  are  minimal  for  good  testing  of  a 
model.  Rainfall  varies  greatly  from  year  to  year,  and  a  short  period 
of  record  can  include  only  a  limited  sampling  of  the  range  of 
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conditions  that  may  occur.  It  is  desirable  to  test  a  model  with  the 
full  range  of  input  conditions  under  which  it  may  be  used.  The  Upland 
Detention/Retention  Demonstration  Project  unfortunately  occurred  during 
a  period  of  record  drought.  While  this  may  check  the  response  of  the 
model  to  low  rainfall  conditions,  this  is  generally  the  condition  of 
least  concern.  Runoff  prediction  for  average  to  high  flow  conditions 
are  usually  of  greater  interest. 

Hydrologic  data  have  been  collected  on  the  Upper  Taylor  Creek 
watershed  for  over  fifteen  years.  The  CREAMS  model  was  designed  to 
simulate  field  sized  areas,  thus  cannot  be  used  directly  to  simulate 
this  watershed.  However,  long  term  studies  of  this  watershed  have 
yielded  some  consistent  relationships  between  rainfall,  runoff,  and 
evapotranspiration,  the  major  components  of  the  water  balance  for  this 
area  (Speir  et  al . ,  1969;  Allen  et  al.,  1982a). 

As  a  means  of  testing  the  long  term  response  of  the  model,  the 
parameters  for  the  Bass  West  site  were  used  with  twenty  years  of 
rainfall  records  from  this  area  (Okeechobee  Hurricane  Gate-6,  1960- 
1979).  Actual  rainfall  records  from  the  region  were  used  so  that 
natural  variability  in  precipitation  patterns  and  amounts  would  be 
represented.  Annual  rainfall  amounts  over  the  twenty  year  period  of 
record  used  for  simulation  have  a  mean  of  46.8  inches  and  standard 
deviation  of  8.20.  Over  the  more  than  50  year  period  of  record  at  that 
station,  similar  statistics  of  47.59  for  the  mean,  and  8.94  for  the 
standard  deviation  (Knisel  et  al . ,  1985)  show  that  the  twenty  year 
period  used  is  representative  of  the  long  term  record  in  both  amount 
and  annual  distribution.  Thus,  the  range  of  dry  to  wet  years  were 
represented  in  this  time  span. 
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Table  3-3.  Annual  water  balance  summary  for  20-year  CREAMS-WT 
simulation. 


STORAGE 
YEAR    PRECIP    RUNOFF     ET     CHANGE 
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1.303 

39.687 
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0.552 
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4.492 
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41.000 
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35.186 

1.334 

78 
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79 
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-2.296 

Average  46.790    8.545    38.299     2.370 
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The  predicted  water  balance  for  the  20-year  period  is  shown  in 
Table  3-3.  The  results  of  this  long  term  simulation  show  that  the 
CREAMS-WT  model  is  stable  and  that  long  term  predictions  are 
reasonable. 

Keeping  in  mind  the  possible  differences  between  the  hydrologic 
response  of  the  field  and  watershed,  limited  comparison  can  be  made 
with  observed  data.  On  a  long  term  basis,  the  water  balance  for  the 
Taylor  Creek  watershed  is  represented  by, 

Rainfall  =  Runoff  +  Evapotranspiration 

Of  the  other  components  usually  found  in  a  water  balance,  change  in 
storage  is  considered  negligible  over  a  long  period.  Recharge  to  deep 
aquifers  as  deep  percolation  is  negligible  in  this  area  and  can  be 
ignored  (Stewart,  1980).  For  an  annual  rainfall  of  47.9  inches, 
estimates  of  ET  for  Upper  Taylor  Creek  (W-2)  and  Northwest  Taylor  Creek 
(W-3)  are  35.4  and  37.1  inches  respectively  (Speir  et  al . ,  1969).  This 
compares  closely  with  the  predicted  annual  ET  of  38.3  inches,  which  was 
averaged  over  the  20  years  simulated.  A  slightly  higher  ET  from  the 
Bass  Pasture  simulation  is  to  be  expected  for  two  reasons.  First, 
runoff  would  be  expected  to  increase  (and  ET  decrease)  in  moving  from  a 
pasture,  to  a  headwaters  watershed  (W-3),  to  a  larger  watershed  (W-2). 
Second,  higher  ET  amounts  may  result  from  the  conversion  of  unimproved 
pasture  and  range  to  improved  pasture  (Yates  et  al.,  1982). 

This  twenty  year  simulation  provides  some  important  verification 
of  the  operation,  appropriateness,  and  applicability  of  CREAMS-WT.  In 
summary:  the  model  is  stable,  it  performs  predictably  over  the  range 
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of  natural  rainfall  variability,  and  long  term  average  annual  predic- 
tions compare  favorably  with  observed  values. 

Sensitivity  Analysis 

The  twenty  year  simulation  described  above  was  used  as  the  base 
simulation  for  the  parameter  sensitivity  analysis.  Sensitivity  was 
computed  using  the  twenty  year  average  annual  value  of  the  desired 
output  variables.  The  effect  of  a  parameter  on  model  response  may  vary 
greatly  between  a  wet  and  a  dry  year.  Simulation  using  this  period  of 
historical  rainfall  averages  the  parameter  effects  over  the  range  of 
rainfall  conditions,  minimizing  the  effects  of  individual  years. 

Parameters  chosen  for  analysis  were  those  known  to  be  sensitive, 
those  that  may  be  more  difficult  to  determine,  and  the  new  parameters 
introduced  in  CREAMS-WT.  Because  the  magnitudes  and  the  range  of 
variation  of  various  parameters  may  be  greatly  different,  a  standard 
percentage  change  in  parameters  was  not  used.  Rather,  parameters  were 
varied  from  the  base  value  according  to  the  reasonable  range  of  values 
for  that  parameter. 

Relative  sensitivity  is  calculated  as, 


Onew  -  Ob     Pb 
Sens  =   * — -  [3-17] 

Ob     Pnew  -  Pb 


where  Onew  and  Ob  are  the  values  of  the  output  variable  for  the  new  and 
the  base  simulations,  and  Pnew  and  Pb  are  the  corresponding  parameter 
values.  This  method  of  calculating  sensitivity  provides  a  consistent 
measure  for  comparison  between  parameters.  A  value  of  zero  means  that 
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a  change  in  the  parameter  had  no  effect  on  the  output.  Higher  values 
of  the  relative  sensitivity  indicate  increased  sensitivity  of  the  model 
output  to  the  parameter  in  question.  A  value  of  one  means  the  output 
changes  by  the  same  percentage  the  parameter  is  changed. 

The  results  of  the  sensitivity  analysis  are  shown  in  Table  3-4. 
Note  that  sensitivities  cannot  be  compared  between  different  output 
variables,  for  instance  between  runoff  and  ET,  since  the  relative 
sensitivity  is  inversely  related  to  the  base  magnitude  of  the  variable. 
The  actual  values  of  the  variables  have  also  been  included  in  the  table 
for  a  reference.  The  parameters  evaluated  fall  into  three  categories: 
soil  characterization,  runoff,  and  ET.  The  parameters,  RC,  FUL,  POROS, 
PORSUB  are  generally  insensitive  and  will  be  discussed  individually  in 
the  following  section  on  parameter  estimation. 

The  curve  number  (CN)  and  the  total  free  water  storage,  plant 
available  and  drainable,  in  the  profile  (UL)  are  parameters  used  in  the 
runoff  algorithms.  The  following  observations  can  be  made. 

1.  The  curve  number  is  much  more  sensitive  for  higher  values 
(lower  storage)  than  it  is  for  lower  values. 

2.  UL  is  sensitive  for  both  positive  and  negative  changes, 
although  there  is  some  increase  in  sensitivity  with  decreasing 
storage  as  is  the  case  with  CN. 

3.  UL  is  more  sensitive  than  CN.  If  the  curve  number  is  trans- 
lated to  its  equivalent  storage  value,  the  difference  in 
sensitivity  is  magnified.  This  is  reasonable  as  CN  appears 
only  in  the  Smx/ULE  ratio,  whereas  UL  is  a  part  of  that  ratio 
and  is  also  a  key  component  in  the  soil  moisture  accounting. 
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Table  3-4 

.  Sensi 

tivity 

analysis 

of  selec 

ted  CREAMS-WT  parameters 

• 

Parameter 

Param. 

%  of 

---  Run 

off  — 

-—  —  -  - 

ET 

Avg  Soil 

Water 

(Base) 

Value 

Base 

Value 

Sens 

Value 

Sens 

Value 

Sens 

(ins) 

(ins) 

(ins) 

BASE  SIMULATION 

8.54 

38.30 

2.371 

70 

-12.5 

8.48 

0.063 

38.37 

-0.014 

2.232 

0.469 

CN 

75 

-6.25 

8.47 

0.129 

38.37 

-0.029 

2.274 

0.655 

(80) 

85 

6.25 

8.73 

0.345 

38.12 

-0.076 

2.582 

1.424 

90 

12.5 

9.37 

0.772 

37.47 

-0.173 

3.178 

2.723 

3.325 

-30 

10.15 

-0.627 

36.74 

0.136 

1.721 

0.914 

ULE  (ins) 

4.038 

-15 

9.25 

-0.556 

37.59 

0.124 

2.044 

0.919 

(4.75) 

5.460 

15 

7.91 

-0.497 

38.94 

0.111 

2.692 

0.903 

6.175 

30 

7.52 

-0.400 

39.33 

0.089 

3.011 

0.900 

152 

-50 

16.46 

-1.855 

30.30 

0.418 

1.123 

1.053 

LAI-Days 

243 

-20 

11.16 

-1.532 

35.65 

0.346 

1.900 

0.993 

(304) 

364 

20 

6.79 

-1.025 

40.09 

0.233 

2.679 

0.650 

455 

50 

5.10 

-0.806 

41.83 

0.184 

3.092 

0.608 

RAD  I 

-20 

12.80 

-2.492 

34.01 

0.561 

1.726 

1.360 

20 

5.95 

-1.518 

40.93 

0.343 

2.816 

0.938 

TEMP 

-20 

11.18 

-1.543 

35.64 

0.348 

1.940 

0.909 

20 

7.12 

0.834 

39.75 

1.038 

2.600 

1.097 

RC  (in/hr 

0.50 

-50 

8.54 

0.000 

38.30 

0.000 

2.371 

0.000 

(1.0) 

2.00 

100 

8.54 

0.000 

38.30 

0.000 

2.371 

0.000 

FUL 

0.15 

-40 

9.15 

-0.177 

37.70 

0.039 

2.246 

0.132 

(.25) 

0.35 

40 

8.12 

-0.124 

38.73 

0.028 

2.412 

0.043 

POROS 

0.34 

-20 

8.54 

0.000 

38.30 

0.000 

2.371 

0.000 

(.43) 

0.52 

20 

8.54 

0.000 

38.30 

0.000 

2.371 

0.000 

0.04 

-50 

9.51 

-0.225 

37.34 

0.050 

2.156 

0.181 

PORSUB 

0.06 

-25 

9.06 

-0.242 

37.79 

0.054 

2.252 

0.201 

(.08) 

0.10 

25 

7.99 

-0.261 

38.86 

0.058 

2.468 

0.164 

0.12 

50 

7.25 

-0.302 

39.59 

0.067 

2.602 

0.195 
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Leaf  area  index  (LAI),  and  temperature  (TEMP)  and  radiation  (RADI) 
data  are  all  used  in  predicting  ET.  As  can  be  seen,  the  model's 
response  is  very  sensitive  to  each  of  these  parameters. 

Temperature  and  radiation  are  input  as  monthly  values.  This  sen- 
sitivity is  a  detriment  to  comparison  of  model  results  with  historical 
data  as  in  the  testing  and  validation  of  a  model.  CREAMS-WT  uses  aver- 
age monthly  temperature  and  radiation  thus  daily  or  weekly  fluctuations 
in  the  historical  period  being  simulated  cannot  be  represented.  This 
will  produce  inaccuracies  in  the  simulated  ET,  and  result  in  discrepan- 
cies between  observed  data  and  the  simulated  values.  However,  for 
general  predictive  use,  this  sensitivity  presents  no  real  drawback  to 
using  the  model.  Average  temperature  and  radiation  data  are  known  with 
reasonable  accuracy  and  will  not  change  from  one  simulation  to  the 
next.  Temperature  and  radiation  are  physical  values  that  are  generally 
well  defined,  thus  leaf  area  index  becomes  one  of  the  most  important 
parameters  to  be  estimated. 

The  base  simulation  used  in  Table  3-4  had  no  deep  seepage  (DSP=0) . 
To  evaluate  the  sensitivity  of  this  parameter,  a  second  set  of  simula- 
tions were  run  with  a  new  base  parameter  set  having  DSP  =  0.004.  The 
results  of  this  sensitivity  analysis  are  given  in  Table  3-5.  Several 
of  the  parameters  are  evaluated  again  along  with  the  new  parameter, 
DSP.  There  are  only  small  changes  in  the  sensitivity  of  the  basic 
parameters,  with  all  the  previous  observations  still  valid.  The 
parameter  DSP  is  significant  in  the  prediction  of  deep  percolation. 
Runoff  appears  to  be  only  slightly  sensitive  to  DSP,  but  comparison  of 
the  actual  values  shows  there  is  some  effect. 
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Parameter  Estimation 
The  CREAMS  manual  provides  extensive  help  for  coming  up  with  the 
necessary  parameter  values  for  the  model,  and  is  the  primary  source  for 
information  for  the  CREAMS-WT  model.  However,  there  are  several  new 
parameters,  and  the  emphasis  on  others  has  changed.  Some  comments  on 
estimating  these  parameters  for  CREAMS-WT  follows. 
CN  -  curve  number.  The  effect  of  this  comes  in  the  Smx/ULE  ratio.  A 
table  of  curve  numbers  and  corresponding  storage  are  included  in 
Appendix  A.  The  standard  curve  number  may  be  used,  or  alternately, 
when  ULE  has  been  determined,  a  curve  number  can  be  chosen  to  give 
a  desired  Smx/ULE  ratio  with  the  corresponding  runoff  curve  as  in 
Figure  3-2. 
UL( I)  -  profile  storage  by  layer  (ins).  Summing  the  seven  layers  gives 
ULE,  the  storage  in  the  root  zone.  This  should  be  a  realistic 
measure  of  the  actual  storage  in  the  profile  to  the  depth  of  RD. 
Storages  calculated  using  0.1  bar  water  content  appear  to  be 
consistently  high.  Better  estimates  may  be  found  in  the  curves  of 
depth  to  water  table  versus  storage  from  the  ARS  (Figure  2-6)  and 
the  SFWMD  (Figure  2-7).  Capece  (1984)  found  these  to  be  the  best 
storage  estimates  for  predicting  storm  runoff.  These  values  may  be 
multiplied  by  the  factor  1/ ( 1-FUL)  (FUL  is  defined  below)  as  depth 
to  water  table  may  represent  only  the  drainable  water  and  not  the 
water  held  at  "field  capacity".  An  additional  estimate  is  from 
Allen  et  al.  (1982a),  who  found  that  2  to  2.5  inches  of  stored 
water  corresponded  to  a  one-foot  change  in  the  water  table. 

This  is  an  important  parameter,  and  the  best  estimate  of  the 
active  storage  in  the  profile  should  be  used. 


RD  -  root  depth.  Generally  in  the  range  of  18  to  36  inches.  This 
should  be  to  the  top  of  the  spodic  horizon  if  one  exists. 

LAI  -  leaf  area  index.  Since  the  land  use  is  predominantly  pasture, 
this  will  not  have  to  be  changed  much.  Caution  must  be  used  in 
making  changes  because  the  model  is  sensitive  to  this  parameter. 
Calibration  may  be  needed  to  achieve  a  reasonable  annual  water  bal- 
ance. When  the  total  leaf -area-days  (printed  in  the  model  output) 
that  gives  an  appropriate  water  balance  has  been  determined,  the 
distribution  of  that  total  within  the  year  can  be  adjusted 
depending  on  the  crop. 

RC  -  saturated  hydraulic  conductivity  (in/hr).  Use  a  representative 
value.  However,  values  above  1.0  in/hr  generally  have  no  effect  as 
this  corresponds  to  moving  24  inches  of  water  per  day  through  the 
profile  which  exceeds  any  rainfall  capacity. 

FUL  -  fraction  of  ULE  held  at  field  capacity,  that  is,  the  ratio  of 
field  capacity  to  field  capacity  plus  gravitational  water. 

P0R0S  -  Has  no  effect.  Not  needed. 

PORSUB  -  porosity  of  the  lower  zone.  This  is  an  estimate  of  drainable 
or  fillable  porosity.  Can  be  estimated  from  soil  characterization 
data.  Data  from  the  Bass  West  site  indicate  a  value  in  the  range 
of  0.06  to  0.1.  This  parameter  primarily  affects  the  recovery  of 
the  water  table  after  a  long  dry  period.  If  the  value  is  too  high, 
recovery  would  be  slower  and  less  runoff  would  be  predicted.  A  low 
value  would  result  in  greater  runoff,  thus  generally  would  be  the 
conservative  choice.  The  parameter  PORS  in  equation  3-4  is  equal 
to  PORSUB  *  12. 

DSP  -  deep  seepage  rate  (in/day).  User  specified. 


CHAPTER  IV 
THE  CREAMS-WT  FIELD  SCALE  MODEL:  NUTRIENT  COMPONENT 


Development 
Description  of  the  CREAMS  Nutrient  Model 

The  CREAMS  nutrient  model  simulates  soluble  and  sediment-attached 
nitrogen  and  phosphorus  in  runoff,  considering  the  top  one  centimeter 
of  the  soil  as  the  active  zone.  Nitrogen  is  modeled  in  more  detail 
than  is  P.  A  nitrogen  budget  considers  plant  uptake,  denitrif ication, 
mineralization  of  organic  nitrogen  and  leaching  of  nitrate.  Fertilizer 
nitrogen  and  phosphorus  can  be  added  to  the  surface  or  incorporated  in 
the  profile  in  single  or  multiple  applications  during  a  year.  Since 
the  erosion  component  is  being  bypassed  for  simulation  of  flatwoods 
sites,  evaluation  and  added  enhancements  were  limited  to  the  portion  of 
the  model  describing  movement  of  soluble  nutrients. 

Soluble  nitrogen  and  phosphorus.  The  soluble  nutrients  in  the  top 
centimeter  of  the  soil  are  assumed  to  be  available  for  leaching  and 
movement  into  runoff.  Phosphorus  is  assumed  to  be  adsorbed  by  the 
soil,  thus  no  leaching  of  phosphorus  is  considered.  Similar  relation- 
ships are  used  for  modeling  both  nitrogen  and  phosphorus.  The  basic 
assumption  is  that  the  rate  of  change  in  concentration  of  soluble 
nutrient  is  proportional  to  the  difference  between  the  existing  concen- 
tration and  a  base  concentration.  This  can  be  written, 

dC/dt  =  Kl  f(t)  (Cb-C)  [4-1] 
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where  C  is  the  concentration  in  the  soil,  Cb  is  the  base  concentration, 
Kl  is  a  rate  constant  for  downward  movement,  and  f(t)  describes  infil- 
tration as  a  function  of  time.  For  nitrogen,  the  base  concentration  is 
taken  as  the  concentration  of  nitrogen  in  the  rainfall  and  for  phospho- 
rus, the  base  concentration  is  assumed  to  be  the  equilibrium  phosphorus 
concentration  of  the  soil.  Both  values  for  base  concentration  are 
input  parameters. 

For  simplification,  it  is  assumed  that  infiltration  occurs  at  a 
uniform  rate,  after  which  runoff  occurs.  Nutrient  interaction  with 
runoff  is  assumed  to  be  similar  to  that  for  leaching  (as  in  equation 
4-1),  with  a  rate  constant  K2  used  for  movement  into  runoff.  These 
relationships  are  integrated  to  give  the  concentration  at  the  end  of  a 
storm: 

C.+1  =  Cb  +  (Ct  -  Cb)  exp(-Kl*F  -K2*Q)  [4-2] 

where  subscripts  t  and  t+1  indicate  the  beginning  of  the  current  and 
the  next  time  period  respectively,  Q  is  the  total  runoff  (mm),  and 
concentrations  are  in  units  of  mg/L.  F  is  the  effective  infiltration 
(mm),  calculated  as  the  difference  between  rainfall  and  runoff  and 
additionally  subtracting  the  pore  volume  of  the  surface  layer.  Since 
the  model  uses  a  daily  time  step,  this  equation  is  used  to  calculate 
the  concentration  for  the  next  day,  given  the  predicted  depths  of 
infiltration  and  runoff  from  the  hydrology  model.  Similar  equations 
are  derived  for  the  mean  concentration  during  infiltration  and  mean 
concentration  during  runoff. 
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The  soluble  N  and  P  in  runoff  (kg/ha)  are  then  determined  by 

■RON  =  Cr  *  EXKN  *  Q  *  0.01  [4-3] 

ROP  =  Cr  *  EXKP  *  Q  *  0.01  [4-4] 

where  Cr  is  the  mean  concentration  of  the  nutrient  during  runoff.  The 
runoff  "extraction  coefficients"  are  defined  as 

EXKN  =  K2N  *  POR  *  d  [4-5] 

EXKP  =  K2p  *  POR  *  d  [4-6] 

with  similar  expressions  for  the  vertical  extraction  coefficients  EXKN1 
and  EXKP1.  POR  is  the  soil  porosity,  and  d  =  10  mm  is  the  depth  of  the 
surface  layer. 

The  phosphorus  buffering  capacity  of  the  soil  is  assumed  to  keep 
the  soluble  P  concentration  from  dropping  below  a  level  characteristic 
of  the  soil  and  adsorb  the  soluble  P  that  is  "leached"  from  the  surface 
layer.  Since  this  is  the  nutrient  of  greatest  concern  in  the  TCNS 
basin  and  since  the  assumptions  used  in  developing  the  CREAMS  P  compo- 
nent are  less  appropriate  for  flatwoods  watersheds,  the  discussion  here 
will  focus  on  the  P  component  of  the  model. 

Parameter  estimation.  In  predicting  the  movement  of  N  and  P  into 
runoff,  there  are  two  primary  input  parameters  to  be  specified,  the 
runoff  extraction  coefficient  and  the  base  concentration.  For  nitro- 
gen, the  additional  parameters  that  describe  the  other  transformations 
also  come  into  play,  but  are  less  significant.  However,  the  phosphorus 
model  is  almost  completely  described  by  the  extraction  coefficients, 
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EXKP1  and  EXKP,  the  base  buffered  soil  P  concentration,  SOLP,  and  any 
fertilizer  inputs. 

In  the  CREAMS  manual,  Frere  et  al .  (1980)  comment  that  data 
suggest  extraction  coefficients  for  pesticides  and  nutrients  range  from 
0.01  to  0.4.  Since  data  are  not  readily  available  for  defining  these 
parameters,  the  vertical  extraction  coefficients,  EXKP1  and  EXKN1,  have 
been  fixed  in  the  model  with  values  of  0.25.  This  was  done  to  simplify 
the  input  requirements,  letting  the  variability  of  the  extraction 
process  be  represented  by  the  one  remaining  coefficient.  The  value  of 
0.25  was  a  somewhat  arbitrarily  chosen  value  which  gave  reasonable 
results  with  available  data.  Frere  et  al .  (1980)  assume  that  leaching 
efficiency  is  greater  than  extraction  into  runoff,  so  the  constraint  on 
EXKP  and  EXKN  is  that  they  be  less  than  0.25. 

The  CREAMS  manual  discusses  methods  of  estimating  the  parameters 
SOLP  and  EXKP.  The  best  method  is  to  determine  the  equilibrium  P 
concentration  in  soil  samples.  The  next  best  method  indicated  is  to 
use  measured  data  from  several  storms  to  fit  the  equation, 

PPMP  =  EXKP  *  CP  [4-7] 

where  PPMP  is  the  concentration  in  runoff  (ppm),  and  CP  is  the  concen- 
tration in  the  pore  water  of  the  surface  soil  (ppm).  Unfortunately, 
these  data  are  not  readily  available,  and  very  few  model  users  will 
make  the  effort  to  acquire  them.  The  above  equation  is  simply  a  form 
of  the  basic  relationship  in  equation  4-4.  If  some  estimate  of  base  or 
background  concentrations  from  a  watershed  are  available,  the  following 
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expression  may  be  used  to  help  estimate  the  correct  range  of  values  for 
SOLP  and  EXKP, 

RP  =  SOLP  *  EXKP  *  25  [4-8] 

where  RP  (ppm)  is  a  base  concentration  in  runoff,  SOLP  has  units  of 
kg/ha,  and  the  25  is  for  units  conversion. 

Applicability  of  P  Model  to  Flatwoods  Watersheds 

The  performance  of  the  phosphorus  model  is  clearly  a  function  of 
making  good  estimates  of  EXKP  and  SOLP.  The  nutrient  load  to  the 
watershed  from  fertilizer  and  animal  waste  is  extremely  important  but 
it  is  assumed  that  these  inputs  can  be  adequately  defined.  The  soluble 
N  and  P  model  in  CREAMS  as  described  above,  especially  with  the  neces- 
sity of  estimating  the  runoff  extraction  coefficient,  is  acknowledged 
by  its  developers  as  being  the  weakest  part  of  the  nutrient  model 
(Frere  et  al.,  1980).  This  is  clearly  a  concern  for  use  in  flatwoods 
watersheds  where  the  P  nutrient  model  is  defined  solely  in  terms  of  the 
soluble  component  since  there  is  no  significant  erosion  and  sediment 
transported  P.  Several  questions  must  be  raised  in  regards  to  the 
appropriateness  of  the  P  model  and  the  reasonableness  of  the  assump- 
tions for  flatwoods  watersheds.  Is  the  basic  model  appropriate?   Is 
the  assumption  of  a  highly  buffered  soil  reasonable?  Is  the  assumption 
of  EXKP1  =  0.25  reasonable? 

The  concept  of  a  thin  surface  layer  being  the  active  zone  for 
movement  of  soluble  P  into  runoff  is  a  reasonable  one,  and  has  been 
demonstrated  (Ahuja  et  al.,  1981;  Sharpley  et  al.,  1981;  Sharpley, 


1985).  There  has  also  been  evidence  of  a  relationship  between  this 
active  depth  and  the  extraction  coefficient,  demonstrating  that  there 
may  be  some  physical  basis  to  the  simple  algorithm  used  in  CREAMS  and 
that  it  can  provide  accurate  predictions  of  mean  annual  P  concentra- 
tions (Sharpley  et  al.,  1982).  However,  this  work  has  been  done  with 
silt  loam  to  clay  soils  and  the  general  applicability  and  procedures 
for  determining  EXKP  have  not  been  established.  Thus,  while  the 
soluble  P  model  has  some  physical  basis,  its  use  generally  involves 
empirical  fitting  of  some  form  of  the  basic  relationship  shown  in 
equation  4-4. 

The  assumption  of  a  highly  buffered  soil  which  adsorbs  soluble  P 
and  maintains  a  constant  equilibrium  phosphorus  content  (the  value  of 
SOLP)  is  questionable  for  the  sandy  flatwoods  soils.  As  discussed  in 
the  literature  review  section,  data  from  many  Florida  sands  and  from 
flatwood  spodosols  in  particular,  indicate  that  P  readily  moves  down 
through  the  profile.  The  P  adsorption  isotherms  for  these  soils  have 
very  low  partition  coefficients.  Thus,  the  addition  or  removal  of  P 
from  the  soil  readily  changes  the  EPC  (equilibrium  phosphorus  concen- 
tration), whereas  for  highly  buffered  soils,  large  changes  in  P  will 
have  only  a  minor  effect  on  the  EPC.  The  assumption  of  a  highly  buf- 
fered soil  is  clearly  a  weakness  of  the  model  for  flatwoods  watersheds, 

The  extraction  coefficients  are  parameters  that  describe  the 
efficiency  of  nutrient  extraction  into  percolating  water  and  surface 
runoff.  The  phosphorus  model  further  assumes  that  there  is  no  actual 
movement  of  P  down  through  the  soil  profile,  as  all  phosphorus  extrac- 
ted by  vertical  flow  is  rapidly  bound  up  by  the  soil.  Phosphorus 
removed  from  the  soluble  pool  by  this  process  is  not  considered  any 
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further  by  the  model.  Thus,  while  the  derivation  of  the  above  algo- 
rithms considers  the  extraction  coefficients  to  be  descriptive  of  the 
flow  process,  in  the  case  of  P,  what  is  being  modeled  is  actually  the 
processes  of  adsorption  and  precipitation. 

In  flatwoods  soils,  the  removal  of  soluble  P  from  the  active  run- 
off extraction  layer  in  the  profile  is  a  combination  of  physical  move- 
ment and  chemical  adsorption.  Due  to  the  low  buffering  capacity  of  the 
soils,  it  would  be  expected  that  the  set  value  of  EXKP1  =  0.25  may  be 
too  high.  Initial  simulations  showed  that  soil  P  was  rapidly  reduced 
to  the  base  P  level  (specified  by  SOLP)  after  a  fertilizer  application, 
while  data  indicate  that  fertilizer  applications  may  contribute  to 
runoff  loads  for  several  months  (Goldstein,  1986). 

CREAMS-WT  Enhancements  for  Flatwoods  Watersheds 

Denitrification.  The  CREAMS  denitrification  algorithm  is  a 
function  of  the  number  of  days  of  drainage  from  the  bottom  layer  of  the 
profile.  This  follows  the  assumption  that  while  drainage  is  occurring 
and  the  soil  profile  is  above  field  capacity,  anaerobic  conditions  will 
be  present  in  the  soil  thus  enabling  denitrification.  Under  the  modi- 
fied hydrologic  model  of  CREAMS-WT,  days  when  saturated  conditions 
exist  in  the  upper  profile  are  determined  explicitly  as  a  result  of 
following  the  level  of  the  water  table.  The  number  of  days  the  water 
table  is  in  the  upper  profile  is  now  passed  for  use  by  the  denitrifica- 
tion algorithm  in  place  of  the  original  CREAMS  value.  This  modifica- 
tion was  actually  made  in  the  hydrology  model  but  its  effect  is  seen  in 
the  nutrient  model. 


Phosphorus  in  rainfall.  The  ability  to  specify  the  concentration 
of  phosphorus  in  rainfall  has  been  added  to  the  model.  This  input  is 
handled  the  same  as  the  nitrogen  loading  in  rainfall  and  is  the  same  in 
concept  to  a  fertilizer  application.  The  load  in  kg/ha  is  calculated 
knowing  the  amount  of  rainfall.  This  input  load  is  then  added  to  the 
nutrient  content  in  the  active  layer  before  calculating  extraction  by 
leaching  and  runoff. 

Soluble  phosphorus  model.  Phosphorus  removal  from  the  active 
layer  is  controlled  by  two  factors,  the  value  of  the  extraction  coeffi- 
cient, EXKP1 ,  and  the  volume  of  effective  infiltration.  While  directly 
changing  the  value  of  EXKP1  is  possible,  an  indirect  way  of  reducing  P 
removal  was  implemented.  The  depth  of  the  active  layer  was  set  in  the 
CREAMS  model  at  1  cm.  Increasing  this  depth  affects  the  model's  re- 
sponse in  two  ways.  First,  it  results  in  a  larger  active  P  pool,  and 
extraction  by  "vertical"  removal  and  by  runoff  has  a  proportionally 
smaller  effect  on  the  P  pool.  The  result  is  that  the  effective  extrac- 
tion coefficients  are  inversely  proportional  to  the  increased  thickness 
of  the  active  layer.  Second,  an  increase  in  the  depth  of  the  active 
layer  reduces  the  amount  of  "effective  infiltration"  (F  in  equation 
4-2)  that  is  used  to  calculate  P  removal  from  the  active  pool. 

There  are  several  advantages  in  using  this  approach  to  modify  the 
response  of  the  P  model.  1)  It  has  the  equivalent  effect  of  reducing 
EXKP1.  2)  Extraction  coefficient  values  remain  in  the  same  range  as 
those  discussed  in  the  CREAMS  manual.  3)  A  reduction  in  the  volume  of 
effective  infiltration  is  reasonable  because  of  the  effects  of  a  high 
water  table.  An  increase  in  depth  of  the  active  layer  is  also 
physically  reasonable  for  the  flatwoods  system,  where  the  majority  of 


runoff  passes  through  the  upper  soil  profile  on  its  way  to  a  ditch  or 
stream. 

Testing  and  Calibration 

Issues  in  calibration  and  verification.  Measured  data  all  have 
some  degree  of  error  or  uncertainty  associated  with  them.  There  are 
several  sources  of  uncertainty  in  water  quality  data  from  the  UD/RDP 
sites  that  need  to  be  recognized,  especially  in  regard  to  model  cali- 
bration and  verification. 

1.  Dates  and  amounts  of  fertilizer  applications  are  not  known 
exactly.  The  timing  of  fertilizer  application  and  a  rainfall 
event  can  significantly  affect  the  amount  of  nutrients  in 
runoff. 

2.  The  activity  of  cattle  cannot  be  predicted  by  the  model,  but 
can  have  a  dramatic  impact  on  the  nutrient  concentrations  in 
waterways  (Goldstein,  1986;  Ritter  and  Allen,  1982).  Thus, 
measured  in-stream  concentrations  may  at  times  not  strictly 
represent  the  quality  of  runoff  from  the  watershed,  but  may  be 
due  in  part  to  the  activity  of  cattle  in  streams  and  ditches. 

3.  Water  quality  samples  were  collected  biweekly.  The  adequacy 
of  representing  the  true  time  series  concentrations  with  this 
sampling  frequency  must  be  questioned.  Indeed,  on  the  SEZ 
dairy  watershed,  Goldstein  (1986)  noted  that:  "The  sensiti- 
vity of  P  concentrations  to  events  on  the  watershed  began  to 
show  up  vividly  when  the  resolution  of  the  sampling  interval 
was  increased  from  weekly  to  daily."  The  low  sampling  fre- 
quency introduces  uncertainty  in  the  calculated  nutrient  loads 
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and  flow  weighted  concentrations  which  are  the  primary  data 
used  to  compare  with  model  predictions. 
4.  Errors  in  flow  measurement  bias  the  calculated  flow  weighted 
concentrations.  Time  series  concentrations  generally  increase 
with  higher  flow  in  these  flatwoods  watersheds  (Goldstein, 
1986).  If  flow  is  over-estimated,  then  the  calculated  flow 
weighted  concentrations  will  be  inadvertently  high.  Since 
runoff  volumes  for  the  UD/RDP  sites  appear  to  be  high  for 
1982,  flow  weighted  concentrations  may  also  be  elevated. 
Measured  water  quality  data  consist  of  samples,  taken  weekly  or 
biweekly,  from  which  the  concentration  of  various  constituents  is 
determined.  These  time  series  concentrations,  in  combination  with 
measured  flow  data,  are  used  to  estimate  the  nutrient  load  leaving  the 
watershed.  From  the  collected  data,  three  forms  of  the  data  can  be 
used:  time  series  concentrations,  flow  weighted  concentrations,  and 
nutrient  load. 

The  CREAMS-WT  nutrient  model  predicts  flow  weighted  concentrations 
and  nutrient  loads.  Thus,  these  two  types  of  data  are  available  for 
use  in  calibration  and  verification  of  the  model.  Of  course,  nutrient 
load  is  a  direct  product  of  the  flow  weighted  concentration  and  the 
volume  of  flow  over  a  given  period.  Thus,  flow  weighted  concentrations 
and  loads  predicted  by  a  model  can  both  accurately  represent  observed 
values  only  when  the  predicted  runoff  matches  observed  runoff.  While 
the  ultimate  objective  in  modeling  is  to  have  all  aspects  of  the  model 
match  perfectly,  this  is  realistically  not  possible.  The  emphasis  in 
model  calibration  is  to  focus  on  the  predicted  variable  that  is  of 
greatest  interest.  Thus,  whether  to  focus  on  concentrations  or  loads 
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will  depend  on  the  objectives  of  the  particular  modeling  exercise.  The 
following  reasons  suggest  an  emphasis  on  concentrations. 

1.  Measurement  of  nutrient  concentrations  is  much  more  accurate 
than  runoff  measurement.  This  is  especially  true  for  flat- 
woods  watersheds  where  runoff  measurement  is  plagued  with 
problems  in  flow  gaging  and  identifying  watershed  boundaries. 
Concentrations  will  thus  be  a  more  reliable  basis  for 
comparison  with  simulated  output. 

2.  In  calibration  or  verification,  errors  in  simulated  and 
observed  runoff  may  compound  differences  in  computed  loading. 
When  the  model  is  used  for  predictive  purposes,  and  the  simu- 
lated concentrations  and  runoff  are  both  representative,  then 
computed  loadings  will  also  be  representative.  The  represent- 
ativeness of  the  hydrology  model  has  been  demonstrated,  in 
both  short  term  comparisons  and  long  term  water  balance  pre- 
dictions. Thus  emphasis  on  concentrations  in  the  nutrient 
model  will  assure  that  predicted  loadings  are  representative. 
While  it  would  be  most  desirable  to  use  time  series  concentra- 
tions, these  are  not  available  from  CREAMS-WT,  so  comparison 
is  made  with  the  next  best  choice,  flow-weighted 
contentrations. 

The  approach  to  calibration  and  verification  is  in  part  directed 
by  the  objectives  for  the  model.  CREAMS-WT  will  be  used  for  a  priori 
simulation  of  ungaged  watersheds,  and  is  intended  to  be  used  in  evalua- 
ting the  relative  differences  between  management  alternatives  applied 
to  the  watershed.  Thus  calibration  was  limited  to  determining  the 
proper  magnitude  of  the  critical  input  parameters.  Demonstrating  that 


92 

the  model  predictions  fall  in  the  range  of  observed  data  is  more  impor- 
tant than  fine  tuning  the  parameters  to  obtain  a  "best  fit". 

Simulation  of  Bass  West  pasture.  Data  from  the  Bass  West 
watershed  were  used  in  the  initial  evaluation  and  calibration  of  the 
modified  model.  The  objective  of  this  process  was  to  determine  the 
depth  of  the  active  layer  and  the  relative  magnitude  of  the  runoff 
extraction  coefficients.  Because  of  the  uniformity  of  the  topography 
and  soils  of  the  region,  it  is  expected  that  these  values  will  be 
characteristic  of  the  entire  basin. 

Parameter  values  for  simulating  the  Bass  West  site  were  estimated 
from  the  CREAMS  manual  and  from  knowledge  of  the  area.  The  parameter 
file  is  listed  in  Appendix  C.  The  predicted  flow  weighted  concentra- 
tions are  shown  in  Table  4-1  along  with  the  observed  values.  The  most 
important  parameters  for  the  nutrient  simulation  are  the  values  of 
SOLP,  the  runoff  extraction  coefficients,  EXKP  and  EXKN,  and  the  ferti- 
lizer and  animal  waste  applied.  Exact  dates  and  amounts  of  fertilizer 
applied  to  this  pasture  are  unknown,  but  typical  practice  is  to  apply 
300  lb/acre  of  20-10-10  or  similar  N-P-K  formulation  fertilizer  in  the 
fall  (Goldstein,  1986).  This  is  equivalent  to  67  kg  N/ha  and  33.6  kg 
P/ha  per  year.  This  pasture  has  been  fertilized  for  many  years,  as  a 
vegetable  field  and  as  an  improved  pasture.  As  a  result,  residual  P  in 
the  soil  is  higher  than  would  be  found  in  native,  or  semi-improved, 
unfertilized  pastures.  The  chosen  value  for  SOLP  =  0.15  is  within  the 
range  of  suggested  values  in  the  CREAMS  manual. 

For  the  Bass  West  watershed,  a  value  of  0.12  for  the  runoff  ex- 
traction coefficients  gave  reasonable  results.  Background  concentra- 
tions of  P  in  runoff  from  Bass  West  are  around  0.5  mg/L.  Thus  the 
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Table  4-1.  Bass  West  pasture  simulated  and  observed  average  annual 
flow  weighted  concentrations. 


FLOW  WEIGHTED  CONCENTRATION 

TP TN 

Simul   Observ    Simul   Observ 


(mg/L) 


1980  1.03  1.25  1.88  2.07 

1981  1.24  0.68  1.33  3.07 

1982  0.66  0.85  1.58  2.20 

Average  0.96  0.93  1.80  2.45 


parameter  values,  when  substituted  in  equation  4-8,  indicate  a  back- 
ground concentration  similar  to  that  observed.  Values  for  extraction 
coefficients  that  appear  in  the  sample  simulations  in  the  CREAMS  manual 
range  from  0.075  -  0.10.  The  AGNPS  model  (Young  et  al.,  1985)  which 
uses  the  CREAMS  nutrient  algorithms,  has  the  extraction  coefficients 
fixed  with  values  of  0.075.  Higher  extraction  efficiency  might  be 
expected  in  these  sandy  soils,  but  direct  comparison  is  difficult  to 
make  because  of  the  changes  in  the  model. 

While  predicting  average  annual  flow  weighted  concentration  is  the 
objective  of  the  nutrient  model,  comparison  of  monthly  values  will  help 
indicate  how  well  the  model  is  representing  the  processes  that  control 
the  response  of  the  watershed.  In  the  calibration  process,  average 
monthly  flow  weighted  concentrations  were  helpful  in  determining  what 
value  to  use  for  the  depth  of  the  active  layer  (DEPTH).  This  depth 
controls  the  rate  of  phosphorus  removal  by  fixation  and  percolation. 
The  unmodified  CREAMS  removed  phosphorus  so  rapidly  that  predicted  soil 
phosphorus  was  rapidly  returned  to  the  equilibrium  concentration. 
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Figure  4-1  compares  CREAMS-WT  and  CREAMS  predictions  (using  the 
same  parameter  set)  with  data  from  Goldstein  (1986).  The  fall 
application  of  fertilizer  is  during  the  dry  season.  The  data  clearly 
indicate  that  some  of  this  phosphorus  remains  in  the  soil  in  available 
form,  and  is  washed  out  in  runoff  with  the  onset  of  the  wet  season. 
The  difference  between  the  predictions  of  the  two  models  is  primarily 
due  to  the  reduction  of  the  effective  "vertical"  extraction  coefficient 
in  CREAMS-WT,  which  is  a  function  of  the  depth  of  the  active  layer 
(DEPTH).  Values  of  DEPTH  from  2  to  4  cm  were  simulated,  with  the  value 
of  3  selected  for  inclusion  in  the  model.  The  changes  made  in  the 
nutrient  model  clearly  enable  it  to  better  represent  the  degree  of 
soil/phosphorus  interaction  occurring  in  these  poorly  buffered  soils. 

Verification 

Two  additional  sites  from  the  UD/RDP  were  simulated.  The  Bass 
East  pasture  is  used  for  grazing  beef  cattle  and  is  not  fertilized, 
thus  has  low  nutrient  loading.  The  SEZ  dairy  watershed  has  the 
heaviest  nutrient  loading  of  all  the  UD/RDP  sites,  with  waste  from  the 
dairy  operation  and  annual  fertilization  of  pasture  and  hayland.  Thus, 
these  two  sites  bracket  the  nutrient  loading  of  the  Bass  West  pasture. 
This  provides  some  indication  of  the  applicability  of  the  model  over 
the  range  of  conditions  which  are  encountered  in  the  TCNS  basin. 

The  following  simulations  were  made  without  calibration.  Values 
for  SOLP,  EXKP,  and  EXKN  were  estimated  by  comparison  with  those  used 
for  the  Bass  West  simulation. 
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SEZ  Dairy 

The  SEZ  dairy  watershed  is  divided  between  two  main  land  uses, 
holding  and  staging  pastures  for  the  dairy  herd  (56  ha),  and  other 
pasture  used  for  hay  and  grazing  dry  cows  and  calves  (217  ha).  To 
simulate  this  entire  area  as  a  unit  requires  assuming  some  average 
condition  that  represents  the  whole.  The  extraction  coefficients  could 
reasonably  be  expected  to  be  fairly  uniform  over  both  land  uses  since 
there  are  no  clear  differences  in  topography  and  drainage  between  the 
two.  Ditching  density  is  lower  than  Bass  West  pasture,  so  a  value  of 
0.10  was  selected  for  the  extraction  coefficients. 

Nutrient  loading  as  fertilizer  and  animal  waste,  as  well  as  resi- 
dual phosphorus  in  the  soil  varies  between  the  two  land  uses.  In 
general,  the  holding/staging  areas  receive  a  significant  waste  load, 
and  the  hay/pasture  land  is  fertilized.  Typical  fertilization  on  the 
hay/pasture  land  is  the  application  of  300  lb/acre  of  16-8-4  or  similar 
N-P-K  formulation  fertilizer  in  April  or  May.  This  is  roughly  equiva- 
lent to  27  kg  P  and  54  kg  N  per  ha  per  year.  In  the  simulation,  ferti- 
lizer is  applied  on  May  10  of  each  year,  a  date  selected  arbitrarily 
since  exact  dates  of  application  are  unknown. 

Nutrients  from  the  manure  of  the  dairy  herd  must  also  be  consi- 
dered. Dairy  cows  receive  much  of  their  feed  from  sources  external  to 
the  pastures  in  which  they  are  kept.  Thus  the  manure  is  a  nutrient 
application  with  a  source  external  to  the  pasture.  This  is  in  contrast 
with  a  beef  pasture  where  the  cattle  receive  little  extra  feed.  Here, 
the  manure  is  simply  one  stage  in  the  nutrient  cycling  of  the  pasture 
ecosystem  and  should  not  be  considered  as  a  nutrient  input  to  the 
watershed. 
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There  are  about  500  cows  in  the  SEZ  dairy  milking  herd.  Goldstein 
(1986)  estimates  phosphorus  loads  of  0.096  lb/day  per  cow  which  is 
comparable  to  other  estimates  (ASAE,  1983).  Distributing  this  load 
over  the  entire  area  gives  a  P  load  of  2.38  kg/ha/month.  Nitrogen  in 
dairy  manure  is  estimated  at  75  kg/yr  per  cow  (ASAE,  1983).  Distrib- 
uted uniformly  over  the  watershed  this  converts  to  a  N  load  of  11.4 
kg/ha/month.  In  the  model,  nutrient  load  from  animal  waste  is  applied 
on  a  monthly  basis  as  a  compromise  in  representing  what  is  really  a 
continuous  application.  (This  issue  will  be  addressed  in  the  sensi- 
tivity analysis  section.) 

The  parameter  file  for  the  SEZ  site  is  included  in  Appendix  C. 
The  measured  and  predicted  average  annual  flow  weighted  concentrations 
are  shown  in  Table  4-2.  Predictions  are  again  within  the  range  of 
observed  values.  One  additional  nutrient  source  in  the  SEZ  dairy 
watershed  which  is  not  directly  accounted  for  by  the  model  is  the  dis- 
charge from  the  dairy  lagoon.  Goldstein  (1986)  estimates  that  lagoon 
discharge  comprises  only  5%  of  the  total  animal  waste  load  on  the 
watershed  which  may  not  be  significant  in  years  with  normal  rainfall. 
However,  in  dry  years,  runoff  from  the  pasture  is  greatly  decreased 
while  the  volume  of  water  flushed  into  the  lagoon  will  remain  fairly 
constant.  Thus  lagoon  discharge  could  contribute  significantly  to  the 
nutrient  load  observed  at  the  watershed  outlet  in  dry  years. 

Bass  East  Pasture 

The  residual  phosphorus  in  the  surface  horizon  at  this  site  is 
expected  to  be  lower  than  that  of  the  neighboring  West  pasture,  which 
has  been  fertilized  for  many  years.  Correspondingly,  a  lower  value  of 
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Table  4-2.  SEZ  dairy  simulated  and  observed  average  annual  flow 
weighted  concentrations. 


FLOW  WEIGHTED  CONCENTRATION 

TP --  TN 

Simul   Observ    Simul   Observ 


(mg/L) 


1980  2.56  3.10  4.21  6.09 

1981  1.69  2.27  5.27  3.70 

1982  1.84  1.71  4.87  4.27 

Average  1.84  2.36  4.93  4.69 


Source:  Goldstein  (1986) 


Table  4-3.  Bass  East  pasture  simulated  and  observed  average  annual 
flow  weighted  concentrations. 


FLOW  WEIGHTED  CONCENTRATION 

TP TN 

Simul   Observ    Simul   Observ 




[mg/L) 

1980 
1981 
1982 

0.30 
0.29 
0.28 

0.38 
0.15 

1.32 
1.22 
1.21 

1.93 
1.78 

Average 

0.29 

0.27 

1.29 

1.86 

Source:  Goldstein  (1986) 
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SOLP  =  0.10  was  chosen.  The  extraction  coefficient  is  also  expected  to 
be  lower  since  ditching  density  is  lower.  A  value  of  EXKP  =  0.10  was 
chosen.  The  parameter  file  is  included  in  Appendix  C,  and  the 
simulation  results  are  shown  in  Table  4-3.  Measured  runoff  from  this 
site  was  negligible  in  1981  (thus  no  "observed"  value  in  Table  4-3), 
but  a  very  small  amount  (less  than  1  cm)  was  predicted  by  the  model. 

Discussion 

The  Bass  East  pasture  and  SEZ  dairy  sites  were  simulated  without 
calibration.  The  results  shown  are  from  the  initial  simulation  made 
after  selecting  parameters  as  described  earlier.  The  ability  of  the 
model  to  simulate  these  two  sites  without  calibration  demonstrates  its 
applicability  for  use  in  a  priori  simulation  of  sites  for  the  purpose 
of  evaluating  the  relative  effects  of  changes  in  management  practices. 

The  reason  the  model  responds  predictably  is  that  the  critical 
parameters  can  be  determined  with  some  certainty.  The  nutrient  loading 
to  the  watershed  from  animal  waste  or  fertilizer  is  very  important  and 
is  generally  well  defined.  The  extraction  coefficient  is  not  expected 
to  vary  much  between  flatwoods  watersheds.  Values  should  range  from 
0.08  for  low  ditching  to  0.12  for  fields  with  high  ditching  densities. 
The  available  soil  phosphorus  (and  nitrogen)  is  harder  to  determine. 
If  data  or  estimates  of  base  concentrations  in  runoff  are  available 
simple  calibration  can  be  done.  Fortunately,  this  parameter  has  its 
greatest  effect  on  the  response  of  the  model  when  there  is  no  nutrient 
load  applied.  Since  areas  with  low  nutrient  loading  are  of  least 
importance  in  pollution  abatement  programs,  the  effects  of  uncertainty 
in  this  parameter  are  minimized. 
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Sensitivity  Analysis 

The  sensitivity  analysis  of  selected  nutrient  model  parameters  was 
performed  using,  as  the  hydrology  input,  the  twenty  year  simulation 
that  was  the  base  for  the  hydrology  sensitivity  analysis.  Relative 
sensitivity  was  calculated  using  equation  3-17.  A  relative  sensitivity 
of  1.0  means  that  a  percentage  change  in  the  parameter  produces  an 
equivalent  percent  change  in  the  output.  The  CREAMS  manual  contains  a 
sensitivity  analysis  of  all  the  CREAMS  parameters.  This  can  serve  as  a 
general  reference,  although  each  parameter  set  will  give  slightly  dif- 
ferent results.  The  primary  parameters  affected  by  the  changes  made  to 
CREAMS-WT  are  evaluated  here. 

The  base  parameter  set  used  for  the  sensitivity  analysis  (included 
in  Appendix  C)  is  similar  to  the  Bass  West  pasture  parameter  set.  This 
was  used  since  it  is  representative  of  a  typical  land  use.  The  sensi- 
tivity analysis  results  are  given  in  Table  4-4. 

As  would  be  expected  from  observing  the  controlling  relationship 
(equation  4-5),  nutrient  yield  is  directly  proportional  to  EXKP  as 
indicated  by  the  constant  value  of  the  relative  sensitivity  for 
different  parameter  values.  The  relative  sensitivity  of  EXKP  is  the 
highest  of  all  the  P  parameters.  The  N  extraction  coefficient,  EXKN, 
is  much  less  sensitive,  but  is  still  significant  in  the  N  prediction. 

The  nutrients  added  as  fertilizer  are  also  significant  to  the 
model  predictions,  more  so  for  P  than  for  N.  The  soil  equilibrium 
phosphorus  content,  SOLP,  is  less  sensitive  than  EXKP  and  fertilizer 
amounts.  Note  that  the  equivalent  parameter  for  N  has  no  effect  on 
model  results.  The  importance  of  SOLP  will  depend  on  the  amount  of 
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Table  4-4.  Relative  sensitivity  (%  change  in  result  divided  by  I 

change  in  parameter)  of  selected  CREAMS-WT  nutrient  model 
parameters,  based  on  the  predicted  average  annual  N  and  P 
in  runoff. 


Parameter 

%   of 

P  in 

runoff 

N  in 

runoff 

(Base  value) 

Value 

Base 

Value 

Sens. 

Value 

Sens. 

(kg/ha) 

(kg/ha) 

Base  simulat 

ion 

2.46 

5.16 

0.06 

-50 

1.34 

0.91 

4.28 

0.34 

EXKP.EXKN 

0.09 

-25 

1.92 

0.87 

4.76 

0.32 

(0.12) 

0.15 

25 

2.96 

0.81 

5.51 

0.27 

0.18 

50 

3.43 

0.79 

5.82 

0.25 

0.06 

-50 

2.10 

0.29 

5.16 

0.00 

S0LP.S0LN 

0.09 

-25 

2.28 

0.29 

5.16 

0.00 

(0.12) 

0.15 

25 

2.64 

0.29 

5.16 

0.00 

0.18 

50 

2.82 

0.29 

5.16 

0.00 

2.00 

-33.3 

1.97 

0.59 

4.25 

0.53 

DEPTH  (cm) 

2.50 

-16.7 

2.21 

0.59 

4.70 

0.54 

(3.0) 

3.50 

16.7 

2.74 

0.70 

5.72 

0.64 

4.00 

33.3 

3.09 

0.77 

6.18 

0.59 

0.00 

-100 

2.26 

0.08 

5.16 

0.00 

RCP  (ppm) 

0.05 

-50 

2.36 

0.08 

5.16 

0.00 

(0.10) 

0.15 

50 

2.55 

0.08 

5.16 

0.00 

0.20 

100 

2.65 

0.08 

5.16 

0.00 

Fertilizer 

-100 

0.92 

0.62 

3.16 

0.39 

(kg/ha) 

-50 

1.69 

0.63 

4.16 

0.39 

(N  -  67.2) 

50 

3.22 

0.63 

6.17 

0.39 

(P  -  33.6) 

100 

3.99 

0.63 

7.17 

0.39 
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fertilizer  applied.  If  no  fertilizer  is  applied,  the  relative  sensi- 
tivity of  SOLP  would  be  close  to  a  value  of  1.0.  As  fertilization 
rates  increase,  the  model's  sensitivity  to  SOLP  will  decrease.  This 
same  response,  though  not  nearly  as  dramatic,  is  expected  with  EXKP. 
A  base  fertilizer  amount  was  used  in  this  analysis  because  high 
nutrient  loading  is  the  land  use  of  greatest  concern. 

Both  N  and  P  yield  are  sensitive  to  the  depth  of  the  active  layer, 
DEPTH.  This  parameter  is  not  one  that  will  be  changed  by  the  model 
user  but  was  included  here  to  help  evaluate  its  effect  on  the  model 
results.  The  effect  of  phosphorus  in  rainfall  (RCP)  is  also  minimal. 
However,  in  the  absence  of  added  fertilizer,  this  parameter  would  also 
be  expected  to  be  more  significant. 

In  addition  to  analyzing  the  sensitivity  of  the  above  parameters, 
additional  aspects  of  the  model's  response  were  evaluated.  In  the 
model,  animal  waste  must  be  added  to  the  soil  as  fertilizer  applica- 
tions. Thus,  while  animal  waste  is  actually  distributed  continuously, 
it  must  be  applied  in  the  model  in  discrete  events.  Several  simula- 
tions were  run  to  determine  how  many  applications  were  required  during 
a  year  for  the  model  response  to  approximate  the  actual  uniform 
distribution. 

The  model's  response  to  distributing  a  constant  annual  load  from  1 
to  40  times  through  the  year  can  be  seen  in  Table  4-5.  Multiple  appli- 
cations were  evenly  spaced  through  the  year.  Single  applications  were 
made  on  the  day  indicated.  The  most  important  factor  affecting  the 
runoff  nutrient  load  is  the  timing  of  rainfall/runoff  events  following 
fertilizer  application.  Due  to  the  variability  of  runoff  throughout 
the  year,  the  fewer  the  fertilizer  applications,  the  more  that 
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variability  will  be  reflected  in  the  predicted  nutrient  loads.  This  is 
demonstrated  by  the  variability  in  the  response  of  the  single  applica- 
tion (Table  4-5).  This  variability  in  response  is  expected  to  decrease 
with  an  increase  in  the  number  of  times  the  fertilizer  load  is  distrib- 
uted throughout  the  year. 

While  the  actual  variability  in  model  response  has  not  been  defin- 
ed here,  the  general  response  of  the  model  is  indicated.  In  choosing 
the  number  of  applications  per  year  we  have  on  one  hand  the  desirabil- 
ity of  adequately  representing  the  continuous  distribution,  and  on  the 
other,  the  desirability  of  keeping  the  model  efficient,  both  in  terms 
of  input  time  and  computation.  Based  on  this,  twelve  applications  per 
year  are  recommended.  The  total  annual  load  is  also  a  factor,  with 
numerous  applications  less  important  for  smaller  loads. 


Table  4-5.  Model  response  to  the  time  and  number  of  applications  for 
a  uniform  annual  load  of  67.2  kg  N/ha  and  33.6  kg  P/ha. 


Number  of 

Phosphorus 

Nitrogen 

Applications 

in  runoff 

in  runoff 

(kg/ha) 

(kg/ha) 

1  (Apr  20) 

1.65 

3.79 

1  (Sep  1) 

7.20 

12.15 

1  (Nov  1) 

2.46 

5.16 

3 

3.81 

6.85 

6 

3.61 

6.46 

12 

3.52 

6.56 

20 

3.61 

6.87 

40 

3.63 

6.94 

CHAPTER  V 
BASIN  DELIVERY  MODEL 


Concepts 
The  objective  of  the  basin  scale  model  is  to  provide  some  estimate 
of  the  water  quality  impact  of  a  basin-wide  BMP  program  at  the  basin 
outlet.  In  combination  with  CREAMS-WT  as  the  field  scale  model,  the 
basin  model  must  provide  the  framework  for  combining  the  individual 
field  scale  predictions,  and  provide  additional  functions  that  are 
beyond  the  capacity  of  CREAMS-WT.  The  BASIN  model  was  developed  to 
perform  the  following  functions: 

1)  Predict  nutrient  delivery  to  sub-basin  and  basin  outlets  using 
the  simulated  CREAMS-WT  edge-of-stream  loads; 

2)  Model  the  effects  of  BMPs  that  cannot  be  simulated  in  CREAMS-WT 
(fencing,  impoundments); 

3)  Compute  nutrient  loads  due  to  cattle  in  streams  and  wetlands,  and 
for  wetlands,  attenuate  those  loads  in  flow  through  the  wetland 
to  streams; 

4)  Provide  background  loads  for  forests,  native  range,  and  other 
non-agricultural  land  use  in  the  basin. 

As  with  the  CREAMS-WT  model,  the  BASIN  model  was  designed  with  the  goal 
of  using  minimal  calibration. 
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Spatial  Discretization 

The  U.S.  Army  Corps  of  Engineers  has  developed  an  extensive 
geographic  database  of  the  TCNS  and  Lower  Kissimmee  River  basins.  A 
200  x  250  meter  grid  defines  five-hectare  (12.35  acre)  cells,  giving  a 
total  of  9280  in  the  TCNS  basin.  Data  for  each  cell  include  basin  and 
sub-basin  codes,  general  soil  group,  hydrologic  soil  group  (A/B/C/D), 
presence  of  a  stream  in  the  cell,  and  land  use  code.  This  'cell'  is 
the  basic  land  unit  used  by  the  model  for  applying  BMPs  and  in  water 
quality  modeling.  This  database  provides  the  basic  information  used 
for  "field"  scale  predictions  and  also  provides  the  spatial  relation- 
ship used  by  the  BASIN  model. 

Time  Scale 

The  BASIN  model  is  designed  to  give  an  estimate  of  the  expected 
long  term  average  annual  response.  The  average  annual  nutrient  load 
from  each  agricultural  cell  is  determined  from  a  20-year  CREAMS-WT 
simulation.  Using  20  years  of  precipitation  that  includes  the  natural 
variability  helps  ensure  that  the  simulation  is  representative  of  the 
average  long  term  response  of  the  system. 

The  BASIN  model  operates  on  average  annual  values  in  estimating 
the  load  at  the  basin  outlet  and  thus  can  best  be  described  as  a 
"delivery"  model.  Assumptions  made  in  the  development  of  the  model 
components  are  designed  to  reflect  average  annual  conditions,  thus  will 
less  accurately  represent  the  actual  annual  variability  in  the  system 
response.  This  approach  was  used  in  the  BASIN  model  because  it:  1) 
utilizes  available  data,  2)  is  computationally  efficient,  3)  enables 
changes  in  individual  fields  to  be  evaluated  quickly,  4)  produces 
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results  that  can  be  used  readily  in  other  types  of  analysis  (cost- 
effectiveness/optimization)  since  interactions  between  loads  from 
individual  cells  are  not  considered. 

Components 
A  basic  flow  chart  of  the  BASIN  model  is  shown  in  Figure  5-1.  A 
description  of  the  concepts  and  the  algorithms  used  in  implementing  the 
different  model  components  follows. 

Nutrient  Attenuation 
Algorithm 

A  first-order  rate  equation  was  used  for  modeling  the  attenuation 
of  wetland  animal  waste  load  in  flow  through  the  wetland  to  the  nearest 
stream,  edge-of-stream  loads  to  sub-basin  outlets,  and  channel  flow 
from  sub-basin  outlets  to  the  basin  outlet.  The  form  of  the  equation 
used  is 


C  =  (CQ-Cb)  e"k0  +  Cb  [5-1] 


where: 

C  =  concentration  after  attenuation,  mg/1 

C  =  input  concentration  from  cell,  mg/1 
o    r 

Ch  =  background  concentration,  mg/1 

D  =  distance,  m 

k  =  rate  coefficient 
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Read  Cell /BMP  data 


forest /range/ urban 


wetland 


I  Load  is  func. 
of  bordering 
pasture  and 
cattle  density 


Attenuate  to 
edge  of  stream 


Background 
load 


agricultural 


Load  from  CREAMS 


Lot  retention/ 
Impoundment 


Add  load  due  to 
cattle  in  stream 


Attenuate  to  w/s 
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Attenuate  to 
basin  outlet 


Summarize  output 


Figure  5-1.  Flow  chart  of  the  BASIN  model 
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The  first  order  rate  equation  has  often  been  used  in  describing 
natural  processes.  This  is  true  for  nutrient  uptake  in  flow  through 
wetlands  (Kadlec  and  Hammer,  1982;  Davis,  1981),  in-stream  nutrient 
interactions  (Brown  and  Barnwell,  1985;  Huber  et  al . ,  1976),  denitri- 
fication  (Reddy  et  al.,  1980),  and  pollutant  removal  from  overland  flow 
(Smith,  1982;  Abernathy  et  al.,  1985).  The  rate  equation  is  most  often 
expressed  as  a  function  of  time,  which  in  a  hydraulic  system  represents 
residence  or  reaction  time.  However,  with  the  inclusion  of  a  velocity 
relationship,  the  expression  can  be  written  in  terms  of  distance. 

The  inclusion  of  background  values  (Cu  in  equation  5-1)  in  the 
model  is  important  as  these  values  represent  the  minimum  levels  that 
can  be  achieved  under  any  reasonable  set  of  land  use  practices.  Note 
that  a  first  order  model  with  the  base  concentration  subtracted  will 
have  a  higher  rate  coefficient  than  the  similar  equation  without  the 
base  subtracted.  Values  for  background  concentration  and  for  the 
coefficient,  k,  were  estimated  from  the  literature,  primarily  from  data 
compiled  in  water  quality  studies  done  in  South  Florida  flatwoods 
watersheds. 
Background  concentrations 

Although  baseflow  concentrations  are  not  constant  through  time, 
data  from  watersheds  in  the  Lower  Kissimmee  River  basin  indicate  pre- 
dictable averages  in  wet  season  flow  from  undeveloped  watersheds. 
Phosphorus  concentrations  in  tributaries  to  the  Lower  Kissimmee  River 
range  from  0.253  to  0.026  mg/L,  mean  wet  season  flow  (Federico,  1982). 
The  low  value  was  from  a  watershed  that  was  primarily  in  native 
rangeland.  Other  watersheds  with  similar  land  use  were  among  the  lower 
values  in  the  range  of  P  values  observed.  The  Boney  Marsh  study 
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(Davis,  1981)  involved  a  managed  wetland  area  in  the  Kissimmee 
floodplain.  A  flow-through  marsh  received  canal  inflow  which  averaged 
0.043  mg/L  P,  and  rainfall  averaged  0.174  mg/L  P.  Outflow  averaged 
0.020  mg/L,  with  most  of  the  reduction  taking  place  in  the  upper  one 
third  of  the  flow  path. 

Nitrogen  concentrations  in  the  Lower  Kissimmee  tributaries  were 
found  to  range  from  1.20  to  1.85  mg/L,  mean  wet  season  flow  (Federico, 
1982).  As  with  P,  the  lower  values  were  observed  in  the  less  developed 
watersheds.  In  the  Boney  Marsh  study,  outflow  contained  1.34  mg/L  N. 
Based  on  these  data,  background  concentrations  used  in  equation  5-1  are 
0.02  mg/L  for  P,  and  1.40  mg/L  for  N. 
Rate  coefficients 

Overland  flow.  Overland  flow  has  been  found  to  be  very  effective 
in  removing  pollutants  from  runoff  (Abernathy  et  al . ,  1985).  Most 
studies  have  examined  the  performance  of  filter  strips  for  purifying 
wastewater  discharge,  or  for  removing  nutrients  and  sediment  from  storm 
runoff.  Thompson  et  al.  (1979)  found  greater  than  80%  P,  and  100%  N 
concentration  reductions  after  runoff  from  a  manured  plot  flowed  over  a 
36m  filter  strip.  Draper  et  al .  (1979)  summarize  studies  which  found 
that  vegetative  filters  reduce  P  to  background  levels  in  30  to  1000 
feet.  The  main  factors  related  to  nutrient  removal  were  the  ratio 
of  runoff  volume  to  treatment  area  (i.e.  loading  rate),  and  detention 
time.  In  their  model  (for  the  Ontario  Great  Lakes  basin),  Draper  et 
al.  (1979)  assumed  that  complete  nutrient  attenuation  would  occur 
within  a  distance  of  400  feet.  Due  to  the  low  sorption  capacity  of 
sands,  these  P  removal  rates  from  overland  flow  would  not  be  expected 
on  flatwoods  soils.  Wolfe  (1985)  found  negligible  P  reductions  in 
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overland  flow  of  wastewater  on  a  sandy  soil  with  pasture  cover  in 
Florida. 

With  the  grid  scheme  used  in  this  model  and  the  low  channel  den- 
sity defined  in  the  database,  it  appears  that  there  are  long  distances 
of  overland  flow.  However,  attenuating  the  load  predictions  from  cells 
away  from  a  defined  stream  cell  is  not  warranted  for  two  reasons. 
First,  very  little  attenuation  would  occur  on  areas  of  similar  land  use, 
and  most  of  the  agricultural  land  use  is  pasture.  Second,  the  more 
intense  land  uses  (dairies,  row  crops)  require  drainage.  So,  for  the 
more  critical  areas  where  overland  flow  attenuation  may  be  applicable, 
very  little  likely  occurs. 

Nutrient  loading  in  wetlands  is  attenuated  to  the  nearest  stream. 
In  Table  5-1,  coefficients  from  Kadlec  and  Hammer  (1982)  and  Davis 
(1981)  were  from  marshes  where  a  succession  in  the  plant  community 
resulted  in  a  front  of  rapid  nutrient  uptake.  Reasonable  values  might 
be  expected  to  fall  between  those  of  Davis  (1981)  and  rates  occurring 
in  ditches  and  streams.  Coefficient  values  used  in  BASIN  are  shown  in 
Table  5-2. 

Attenuation  to  sub-basin  outlet.  Attenuation  in  flow  through 
wetlands  is  applied  only  to  the  load  that  is  a  direct  result  of  cattle 
in  the  wetland.  Runoff  from  pastures  is  not  attenuated  through 
wetlands  since  flow  paths  are  not  known.  Also,  the  areas  with  the 
highest  nutrient  loads  in  runoff  will  be  ditched,  thus  will  not  produce 
sheet  flow  through  a  receiving  wetland.  However,  some  attenuation 
would  be  expected  in  the  field  ditches  and  drainageways  leading  to  the 
main  channel,  thus  nutrient  loads  are  attenuated  to  the  outlet  of  the 
sub-watersheds.  Coefficient  values  are  expected  to  be  in  the  range  of 


Ill 


Table  5-1.  Rate  coefficients  for  nutrient  uptake. 


System 


Nitrogen    Phosphorus   Reference 


-1. 


Natural  wetland     0.02-0.06     0.06-0.09    Kadlec  and  Hammer  (1982] 
(Michigan) 


Natural  wetland 
(Boney  Marsh) 

South  Florida 
canals 


0.000025 


Common  values  ford  0.000038 
natural  processes 


0.001-0.006   Davis  (1981) 
0.000013    Mierau  (1981) 
0.000038    Huber  et  al.  (1976] 


a  Rate  of  0.1  day"  assuming  a  velocity  of  0.1  ft/s  (Mierau,  1981) 


Table  5-2.  Rate  coefficients  for  nutrient  uptake  used  in  BASIN. 


System 


RATE  COEFFICIENT 
Units     Nitrogen    Phosphorus 


Flow  through  wetlands       m 
Flow  to  sub-watershed  outlet  m 


Main  channel 
Impoundment 


day 


0.0005 

0.0002 

0.00005 

0.00002 

0.00005 

0.0 

0.08 

0.04 
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those  observed  by  Mierau  (1981)  but  could  be  slightly  higher  due  to 
lower  flow  velocities  and  increased  contact  with  sediment  and  wetland 
areas  in  the  upper  reaches  of  the  watershed  drainage  network. 

Main  channel .  On  a  long  term  basis,  no  permanent  uptake  of 
phosphorus  is  expected  in  the  main  Taylor  Creek  channel.  The  sandy 
soils  have  little  capacity  for  adsorbing  phosphorus,  and  it  is  assumed 
that  the  channel  system  is  effectively  saturated  with  respect  to 
phosphorus.  A  reduction  of  phosphorus  load  may  occur  in  the  main 
channel  in  dry  years,  but  limited  data  suggest  that  in  wet  years  this 
temporarily  stored  phosphorus  is  flushed  from  the  channel  (Ritter  and 
Allen,  1982).  However,  reduction  of  nitrogen  loads  can  be  expected  due 
to  the  denitrification  process.  Coefficients  used  in  BASIN  are  shown 
in  Table  5-2. 

Impoundments.  As  discussed  in  the  literature  review,  the  uptake 
of  nutrients  in  impoundments  can  vary  tremendously  depending  on  the 
type  of  system,  type  of  sediment,  level  of  management,  harvesting  of 
macrophytes,  nutrient  concentrations  in  the  inflow,  and  residence  time. 
Thus  it  is  difficult  to  specify  coefficients  which  are  supposed  to  be 
representative  of  all  BMP  impoundments  applied  in  the  TCNS  basin.  The 
"typical"  impoundment  is  assumed  to  be  on  sandy  soils  (thus  low  P  up- 
take), and  have  a  low  level  of  management  with  no  effective  harvesting. 
Coefficients  chosen  are  in  the  range  of  those  found  by  Reddy  et  al. 
(1982b, a)  for  flooded  fields  (N:  0.03-0.20  day"1;  P:  0.02-0.12  day"  ). 
Distance  estimates 


Estimates  of  flow  distance  are  needed  for  computing  the  nutrient 
attenuation  in  flow  through  wetlands  and  streams.  Distance  from  each 
cell  to  the  nearest  identified  stream  is  included  in  one  of  the  input 


113 

data  files  ("BMP  file")  as  will  be  discussed  later.  This  distance  is 
given  as  the  number  of  cells  away  from  the  stream.  A  distance  of  300m 
is  used  for  cell  to  cell  flow  distances.  This  is  a  compromise  between 
the  perpendicular  and  diagonal  distances  across  a  cell. 

The  direct  distance  from  a  cell  to  the  watershed  outlet  is  deter- 
mined using  the  row/column  coordinates  of  the  cell  and  outlet.  This  is 
multiplied  by  the  300m  per  cell  distance  and  an  arbitrary  factor  of  1.5 
to  help  account  for  the  fact  that  flow  paths  are  not  direct.  The  dis- 
tance from  watershed  outlets  to  the  basin  outlet  is  computed  similarly, 
except  a  factor  of  1.2  is  used  in  this  case. 

Runoff  Detention  by  Impoundments 

CREAMS-WT  does  not  consider  retention  or  impoundment  effects  on 
nutrient  loading  in  runoff,  so  this  was  included  in  the  BASIN  model. 
It  is  assumed  that  detention  basins  are  designed  so  that  after  filling, 
additional  runoff  flows  around  the  impoundment  (Bottcher  and  Baldwin, 
1983).  With  outflow  restricted,  a  slow  drain-down  then  occurs.  To 
simplify  inclusion  in  the  BASIN  model  and  to  minimize  computation,  an 
explicit  relationship  was  sought  between  annual  runoff  volume,  Q, 
detention  basin  volume,  S,  and  the  annual  volume  of  water  detained,  DV. 

A  program  was  written  to  take  daily  runoff  amounts  from  20  years 
of  CREAMS-WT  output,  and  simulate  a  ten  day  exponential  drain-down  to 
determine  the  total  volume  that  is  detained  by  the  impoundment.  This 
was  done  for  five  detention  storage  volumes  for  each  of  eleven  differ- 
ent hydrologic  conditions,  yielding  (20-year)  average  annual  values  of 
DV  for  each  Q  and  S  combination.  The  following  equation  was  fit  to 
these  data  by  non-linear  regression  (coeff.  of  determination  0.997). 
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DY  =  q0-994  *  (i  _  3-0.563*5^  ^5_2^ 

This  relationship  is  used  in  BASIN  to  determine  average  annual 
volume  detained  given  the  design  storage  and  the  average  annual  runoff 
predicted  by  CREAMS-WT.  Average  detention  time  is  found  as  the  cen- 
troid  of  the  area  under  the  drawdown  curve  and  is  approximately  three 
days  for  a  ten  day  draw  down  period.  Nutrient  loads  in  the  detained 
volume  were  then  attenuated  as  in  equation  5-1,  with  the  variable  D  in 
equation  5-1  reflective  of  detention  time  rather  than  distance,  and 
other  parameters  as  previously  defined. 

Direct  Loads  to  Wetlands  and  Streams 

Cattle  with  free  access  to  wetlands  may  spend  6  to  8  hours  a  day 
in  them  to  escape  the  heat.  With  fencing,  the  portion  of  animal  waste 
that  would  have  been  added  directly  to  the  wetland  is  instead  deposited 
on  the  surrounding  pasture.  Nutrient  yield  from  the  pasture,  both 
before  and  after  fencing,  is  handled  by  CREAMS-WT.  However,  CREAMS-WT 
is  not  designed  for  simulating  wetland  conditions,  so  an  estimate  of 
nutrient  yield  from  wetlands  for  the  no-fencing  condition  was  needed. 

Gross  nutrient  load  in  a  wetland  cell  or  stream  is  the  number  of 
animals  in  the  wetland  cell  times  the  nutrient  content  of  an  animal's 
waste.  The  number  of  animals  is  a  function  of  the  animal  density  in 
the  surrounding  pasture,  the  estimated  fraction  of  time  spent  in  the 
wetland,  and  the  degree  of  fencing.  The  number  of  "full  time  equi- 
valent" (FTE)  animals  in  a  wetland  cell  or  stream  is  specified  as  a 
model  input  in  the  "BMP  file"  as  discussed  in  a  following  section. 
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Table  5-3.  Annual  nutrient  content  of  cattle  waste  (per  1000  lb  live 
weight) ,  by  land  use. 


Land  use  Nitrogen    Phosphorus 


(kg/yr)  

Intensively  managed  beef  pasture 

and  dairy  pasture  70.0  13.0 

Improved  pasture  60.0  11.5 

Semi-improved  pasture  50.0  10.0 

Sources:  ASAE  (1983);  Personal  communication  (1984),  Dr.  B.  Kunkle, 

Assoc.  Prof,  Animal  Science  Dept,  Univ.  of  Fla. 


Nutrient  content  in  cattle  waste  is  a  function  of  the  type  of 
animal  and  the  type  of  feed  (ASAE,  1983).  Three  different  nutrient 
contents  were  used,  depending  on  the  land  use  of  the  surrounding 
pasture  (Table  5-3).  The  highest  values  were  for  dairy  cows,  and  the 
lowest  for  beef  cattle  on  semi-improved  pasture. 

All  of  the  animal  waste  excreted  in  a  wetland  or  stream  channel 
would  not  be  expected  to  be  mobile.  The  mobile  nutrients  would  depend 
on  the  fraction  of  readily  soluble  nutrients  in  the  waste,  and  the 
degree  to  which  the  cattle  mix  or  stir  up  the  bottom  sediments  includ- 
ing the  manure.  As  a  means  of  estimating  this  mobile  fraction,  it  was 
assumed  that  a  lower  limit  would  be  the  fraction  of  nutrients  coming 
off  pasture  land  as  predicted  by  CREAMS-WT.  Percentages  of  the  animal 
waste  applied  to  the  land  and  lost  in  the  runoff  for  semi-improved  pas- 
ture with  a  high  runoff  condition  (D  soil  group)  were  approximately  15% 
for  both  N  and  P.  For  wetlands,  then,  the  mobile  fraction  of  animal 
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waste  is  assumed  to  be  20%  of  N  and  P,  which  is  then  attenuated  to 
edge-of-stream.  The  mobile  fraction  in  streams  would  be  expected  to  be 
higher,  and  was  arbitrarily  set  at  50%. 

Background  Loading  Rates 

A  background  nutrient  load  was  simulated  for  all  non-agricultural 
cells.  This  primarily  includes  land  designated  as  urban,  wetlands, 
forests,  and  native  range.  Together  the  non-agricultural  land  use 
covers  about  30%  of  the  TCNS  basins.  This  is  a  large  enough  area  that 
the  natural  background  nutrient  yield  from  these  areas  could  comprise  a 
significant  portion  of  the  total  basin  nutrient  load.  Although  includ- 
ing this  load  will  not  affect  the  absolute  difference  in  loading 
between  various  BMP  implementations,  it  will  affect  the  percentage  dif- 
ference in  water  quality  improvement  between  BMP  scenarios.  The  annual 
N  and  P  loads  used  for  these  background  loads  were  taken  from  data  on 
several  undeveloped  sub-watersheds  in  the  Kissimmee  Basin  as  discussed 
previously  in  the  nutrient  attenuation  section. 

Implementation 

The  flow  of  data  in  the  implementation  of  the  BASIN  model  is  shown 
in  Figure  5-2.  There  are  four  files  which  supply  data  to  BASIN,  and 
these  will  be  described  here. 
Cell  Database 

This  is  the  cell  data  assembled  by  the  Corps  of  Engineers.  Data 
used  by  BASIN  includes:  cell  position  coordinates,  basin  and  sub-basin 
codes,  hydrologic  soil  group  code,  land  use  code,  and  a  code  indicating 
whether  or  not  a  stream  channel  passes  through  the  cell. 
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Figure  5-2.  Data  flow  in  the  implementation  of  the 
BASIN  model. 
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Cell  BMP  File 

This  file  contains  a  record  for  each  cell.  Data  used  by  the  BASIN 
model  includes:  code  for  BMP(s)  applied,  distance  to  the  nearest 
stream  (number  of  cells),  cattle  density  (cattle/ha),  number  of  cattle 
in  stream,  length  of  fencing  (km/km  ),  impoundment  volume  (ha-cm/ha), 
and  field  detention  volume  (ha-cm/ha).  While  this  file  could  be 
created  by  hand,  it  would  be  tedious  for  anything  other  than  a  small 
simulation  area,  so  the  BMP  application  model  developed  by  Bottcher  et 
al.  (1984)  was  used.  This  program  creates  the  BMP  file  based  on  user 
specified  selection  criteria  in  conjunction  with  the  cell  data  from  the 
primary  database. 

BMPs  are  assigned  to  individual  cells  by  considering  the  land  use, 
cattle  density  (for  pasture  land),  hydrologic  soil  group,  and  the 
distance  of  the  cell  from  a  stream  or  wetland.  Cells  neighboring  a 
stream  cell  would  be  expected  to  deliver  greater  amounts  of  N  and  P  to 
streams  than  those  not  bordering  the  stream.  Thus  BMP  application 
levels  are  decreased  with  increasing  distance  of  a  cell  from  a  stream 
or  wetland.  For  fencing  of  pasture  from  wetlands,  decreasing  amounts 
of  fencing  are  applied  as  the  distance  of  the  adjoining  wetland  from  a 
stream  cell  increases. 

All  land  under  dairy  ownership  was  identified  with  the  same  land 
use  code  (18),  but  may  include  a  variety  of  uses  and  cattle  densities. 
For  cells  identified  as  'dairy1,  a  method  for  estimating  cattle  density 
was  needed.  The  location  and  number  of  milking  cows  at  each  dairy  was 
known,  but  additional  breakdown  of  the  'dairy'  land  into  pasture  or 
hayland  was  not  available.  Dairy  cells  were  divided  into  three  groups 
of  decreasing  cattle  density  ('cow  pasture',  'dry  cow  pasture',  and 
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'hayland')  based  on  distance  from  the  barn  cell.  Cattle  density  was 

then  computed  as  a  function  of  the  assigned  dairy  land  use  group  and 

the  size  of  the  milking  herd.  For  purposes  of  categorizing  the 

intensity  of  land  use,  these  three  dairy  land  uses  were  identified, 

respectively,  with  intensively  managed,  improved,  and  semi-improved 

pasture. 

BASIN  Parameter  File 

The  parameter  file  contains  additional  data  required  by  the  model. 
This  includes  the  array  of  discrete  cattle  densities  used  in  the 
CREAMS-WT  simulations,  the  cell  coordinates  of  the  basin  outlet,  and 
the  cell  coordinates  of  the  different  sub-basin  outlets.  The  parameter 
file  used  for  the  TCNS  simulation  is  included  in  Appendix  D. 
CREAMS-WT  Cell  Loads 

The  long  term  average  annual  nutrient  loads  for  different  cells 
are  contained  in  this  file.  Data  in  the  file  are:  cell  access  code, 
annual  nitrogen  load  (kg/ha),  annual  phosphorus  load  (kg/ha),  and 
annual  runoff  (cm).  Putting  this  file  together  is  the  biggest  task  in 
preparing  to  run  BASIN.  The  procedure  used  in  developing  this  file  for 
the  TCNS  basin  will  be  described,  with  the  resulting  file  included  in 
Appendix  D. 

Because  several  soil  associations  and  land  uses  predominate  in  the 
TCNS  basin,  many  of  the  cells  would  have  similar  CREAMS-WT  parameters. 
As  a  result,  a  great  savings  in  computation  time  and  storage  space  can 
be  realized  by  defining  an  access  code  for  each  unique  combination  of 
cell  conditions.  A  single  CREAMS-WT  simulation  can  then  be  run,  and 
the  results  accessed  by  BASIN  for  all  cells  having  similar 
characteristics. 
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Access  code.  The  cell  access  code  for  CREAMS-WT  simulations  is  a 
five  digit  number.  The  first  two  digits  are  the  land  use  code.  The 
third  digit  represents  the  hydrologic  soil  group,  with  the  numbers  1-4 
corresponding  to  A-D.  The  last  two  digits  are  a  code  for  the  cattle 
density,  and  are  really  an  index  into  the  "cattle  density"  array  read 
in  from  the  parameter  file. 

Hydrology  simulations.  The  hydrologic  soil  group  for  flatwoods 
soils  is  frequently  given  a  dual  classification  (A/D,  B/D) ,  the  actual 
classification  being  determined  by  the  drainage  improvements  at  a  site. 
These  dual  classifications  were  converted  to  single  values  by  taking 
the  higher  classification  (A,B)  for  well  drained  land  uses  (row  crops, 
citrus,  improved  pastures).  Other  land  uses  were  converted  to  the 
lower  classification. 

Information  for  defining  hydrologic  parameters  for  CREAMS-WT  is 
limited  to  knowledge  of  the  land  use  and  the  hydrologic  soil  group. 
Accordingly,  nine  different  hydrologic  groups  were  identified  from 
combinations  of  those  two  characteristics  found  in  the  cell  database. 
The  hydrology  parameter  files  are  listed  in  Appendix  D.  Twenty  years 
(1960-1979)  of  daily  rainfall  from  the  Okeechobee  Hurricane  Gate  6 
station  were  used  for  the  precipitation  input.  As  discussed  earlier, 
this  period  is  a  good  representation  of  the  long  term  record  in  both 
mean  and  distribution  of  annual  precipitation  amounts. 

Nutrient  simulations.  Differences  between  nutrient  simulations 
are  based  on  land  use  and  cattle  density.  Fertilization  is  specified 
by  land  use,  in  accordance  with  knowledge  of  typical  practices.  The 
other  major  nutrient  source,  animal  waste,  is  a  function  of  land  use 
(dairy  or  beef)  and  cattle  density.  The  nutrient  loads  used  for  the 
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CREAMS-WT  simulations  are  given  in  Table  5-4.  Animal  waste  loads  are 
multiplied  by  the  cattle  density  being  simulated,  and  are  distributed 
throughout  the  year  in  12  equal  applications.  The  CREAMS-WT  nutrient 
parameter  files  (without  the  fertilizer  parameters  which  are  added  as 
described  earlier  in  this  paragraph)  are  listed  in  Appendix  D.  The 
"fertilizer"  portion  of  the  parameter  files  are  added  to  reflect  the 
land  use  and  the  range  of  cattle  densities  occuring  with  that  land  use. 
Animal  waste  on  beef  and  dairy  pastures  must  be  considered  differ- 
ently because  of  the  source  of  the  nutrients.  Beef  cattle  generally 
receive  no  feed  other  than  grazing,  thus  their  waste  is  largely  a 
recycling  of  nutrients  in  the  pasture  system.  In  contrast,  dairy  cows 
receive  a  large  percentage  of  their  feed  from  sources  outside  the 
pasture.  Thus  while  dairy  waste  can  be  considered  a  "fertilizer" 
application,  it  is  less  reasonable  to  assume  the  same  for  waste  from 
beef  cattle.  In  addition,  as  the  CREAMS  manual  points  out,  all  animal 
waste  applied  to  the  land  is  not  immediately  available.  The  "factor" 
in  Table  5-4  is  to  reflect  the  fraction  of  the  animal  waste  load  that 
is  to  be  simulated  in  CREAMS-WT  as  a  nutrient  load,  taking  into  account 
the  source  of  the  cattle  feed  and  the  availability  of  waste  to  runoff. 

Verification 
Taylor  Creek—Nubbin  Slough 

Combined  flow  and  water  quality  data  in  the  TCNS  basin  are  avail- 
able for  a  six  year  period  (1978-1983)  for  Otter  Creek,  Northwest 
Taylor  Creek,  Williamson  Ditch  and  at  the  basin  outlet  which  is  the 
S-191  structure  at  Lake  Okeechobee.  Records  at  the  Upper  Taylor  Creek 
watershed  outlet  are  also  available  for  this  period,  however  flows  are 
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Table  5-4.  Base 

annua 

1  nutrient  loads  for 

TCNS  simu 

1  at  ion 

LAND  USE 

ANIMAL  WASTE 

FERTILIZER 

Description 

Code 

N 

P 

Factor 

N 

P 

(kg/yr 

/animal 

) 

(kg/ha) 

Beef  Pasture 

Intensive 

15 

70 

13.0 

0.25 

60 

30 

Improved 

16 

60 

11.5 

0.25 

40 

20 

Semi-improv 

ed 

17 

50 

10.0 

0.25 

0 

0 

Dairy  Pasture 

18  * 

Cow  pasture 

(15) 

70 

13.0 

0.75 

0 

0 

Dry-cow  pas 

ture 

(16) 

60 

11.5 

0.50 

40 

20 

Hayland 

(17) 

50 

10.0 

0.25 

40 

20 

*  The  physical  and  management  characteristics  of  these  land 
uses  are  assumed  in  this  breakdown  of  "dairy"  land. 


Table  5-5.  Characteristics  of  watersheds  in  the  Taylor  Creek—Nubbin 
Slough  basin  (area  and  land  use  as  coded  in  cell  database) 


—  CODED 

LAND 

USE 

Name 

Code 

Area 

Beef 

Dairy 

C 

rops 

Wetland 

Other 

(km2) 

(% 

area 



Otter  Creek 

W-13 

22.6 

11 

74 

9 

1 

4 

N.W.  Taylor  Cr 

W-3 

43.3 

55 

29 

1 

11 

5 

Williamson  Ditch 

W-5 

84.9 

51 

3 

3 

7 

36 

Upper  Taylor  Cr 

W-2 

270.0 

40 

24 

3 

9 

25 

Mosquito  Creek 

57.9 

27 

37 

2 

6 

27 

Nubbin  Slough 

56.0 

39 

38 

1 

4 

18 

Henry  Creek 

18.2 

34 

31 

10 

3 

21 

Lettuce  Creek 

62.0 

51 

5 

3 

8 

33 

TCNS  basin 

S-191 

464.0 

39 

25 

3 

8 

25 

(a; 
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Table  5-6.  Observed  annual  loads  and  average  annual  flow  weighted 
concentrations  for  the  Taylor  Creek--Nubbin  Slough  basin 
and  several  sub-watersheds  (1978-1983),  with  corresponding 
predictions  by  the  BASIN  model  for:   (a)  nitrogen,  and 
(b)  phosphorus. 


FLOW  WEIGHTED  N  CONC.  ANNUAL  N  LOAD 

— Observed —   Pred    Observed Pred 

Watershed         Range   Mean   Mean      Range    Mean   Mean 


(mg/1) --(kg/ha) 

Otter  Creek  3.81-8.75  5.86  7.55  8.47-17.64  13.42  21.23 

NW  Taylor  Creek  1.17-2.27  1.75  2.34  0.85-10.88  5.39  6.90 

Williamson  Ditch  1.70-3.07  2.53  2.48  1.90-  9.53  6.95  6.79 

Upper  Taylor  Cr  1.88-2.66  2.14  2.30  1.39-  9.02  5.97  6.54 

TCNS  (S-191)  1.94-3.17  2.57  3.16  4.78-14.78  10.04  8.96 


FLOW  WEIGHTED  P 

CONC. 

ANNUAL  P  LOAD 

— Observed — 

Pred 

Observed 

Pred 

Watershed 

Range   Mean 

Mean 

Range    Mean 

Mean 

i  mn  i  -\  \ 

flen/hal 

ifig/ 1  j- 

Otter  Creek 

1.60-3.54  2.49 

2.27 

2.85-8.75   5.99 

6.39 

NW  Taylor  Creek 

0.45-0.87  0.60 

0.71 

0.26-3.53   1.77 

2.08 

Williamson  D 

'tch 

0.38-1.05  0.59 

0.73 

0.40-2.01   1.51 

1.99 

Upper  Taylor 

Cr 

0.67-1.00  0.86 

0.83 

0.58-4.53   2.42 

2.36 

TCNS  (S-191) 

0.73-1.24  0.91 

1.16 

1.18-5.42   3.64 

3.29 

Sources:  RCWP  (1983,1984);  Unpublished  data,  SFWMD,  West  Palm  Beach. 
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considered  "estimates"  at  this  station.  Water  quality  data  have  been 
collected  at  many  other  locations  within  the  basin  (see  Ritter  and 
Allen,  1982),  but  without  estimates  of  flow.  The  delineation  of  the 
basin,  interior  sub-watersheds  and  stream  channels  as  defined  in  the 
cell  database  are  shown  in  Figure  2-2.  Additional  data  pertaining  to 
these  areas  are  given  in  Table  5-5. 

The  entire  TCNS  basin  was  simulated  with  the  BASIN  model  without 
calibration.  Necessary  data  files  were  developed,  and  parameter  values 
chosen  as  described  in  the  previous  sections.  The  BASIN  parameter  file 
and  CREAMS-WT  cell  loads  file  are  included  in  Appendix  D,  along  with 
the  model  output.  The  initial  simulation  results  are  presented  in  the 
following  discussion.  It  would  have  been  possible  to  do  some  calibra- 
tion after  comparison  with  the  historical  data,  but  this  was  not  done, 
consistent  with  the  model  design  concepts. 

In  line  with  the  discussion  of  the  verification  of  the  CREAMS-WT 
nutrient  model,  flow  weighted  concentrations  are  the  best  data  for 
comparison  with  the  model  output.  BASIN  predicts  long  term  average 
annual  nutrient  yield  so  this  is  compared  with  average  annual  observed 
data  in  Table  5-6. 

Predicted  concentrations  and  loads  compare  very  well  with  observed 
data  for  both  nitrogen  and  phosphorus.  This  is  especially  true  for  an 
uncalibrated  simulation.  Predictions  match  well  with  data  at  the  basin 
outlet  (S-191),  as  well  as  at  the  sub-watershed  outlets.  The  sub- 
watersheds  incorporate  the  range  of  land  use  intensities  that  are  found 
in  the  basin,  from  the  very  heavy  dairy  area  of  Otter  Creek,  to  the 
beef  pastures  of  Williamson  Ditch.  The  differences  in  land  use  between 
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the  watersheds  are  evident  in  the  observed  nutrient  loads  at  the  out- 
lets, and  those  variations  are  reflected  by  the  model. 

Some  secondary  verification  of  model  results  can  be  made  using 
data  from  the  lower  watersheds  in  the  TCNS  basin.  Only  time  series 
data  are  available,  but  this  does  supply  some  indication  of  the  range 
of  nutrient  concentrations  coming  from  these  areas.  The  average  time 
series  concentrations  over  the  period  of  record  are  compared  with  the 
predicted  flow  weighted  concentrations  in  Table  5-7.  Some  sites  and 
data  from  Table  5-6  are  included  to  provide  some  comparison  between 
mean  observed  time  series  and  flow  weighted  concentrations.  Predic- 
tions for  Nubbin  Slough,  Mosquito  and  Henry  Creeks  compare  favorably. 
The  large  difference  in  the  Lettuce  Creek  values  suggest  that  land  use 
definition  in  this  watershed  may  not  be  correct.  The  short  period  of 
record,  and  distinct  problems  with  identifying  the  watershed  boundaries 
(RCWP,  1984)  may  also  contribute  to  this  discrepancy. 

The  appropriateness  of  the  response  of  the  BASIN  model  is  thus 
verified  in  the  prediction  of  the  response  at  the  basin  outlet  and  at 
internal  points  in  the  basin.  This  comparison  covers  the  range  of 
expected  land  usage  for  different  watershed  sizes,  and  was  made  without 
calibration. 

Otter  Creek—Annual  Response 

BASIN  was  designed  with  the  assumption  of  predicting  long  term 
average  response.  As  a  result,  predictions  of  year  to  year  response 
will  be  biased.  Conceptually,  predicted  annual  nutrient  loads  would  be 
expected  to  vary  less  about  the  long  term  mean  than  would  the  corres- 
ponding measured  values.  This  follows  the  assumption  that  during  dry 
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years  a  greater  portion  of  the  nutrient  load  will  be  stored  in  the 
stream  channels  but  will  be  flushed  out  during  wet  years,  a  process  not 
modeled  in  BASIN. 

While  not  directly  within  the  intentions  of  the  model's  use,  the 
annual  predictions  for  Otter  Creek  are  compared  with  observed  values  to 
demonstrate  the  annual  variation  that  does  occur  in  model  predictions. 
Nine  years  of  precipitation  records  (1976-1984)  for  the  Otter  Creek 
watershed  (Thiessen  weighted  data  from  USDA-ARS,  Tifton,  Georgia)  cover 
the  period  of  measured  water  quality  data,  and  provide  a  period  for 
equilibration  of  the  model.  CREAMS-WT  hydrology  and  nutrient  simula- 
tions were  run  for  this  nine  year  period  for  the  different  cell  combi- 
nations in  the  database.  From  these  simulations,  files  of  CREAMS-WT 
cell  loads  were  assembled  for  the  years  1978-1983  (the  only  years  for 
which  observed  data  were  available).  BASIN  was  then  run  for  each  of 
these  years  without  calibration.  The  predicted  annual  flow  weighted 
concentrations  are  shown  in  Table  5-8  with  observed  data. 

Predicted  nitrogen  concentrations  are  high  (as  in  Table  5-6)  but 
do  follow  the  trends  in  the  historical  data.  The  first  two  years  of 
phosphorus  predictions  appear  to  be  affected  by  the  significant  errors 
in  the  runoff  prediction.  (A  check  of  the  precipitation  record  showed 
annual  differences  of  up  to  16cm  between  SFWMD  and  ARS  data).  For  the 
final  four  years,  predicted  phosphorus  concentrations  matched  observed 
values  and  followed  annual  variations. 
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Table  5-7.  Observed  mean  time  series  concentrations  (bi-weekly 
sampling)  with  mean  flow  weighted  concentrations 
predicted  by  BASIN. 


NITROGEN 

PHOSPHORUS 

Period  of 

Watershed 

Observed 

Pred 

Observed 

Pred 

Record 

(mg/1) 

Otter  Creek 

6.38 

7.55 

2.43 

2.27 

78-83 

N.W.  Taylor  Creek 

1.50 

2.34 

0.43 

0.71 

" 

Williamson  Ditch 

2.18 

2.48 

0.44 

0.73 

" 

TCNS  basin 

2.88 

3.16 

0.97 

1.16 

Mosquito  Creek 

6.23 

5.92 

2.35 

1.65 

H 

Nubbin  Slough 

8.42 

7.44 

2.33 

2.02 

11 

Henry  Creek 

6.48 

7.58 

2.04 

1.92 

80-83 

Lettuce  Creek 

1.73 

4.32 

0.26 

1.15 

Sources:  RCWP  (1983,1984);  Unpublished  data,  SFWMD,  West  Palm  Beach, 


Table  5-8.  Observed  and  predicted  runoff  and  mean  annual  flow 
weighted  concentrations  for  Otter  Creek  (W-13). 


YEAR 

RUNOFF 

NITROGEN 

PHOSPHORUS 

Obsv   Pred 

Obsv 

Pred 

Obsv 

Pred 

(cm) 

i 

mn/n 

l'"y  ' ) 

1978 

17.08  26.98 

8.75 

7.59 

3.54 

2.51 

1979 

28.16  66.55 

5.69 

4.64 

3.03 

1.60 

1980 

11.38  14.15 

7.44 

10.97 

2.51 

3.24 

1981 

21.57  15.24 

5.11 

11.28 

2.09 

2.73 

1982 

32.65  45.60 

3.81 

5.36 

1.60 

1.75 

1983 

40.44  34.17 

4.36 

7.56 

2.16 

2.56 

Sources:  RCWP  (1983,1984);  Unpublished  data,  SFWMD,  West  Palm  Beach. 
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Table  5-9.  Predicted  loads  for  TCNS  basin  for  base  simulation  and  for 
simulation  with  a  high  rate  BMP  scenario,  with  percent 
reductions  attributed  to  BMPs. 


NITROGEN  PHOSPHORUS  

Base     BMP    %       Base     BMP    % 
Load    Load  Reduc     Load    Load  Reduc 


— (kg/h 

a)  — 

% 

---(kg/h 

a)  — 

% 

Otter  Creek 

21.23 

9.12 

57 

6.39 

3.40 

47 

N.W.  Taylor  Cr 

6.90 

5.51 

20 

2.08 

1.66 

20 

Williamson  Ditch 

6.79 

4.51 

34 

1.99 

1.40 

30 

Upper  Taylor  Cr 

6.54 

4.77 

27 

2.36 

1.66 

30 

Mosquito  Creek 

17.00 

5.43 

68 

4.75 

1.65 

65 

Nubbin  Slough 

20.91 

5.79 

72 

5.68 

1.82 

68 

Henry  Creek 

21.07 

6.87 

67 

5.33 

2.35 

56 

Lettuce  Creek 

12.07 

4.77 

60 

3.21 

1.38 

57 

TCNS  (S-191) 

9.04 

4.67 

48 

3.29 

1.67 

49 
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BMP  Simulation 

Data  for  verification  of  the  accuracy  of  the  model  in  simulating 
BMPs  are  not  available.  However,  a  simulation  was  run  using  the 
"maximum"  BMP  scenario  considered  by  Bottcher  et  al.  (1984)  to  provide 
some  indication  of  how  the  model  responds.  The  "maximum"  scenario 

applies  high  rates  of  BMPs  to  all  applicable  cells.  This  includes: 

2  2 

fencing  of  streams  (0.099  km/km  ),  fencing  of  wetlands  (0.04  km/km  ), 

impoundment  (5.1  ha-cm/ha),  and  lot  retention  (2.5  ha-cm/ha).  The 
model  output  for  the  BMP  simulation  is  included  in  Appendix  D.  Table 
5-9  compares  the  base  simulation  (no  BMPs)  to  that  of  the  BMP  simula- 
tion, with  the  percentage  reduction  indicated. 

While  load  reductions  predicted  by  this  BMP  scenario  cannot  be 
verified  by  measured  data,  they  do  appear  conceptually  sound.  Signif- 
icant reductions  in  nutrient  load  would  be  expected  were  such  a  BMP 
program  implemented.  However,  lacking  the  complete  removal  of  cattle 
and  cessation  of  current  land  use,  nutrient  loads  from  the  watershed 
would  be  expected  to  be  above  native  background  levels.  Predicted 
loads  fall  within  this  range  and  appear  reasonable. 

Sensitivity  Analysis 
A  numerical  sensitivity  analysis  of  several  parameters  and  of  the 
input  cell  nutrient  loads  was  performed,  with  results  listed  in  Tables 
5-10  and  5-11.  The  relative  sensitivity  (defined  in  equation  3-17)  was 
calculated  for  two  different  model  results:  average  loads  delivered  to 
the  sub-basin  outlets,  and  to  the  basin  outlet.  For  the  sub-basin 
outlet  sensitivity,  the  loads  at  all  outlets  were  summed  for  a  single, 
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composite  comparison.  Delivery  to  the  basin  outlet  is  directly  given 
in  the  model  output. 

The  difference  between  the  predictions  at  sub-basin  outlets  and 
the  basin  outlet  is  simply  the  channel  nutrient  attenuation.  Since  no 
such  attenuation  is  modeled  for  phosphorus,  there  is  no  difference 
between  relative  sensitivities  for  predictions  at  sub-basin  and  basin 
outlets.  For  nitrogen,  which  is  attenuated  in  delivery  to  the  basin 
outlet,  model  sensitivity  at  the  basin  outlet  is  proportionally  reduced 
for  parameters  which  affect  only  the  yield  at  the  sub-basin  outlet 
(cell  loads,  BCON) .  Nitrogen  model  parameters  which  apply  to  channel 
attenuation  (KCHN,  BCDN)  have  increased  sensitivity. 

Model  results  are  sensitive  to  the  cell  nutrient  loads,  more  so 
for  phosphorus  than  nitrogen.  For  phosphorus,  cell  nutrient  loads  are 
clearly  the  most  important  (of  the  inputs  evaluated)  in  determining 
predicted  results.  The  rate  coefficient  for  stream  attenuation  to  the 
sub-basin  outlet  (BCOP)  is  moderately  sensitive.  Interestingly,  the 
background  concentration  for  phosphorus  (BCDP)  has  little  effect  on 
model  results.  The  nitrogen  predictions  are  affected  at  the  outlet  as 
much  by  rate  coefficients  and  the  background  concentration  as  by  the 
cell  loads. 

This  sensitivity  analysis  is  appropriate  only  to  this  basin  as  the 
model  predictions  are  also  a  function  of  distances  to  the  outlets,  the 
spatial  relationship  of  the  watersheds  to  the  basin  outlet,  basic 
nutrient  concentrations,  and  other  characteristics  of  the  basin.  For 
this  basin,  it  does  indicate  that  for  phosphorus,  effort  should  be 
spent  in  improving  the  cell  load  predictions,  while  for  nitrogen,  para- 
meters in  the  BASIN  model  are  equally  as  important. 
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Table  5-10.  Sensitivity  of  predicted  nitrogen  loads  at  sub-basin 
and  basin  outlets  to  selected  model  parameters. 


SUB-BASIN  OUTLET     BASIN  OUTLET 
Parameter   Fraction    Load  Relative    Load  Relative 
(base)     of  base    Fract    Sens     Fract    Sens 


CREAMS-WT  cell  loads 

0.50       0.84    0.31  0.90  0.20 

0.75       0.92    0.31  0.95  0.20 

1.25       1.08    0.31  1.05  0.20 

1.50       1.16    0.31  1.10  0.20 

Rate  coeff.  for  stream  attenuation  to  sub-basin  outlet 

0.10       1.32   -0.36  1.29  -0.32 

BCON        0.50       1.16   -0.32  1.14  -0.29 

(.015)      2.00       0.79   -0.21  0.81  -0.19 

10.00       0.41   -0.07  0.49  -0.06 

Rate  coeff.  for  main  channel  attenuation 


0.10       1.00 

KCHN        0.50       1.00 

(.015)      2.00       1.00 

10.00       1.00 

0 
0 
0 
0 

Background  concentration 

0.50       0.95 

BCON        0.75       0.98 

(1.4)       1.25       1.02 

1.50       1.05 

0.10 
0.10 
0.10 
0.10 

1.26 

-0.29 

1.11 

-0.22 

0.88 

-0.12 

0.71 

-0.03 

0.86 

0.29 

0.93 

0.29 

1.07 

0.29 

1.14 

0.29 
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Table  5-11.  Sensitivity  of  predicted  phosphorus  loads  at  sub- 
basin  and  basin  outlets  to  selected  model  parameters, 


SUB-BASIN  OUTLET     BASIN  OUTLET 
Parameter  Fraction     Load  Relative    Load  Relative 
(base)     of  base     Fract    Sens     Fract    Sens 


CREAMS-WT  cell  loads 

0.50  0.76  0.47  0.76  0.47 

0.75  0.88  0.47  0.88  0.47 

1.25  1.12  0.47  1.12  0.47 

1.50  1.24  0.47  1.24  0.47 

Rate  coeff.  for  stream  attenuation  to  sub-basin  outlet 

0.10  1.15  -0.17  1.16  -0.17 

BC0P        0.50  1.08  -0.16  1.08  -0.16 

(.006)      2.00  0.86  -0.14  0.86  -0.14 

10.00  0.40  -0.07  0.40  -0.07 

Rate  coeff.  for  main  channel  attenuation  * 

KCHP     value=  0.001  1  0  0.96  * 

(0.0)     value=  0.01  1  0  0.74 

Background  concentration 

0.10  0.98  0.02  0.98  0.02 

BCDP        0.50  0.99  0.02  0.99  0.02 

(.2)        2.00  1.02  0.02  1.02  0.02 


Cannot  be  computed  since  base  parameter  value  was  0. 


CHAPTER  VI 
SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


Summary 
The  South  Florida  Flatwoods,  is  a  very  flat  region,  with  sandy 
soils  and  a  high  water  table  that  fluctuates  from  the  surface  to  two 
meters  deep.  Several  characteristics  of  this  region  are  rather  unique 
compared  to  the  concepts  followed  in  the  design  of  most  water  quality 
models. 

1.  The  hydrologic  response  is  dominated  by  the  flat  terrain, 
sandy  soils,  and  high  water  table.  Runoff  occurs,  not  by  the 
generation  of  an  infiltration  limited  "rainfall  excess",  but 
when  the  available  storage  in  the  profile  (depth  to  the  water 
table)  is  filled  and  the  water  table  rises  close  to  the 
surface. 

2.  The  phosphorus  buffering  capacity  of  the  sandy  soils  is  very 
low  and  phosphorus  will  leach  through  the  profile.  The  phos- 
phorus content  of  the  soil  may  build  up  in  association  with 
liming  and  animal  waste,  but  is  also  readily  released  into 
runoff. 

3.  Wetland  areas  can  be  a  significant  fraction  of  the  land  area, 
and  commonly  border  streams  and  drainageways.  These  wetlands 
receive  runoff  water  and  direct  waste  loading  from  cattle,  and 
have  been  proposed  as  "filtering"  systems  for  removing 
nutrients  from  runoff. 
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4.  Agricultural  land  use  is  predominantly  pasture  for  beef  and 
dairy  cattle.  There  is  very  little  tilled  land.  Management 
practices  for  the  reduction  of  nonpoint  pollution  are  thus 
oriented  towards  pasture  and  cattle  management  instead  of  the 
more  common  emphasis  on  reducing  erosion. 
The  CREAMS-WT  and  BASIN  water  quality  models  are  designed  to  re- 
flect the  hydrologic,  nutrient,  and  management  characteristics  typical 
of  flatwoods  watersheds.  These  models  will  provide  a  means  of  evalua- 
ting the  relative  effects  of  alternative  management  practices  at  both 
the  field  and  basin  level.  Since  management  practices  are  implemented 
on  a  field-by-field  basis,  field-scale  simulation  can  provide  the  best 
assessment  of  the  effectiveness  of  those  practices  in  reducing  nonpoint 
pollution.  However,  an  estimate  of  the  impact  at  the  watershed  outlet 
is  also  valuable  for  evaluating  basin-wide  programs. 

CREAMS-WT  Field  Scale  Model 

CREAMS-WT  consists  of  hydrology  ana  nutrient  sub-models,  and 
predicts  runoff  and  nutrient  yield  (nitrogen  and  phosphorus)  for  field- 
sized  areas.  It  is  a  modified  version  of  CREAMS  (Chemicals,  Runoff, 
and  Erosion  from  Agricultural  Management  Systems  [Knisel,  1980])  that 
is  designed  to  represent  the  characteristics  of  flatwoods  watersheds. 
CREAMS-WT  was  designed  to  1)  simulate  the  major  physical  processes 
that  control  water  balance  and  movement  of  plant  nutrients,  2)  use 
physically-based  parameters  that  can  reflect  changes  in  management 
systems,  3)  operate  on  a  daily  time  step  for  computational  efficiency, 
and  4)  be  field  scale.  Because  the  parameters  of  the  model  are 
physically  based,  they  can  be  estimated  from  site  visits,  maps,  soil 
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survey  reports,  and  the  literature,  with  the  additional  advantage  that 
the  need  for  calibration  is  minimized.  Thus  the  model  can  be  used  for 
simulating  ungaged  areas. 

Hydrology.  Flatwoods  watersheds  respond  as  storage-based  rather 
than  as  infiltration-based  hydrologic  systems.  Components  added  to 
CREAMS-WT  simulate  the  dynamic  water  table  while  limiting  the  amount  of 
deep  seepage,  and  use  a  modified  curve  number  method  for  predicting 
runoff. 

The  water  table  is  simulated  by  limiting  the  amount  of  water 
leaving  the  profile  as  "deep  seepage"  and  using  an  accounting  procedure 
to  keep  track  of  the  water  in  the  profile  as  it  is  withdrawn  by  evapo- 
transpiration  and  replenished  by  rainfall.  A  water  table  recession 
curve  is  used  in  simulating  the  water  budget  below  the  root  zone 
considered  by  CREAMS.  The  height  of  the  water  table  is  important 
because  it  is  the  best  indicator  of  the  volume  of  storage  available  in 
the  profile. 

The  modified  SCS  runoff  equation  used  in  CREAMS  was  adapted  to 
better  represent  the  storage-based  hydrologic  system  of  flatwoods 
watersheds.  CREAMS-WT  uses  an  estimate  of  the  actual  available  storage 
in  the  profile  in  place  of  the  storage  parameter  (which  is  commonly 
determined  by  the  "curve  number")  in  the  basic  SCS  runoff  equation.  A 
theoretical  evaluation  of  the  CREAMS-WT  implementation  of  the  SCS 
runoff  equation  shows  that  the  method  conceptually  represents  the 
rainfall /runoff  process  on  flatwoods  watersheds.  In  addition  to 
representing  the  storage-based  system,  the  algorithm  includes,  in 
simple  form,  the  variable  source  area  concept  of  runoff  generation. 
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Water  table  depth  and  annual  water  balance  predicted  by  the 
CREAMS-WT  hydrology  model  compare  favorably  with  data  from  three  flat- 
woods  watersheds  of  65,  291,  and  1456  ha.  CREAMS-WT  predictions  show 
dramatic  improvement  over  the  results  of  CREAMS  simulations  of  the  same 
watersheds.  A  20-year  simulation  using  precipitation  data  from  the 
same  region  shows  that  the  model  is  stable  and  performs  predictably 
over  the  range  of  natural  rainfall  variability.  Long  term  average 
annual  evapotranspiration  predicted  by  the  model  and  estimated  from  the 
Upper  Taylor  Creek  watershed  also  compare  favorably. 

Nutrients.  The  high  phosphorus  buffering  capacity  of  the  soil 
that  is  assumed  in  CREAMS  is  not  realistic  for  the  sandy,  low  clay  con- 
tent, flatwoods  soils.  Changes  were  made  in  the  nutrient  component  to 
represent  the  lower  buffering  capacity  of  flatwoods  soils.  This 
modification  enables  CREAMS-WT  to  better  represent  observed  phosphorus 
dynamics. 

Data  from  three  small  watersheds  were  available  for  model  evalua- 
tion. These  sites  cover  the  spectrum  of  areal  nutrient  loading  in  the 
TCNS  basin,  from  a  beef  pasture  with  no  fertilization,  to  a  dairy  which 
has  fertilized  pastures  and  a  high  animal  waste  nutrient  load.  This 
range  of  watershed  conditions  was  well  represented  by  the  predicted 
annual  nutrient  loads  from  the  CREAMS-WT  simulations  which  were  run 
without  calibration.  In  addition  to  comparison  of  annual  loads,  pre- 
dicted monthly  flow  weighted  concentrations  also  reflected  trends  in 
observed  data.  Thus,  CREAMS-WT  gives  reasonable  predictions  of  annual 
nutrient  load,  and  represents  the  distribution  of  that  load  throughout 
the  year. 
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BASIN  Model 

BASIN  provides  an  estimate  of  the  water  quality  impact  at  the 
basin  outlet.  Individual  field/cells  are  delineated  within  the  basin, 
with  CREAMS-WT  used  to  predict  average  annual  nutrient  loads  for  each 
cell.  The  BASIN  model  then  predicts  the  delivery  of  cell  nutrient 
loads  to  sub-basin  and  basin  outlets,  and  provides  some  functions  that 
are  beyond  the  range  of  CREAMS-WT.  BASIN  incorporates  the  following 
functions: 

1)  Predict  nutrient  delivery  to  sub-basin  and  basin  outlets  using 
the  simulated  CREAMS-WT  edge-of -stream  loads; 

2)  Model  the  effects  of  BMPs  that  cannot  be  simulated  in  CREAMS-WT 
(fencing,  impoundments); 

3)  Compute  nutrient  loads  due  to  cattle  in  streams  and  wetlands, 
and  for  wetlands,  attenuate  those  loads  in  flow  through  the 
wetlands  to  streams; 

4)  Provide  background  loads  for  forests,  native  range,  and  other 
non-agricultural  land  use  in  the  basin. 

A  first-order  rate  equation  is  used  for  modeling  the  attenuation 
of  wetland  animal  waste  load  in  flow  through  the  wetland  to  the  nearest 
stream,  wetland  edge-of -stream  and  CREAMS-WT  cell  loads  to  sub-basin 
outlets,  and  channel  flow  from  sub-basin  outlets  to  the  basin  outlet. 
Nutrient  attenuation  is  a  function  of  distance  of  flow  and  the  rate 
coefficient.  Rate  coefficients  for  each  process  were  determined  from 
nutrient  uptake  studies  reported  in  the  literature. 

Nutrient  delivery  in  the  TCNS  basin  was  simulated  using  a  5  ha 
grid  cell  as  a  "field".  The  simulation  was  run  without  calibration, 
with  parameters  selected  from  the  literature  and  from  knowledge  of  the 
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area.  Predicted  long  term  average  annual  nitrogen  and  phosphorus 
yields  compare  closely  with  mean  observed  values  at  the  basin  outlet 
and  at  several  gaged  sub-basin  outlets.  Sub-basins  encompass  the  range 
of  land  use  intensities,  from  a  predominantly  beef  range  watershed,  to 
the  high  intensity  Otter  Creek  watershed  with  its  numerous  dairies. 
Thus  the  model  gives  reasonable  predictions  not  only  at  the  basin 
outlet,  but  also  at  internal  points  in  the  basin,  with  the  model's 
response  also  verified  over  the  range  of  land  use  intensities  expected. 

Conclusions 

1.  Flatwoods  watersheds  are  storage-based  hydrologic  systems,  and  the 
SCS  runoff  equation,  which  is  a  storage-based  model,  is  a  suitable 
method  for  predicting  runoff  from  these  watersheds  although  an  ap- 
propriate estimate  of  the  storage  parameter,  S,  must  be  determined. 

2.  The  SCS  runoff  equation  as  modified  for  the  CREAMS-WT  model  is 
particularly  appropriate  for  flatwoods  watersheds,  enabling 
realistic  physical  determination  of  the  storage  parameter,  and 
incorporating  in  simple  form,  the  variable  source  area  concept. 

3.  The  CREAMS-WT  hydrology  model  conceptually  represents  the  major 
components  of  the  hydrologic  system  of  flatwoods  watersheds  and 
can  give  reasonable  predictions  of  annual  water  balance  components 
without  calibration. 

4.  The  CREAMS-WT  nutrient  model  can  give  reasonable  predictions  of 
annual  nitrogen  and  phosphorus  loads  from  flatwoods  watersheds 
without  calibration. 

5.  The  BASIN  model  gives  reasonable  estimates  of  average  annual  loads 
from  the  Taylor  Creek—Nubbin  Slough  basin  and  sub-watersheds,  and 
should  be  applicable  to  other  similar  flatwoods  watersheds. 
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6.  Based  on  the  conceptual  design  of  the  models,  the  use  of 

physically-based  parameters,  and  the  verification  of  the  models 
without  calibration:  The  CREAMS-WT  and  BASIN  models  can  be  used  to 
evaluate  the  relative  effectiveness  of  alternative  management 
practices  at  the  field  level  and  at  the  watershed/basin  outlet. 

Recommendations  for  Additional  Research 
There  are  two  major  areas  that  need  further  research  to  advance 
the  state-of-the-art  in  hydrologic/water  quality  models  for  flatwoods 
watersheds. 

Hydrologic  Data 

Additional  evaluation  of  the  SCS  runoff  equation  (and  associated 
curve  number  method)  is  needed  to  clarify  the  parameter  values  appro- 
priate for  flatwoods  watersheds.  This  effort  should  seek  to  relate 
watershed  characteristics,  land  use,  and  management  practices  to  the 
Smx/ULE  ratio  discussed  in  this  study,  or  alternately,  to  the  commonly 
used  curve  number.  In  a  monitoring  program,  attempts  must  be  made  to 
provide  high  quality  data  by  careful  attention  to  the  sources  of  error 
that  plague  hydrologic  data  collection  on  flatwoods  sites. 

Phosphorus  Model 

A  more  detailed  phosphorus  model  is  needed.  The  inclusion  of  a 
phosphorus  isotherm  relationship  will  not  be  adequate  as  this  does  not 
describe  the  relationship  between  the  equilibrium  phosphorus  concen- 
tration (EPC)  and  the  change  in  the  soil  phosphorus  content.  A  phos- 
phorus model  for  flatwoods  must  include  a  phosphorus  budget,  since 
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interactions  with  soil  and  runoff  vary  with  the  soil  phosphorus 
content.  A  phosphorus  model  must  also  be  able  to  reflect  the  manage- 
ment practices  appropriate  to  these  soils.  This  would  include  plant 
uptake,  for  modeling  the  utilization  of  animal  waste,  and  simulating 
the  effects  of  soil  modifiers,  such  as  lime  or  alum. 

The  phosphorus  model  in  EPIC  (Jones  et  al . ,  1984a, b;  Sharpley  et 
al.,  1984)  is  built  around  a  phosphorus  budget,  modeling  the  pathways 
between  the  different  phosphorus  pools.  This  model  is  one  of  the  most 
complete  phosphorus  models  designed  from  a  water  quality  perspective. 
Its  applicability  to  flatwoods  conditions  should  be  evaluated,  and 
appropriate  parameter  values  determined  and/or  modifications 
incorporated. 
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APPENDIX  A 
RELATIONSHIP  BETWEEN  CURVE  NUMBER  AND  AVAILABLE  STORAGE 


In  the  SCS  runoff  equation,  the  storage  parameter,  S,  is  generally 
determined  through  a  second  parameter,  the  curve  number.  Curve  number 
is  related  to  S  by  the  function, 


1000 

10 

CN 


The  CREAMS  version  of  the  curve  number  method  uses  this  relationship  to 
determine  the  maximum  value  of  S  as, 

1000 
Smx  =  10 

CN, 

where  CN.  is  the  curve  number  at  antecedent  moisture  condition  I  (AMC 
I)  which  corresponds  to  dry  conditions.  An  explicit  function  is  used 
in  CREAMS  to  determine  the  CNr  value  from  the  input  curve  number  which 
is  for  average  moisture  (AMC  II).  This  function  was  curve  fitted  to 
the  corresponding  curve  numbers  for  the  different  antecedent  moisture 
conditions  as  tabulated  by  the  SCS  and  is: 

CNj   =  -16.91  +  1.348(CNn)    -  0.01379(CNn  )2  +  0.0001177(CNn  )3 

For  reference,  values  of  the  input  curve  number  (AMC  II)  and  the 
equivalent  maximum  storage  as  calculated  by  the  CREAMS  algorithms  are 
tabulated  here. 
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Table  A-l.  Curve  number  (AMC  II)  and  equivalent  maximum  storage 
(Smx)  as  calculated  by  CREAMS. 


CN  Smx  (ins 


100  0.24 

99  0.46 

98  0.69 

97  0.93 

96  1.17 

95  1.41 

94  1.67 

93  1.93 

92  2.19 

91  2.46 

90  2.74 

89  3.02 

88  3.31 

87  3.61 

86  3.91 

85  4.22 

84  4.54 

83  4.87 

82  5.20 

81  5.54 


80  5.89 

79  6.25 

78  6.61 

77  6.99 

76  7.38 

75  7.77 

74  8.17 

73  8.59 

72  9.02 

71  9.45 

70  9.90 

69  10.68 

68  10.84 

67  11.32 

66  11.82 


65  12.34 

60  15.16 


APPENDIX  B 
CREAMS-WT  HYDROLOGY  PARAMETER  FILES 


Table  B-l.  Hydrology  parameter  files  for  Bass  West  pasture. 


CREAMS  parameter  file: 


BASS  PASTURE  (WEST),  OKEECHOBEE  COUNTY,  FL 

MODIFIED  HYDROLOGY  SIMULATION  --  WATER  TABLE  VERSION 

SOIL  MOISTURE  STORAGE  FROM  ARS  STORAGE  vs.  WATER  TABLE  CURVE 


81.3   80.7   75.1 
447.   405.   370. 


79152 

0      1 

1 

160.0 

1.00    .25 

.50 

3.3 

.43 

0.03 

.20 

80.   .010 

2.0 

30.0 

.056 

.277   0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2   66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340.   436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.50 

35 

.50 

125 

1.00 

225 

1.00 

335 

.80 

366 

.50 

1 

0     0 

62.5 

59.3   69.4 

71.7 

76.9 

80.8 

82.3 

69.6 

61.5 

1 

0     0 

53.3 

63.7   64.5 

74.1 

76.3 

81.9 

82.5 

67.0 

61.8 

1 

0     0 

60.2 

61.2   66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

-1 

0     0 

CREAMS-WT 

parameters: 

82.4  80.9  76.0 
80.8  80.3  76.9 
81.3   80.7   75.1 


1  'IWTFLAG   -   (II)    =   1,   PRINT  WT;   =  0,   DON'T 

0.000  'DSP        (F5.0)      -   DEEP   SEEPAGE  RATE    (IN/DAY] 

0.080  'PORSUB  -   SUB-SURFACE   POROSITY 

0.36  'RXP   -  EXPONENT   FOR  WT  RECESSION 
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Table  B-2.  Hydrology  parameter  files  for  Armstrong  Slough, 


CREAMS  parameter  file: 


ARMSTRONG  SLOUGH 
MODIFIED  HYDROLOGY  SIMULATION  --  WATER  TABLE  VERSION 

SOIL  MOISTURE  STORAGE  FROM  ARS  STORAGE  vs.  WATER  TABLE  CURVE 
79152     0      1      1 

3250.0   1.00    .25    .50    3.3    .43   0.03 
.20    82.   .010    2.0   30.0 

.056  .277   0.73   0.84   .847   1.00   1.00 

60.2  61.2   66.4   74.2   76.1   80.8   82.7   81.3   80.7   75.1 
69.5  65.1 

282.  340.   436.   510.   507.   472.   452.   447.   405.   370. 

302.  281. 
0.5 

1    .50 
35    .50 

125  1.00 

225  1.00 
335  .80 
366    .50 

1      0     0 

62.5  59.3   69.4   71.7   76.9   80.8   82.3   82.4   80.9   76.0 

69.6  61.5 

1      0     0 

53.3  63.7   64.5   74.1   76.3   81.9   82.5   80.8   80.3   76.9 
67.0  61.8 

1      0     0 

60.2  61.2   66.4   74.2   76.1   80.8   82.7   81.3   80.7   75.1 

69.5  65.1 

-10     0 


CREAMS-WT  parameters: 


1  'IWTFLAG  -  (II)  =  1,  PRINT  WT;  =  0,  DON'T 

0.002  'DSP   (F5.0)   -  DEEP  SEEPAGE  RATE  (IN/DAY] 

0.080  'PORSUB  -  SUB- SURFACE  POROSITY 

0.36  'RXP  -  EXPONENT  FOR  WT  RECESSION 
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Table  B-3.  Hydrology  parameter  files  for  SEZ  Dairy. 


CREAMS  parameter  file: 

SEZ  DAIRY 
MODIFIED  HYDROLOGY  SIMULATION  --  WATER  TABLE  VERSION 
SOIL  MOISTURE  STORAGE  FROM  ARS  STORAGE  vs.  WATER  TABLE  CURVE 


81.3   80.7   75.1 
447.   405.   370. 


79152 

0 

1 

1 

675.0 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

81. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.50 

35 

.50 

125 

1.00 

225 

1.00 

335 

.80 

366 

.50 

1 

0 

0 

62.5 

59.3 

69.4 

71.7 

76.9 

80.8 

82.3 

69.6 

61.5 

1 

0 

0 

53.3 

63.7 

64.5 

74.1 

76.3 

81.9 

82.5 

67.0 

61.8 

1 

0 

0 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

-1 

0 

0 

82.4  80.9  76.0 
80.8  80.3  76.9 
81.3   80.7   75.1 


CREAMS-WT  parameters: 


1  'IWTFLAG  -  (II)  =  1,  PRINT  WT;  =  0,  DON'T 

0.000  'DSP   (F5.0)   -  DEEP  SEEPAGE  RATE  (IN/DAY) 

0.080  'PORSUB  -  SUB- SURF ACE  POROSITY 

0.36  'RXP  -  EXPONENT  FOR  WT  RECESSION 
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Table  B-4.  Hydrology  parameter  files  for  20  year  simulation. 

CREAMS  parameter  file: 


LONG  TERM  SIMULATION  FOR  FLATWOODS  WATERSHED 

MODIFIED  HYDROLOGY  SIMULATION  --  WATER  TABLE  VERSION 

SOIL  MOISTURE  STORAGE  FROM  ARS  STORAGE  vs.  WATER  TABLE  CURVE 


81.3   80.7   75.1 
447.   405.   370. 


60001 

0 

0 

1 

160.0 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

80. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

68.5 

61.8 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.50 

35 

.50 

125 

1.00 

225 

1.00 

335 

.80 

366 

.50 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

CREAMS-WT  parameters: 


0  'IWTFLAG   -    (II)    =  1,   PRINT  WT;   =  0,   DON'T 

0.000  'DSP        (F5.0)      -  DEEP  SEEPAGE  RATE   (IN/DAY' 

0.080  'PORSUB  -   SUB-SURFACE   POROSITY 

0.36  'RXP   -  EXPONENT  FOR  WT  RECESSION 


APPENDIX  C 
CREAMS-WT  NUTRIENT  PARAMETER  FILES 


Table  C-l.  Nutrient  parameter  file  for  Bass  West  Pasture. 


BASS  PASTURE  (WEST),  OKEECHOBEE  COUNTY,  FL 
NUTRIENT  SIMULATION 


150   15.0  .0003  .0001   0.12   0.12   99.0   99.0   99. 


79152 

1 

0 

.43 

0.25 

1.0 

2 

.15 

.150 

15.0 

99.0 

1.0 
79365 

1 

01 

365 

450. 

5000. 

5. 

79305 

67.3 

33.6 
80365 

1.0 

1 

01 

365 

450. 

5000. 

5. 

80305 

67.3 

33.6 
81365 

1.0 

1 

01 

365 

450. 

5000. 

5. 

81305 

67.3 

33.6 
82365 

1.0 

0 

01 

365 

450. 

5000. 
0 

5. 

45.    180    90.    64. 


45.    180    90.    64. 


45.    180    90.    64. 


45.    180    90.    64. 
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Table  C-2.  Nutrient  parameter  file  for  Bass  East  Pasture. 


BASS  PASTURE  (EAST),  OKEECHOBEE  COUNTY,  FL 
NUTRIENT  SIMULATION 

79152     0     0     0      1 
.43   0.25    1.0 

2 
.10   .100   15.0  .0003  .0001   0.10   0.10   99.0   99.0   99, 
99.0    1.0 
79365 
1     01    365 
450.  5000.     5.    45.    180    90.    64. 
79305 

0.1    0.1    1.0 
80365 
1     01    365 
450.  5000.     5.    45.    180    90.    64. 
80305 

0.1    0.1    1.0 
81365 
1     01    365 
450.  5000.     5.    45.    180    90.    64. 
81305 

0.1    0.1    1.0 
82365 
0     01    365 
450.  5000.     5.    45.    180    90.    64. 
0 


Table  C-3.  Nutrient  parameter  file  for  SEZ  Dairy. 
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79152 

0 

0 

.43 

0.25 

1.0 

2 

.20 

.200 

15.0 

99.0 

1.0 
79365 

7 

01 

365 

450. 

5000. 

5. 

79152 

11.4 

2.38 

1.0 

79182 

11.4 

2.38 

1.0 

79213 

11.4 

2.38 

1.0 

79244 

11.4 

2.38 

1.0 

79274 

11.4 

2.38 

1.0 

79305 

11.4 

2.38 

1.0 

79335 

11.4 

2.38 
80365 

1.0 

13 

01 

365 

450. 

5000. 

5. 

80001 

11.4 

2.38 

1.0 

80032 

11.4 

2.38 

1.0 

80060 

11.4 

2.38 

1.0 

80091 

11.4 

2.38 

1.0 

80121 

11.4 

2.38 

1.0 

80130 

54.0 

26.9 

1.0 

80152 

11.4 

2.38 

1.0 

80182 

11.4 

2.38 

1.0 

SEZ  DAIRY 

NUTRIENT  SIMULATION 

0      1 
,0003  .0001   0.10   0.10   99.0   99.0 


99. 


45. 


45. 


180 


90, 


100, 


180 


90, 


100 


(continued  on  next  page) 


Table  C-3  continued. 
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80213 

11.4 

2.38 

1.0 

80244 

11.4 

2.38 

1.0 

80274 

11.4 

2.38 

1.0 

11.4 

2.38 

1.0 

80335 

11.4 

2.38 
81365 

1.0 

13 

01 

365 

450. 

5000. 

5. 

81001 

11.4 

2.38 

1.0 

(siim 

lar  to 
ti 

1980) 

81335 

11.4 

2.38 
82365 

1.0 

13 

01 

365 

450. 

5000. 

5. 

82001 

11.4 

2.38 

II 

1.0 

(sirm 

lar  to 

1980) 

82335 

11.4 

2.38 

1.0 

45, 


45, 


180 


90.   100, 


180 


90.   100. 


APPENDIX  D 
BASIN  MODEL  FILES 


Table  D-l.  BASIN  parameter  file  for  the  TCNS  simulation. 


File:  BASIN. PARM 


1 

1 

47   10    0 

20 

118 

85 

0.00 

3.25 

0.50 

0.75  1.00 

1.25 

1.50 

1.75 

2.00 

2.25 

2.50  . 

3.00 

3.50 

4.00  4.50 

5.00 

5.50 

6.00 

7.00 

8.00 

9.00  10.00 

11.00 

12.00  14.00 

16.00 

18.00 

20.00 

22.00 

25.00 

28.0 

31.0 

34.0 

43.0  48.0 

56.0 

60.0 

67.0 

87.0 

95.0 

105.0  120.0 

138.0 

171.0  188.0 

205.0 

240.0 

298 

302 

301.1010 

210 

273 

301.1020 

220 

257 

301.1030 

212 

232 

301.1040 

220 

256 

301.1050 

255 

269 

301.2010 

261 

268 

302.0000 

292 

297 

303.0000 

298 

302 

304.0000 

310 

327 

305.0000 

314 

333 
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Table  0-2.  Average  annual  cell  runoff,  nitrogen  and  phosphorus  loads 
determined  by  CREAMS-WT  simulations,  and  input  to  the 
BASIN  model  for  simulating  the  TCNS  basin. 


ACCESS 

NITROGEN 

PHOSPHORUS 

RUNOFF 

CODE 

(kg/ha) 

(kg/ha) 

(cm) 

13101 

6.232 

6.032 

21.687 

13201 

6.410 

6.406 

21.920 

13301 

6.809 

6.906 

21.987 

14101 

6.232 

6.032 

21.687 

14201 

6.410 

6.406 

21.920 

14301 

6.809 

6.906 

21.987 

15105 

5.515 

2.964 

21.687 

15106 

5.697 

3.009 

21.687 

15107 

5.879 

3.054 

21.687 

15108 

6.062 

3.108 

21.687 

15109 

6.246 

3.152 

21.687 

15110 

6.430 

3.206 

21.687 

15111 

6.622 

3.251 

21.687 

15112 

6.996 

3.350 

21.687 

15113 

7.371 

3.448 

21.687 

15114 

7.755 

3.547 

21.687 

15115 

8.134 

3.636 

21.687 

15205 

5.719 

3.149 

21.920 

15206 

5.906 

3.196 

21.920 

15207 

6.093 

3.243 

21.920 

15208 

6.281 

3.299 

21.920 

15209 

6.469 

3.346 

21.920 

15210 

6.658 

3.402 

21.920 

15211 

6.854 

3.449 

21.920 

15212 

7.235 

3.553 

21.920 

15213 

7.618 

3.656 

21.920 

15214 

8.009 

3.759 

21.920 

15215 

8.397 

3.853 

21.920 

15305 

6.003 

3.343 

21.987 

15306 

6.200 

3.392 

21.987 

15307 

6.397 

3.442 

21.987 

15308 

6.595 

3.501 

21.987 

15309 

6.793 

3.550 

21.987 

15310 

6.992 

3.610 

21.987 

15311 

7.201 

3.659 

21.987 

15312 

7.608 

3.768 

21.987 

15313 

8.016 

3.877 

21.987 

15314 

8.433 

3.986 

21.987 

15315 

8.844 

4.085 

21.987 
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Table  D-2  continued. 


ACCESS   NITROGEN  PHOSPHORUS  RUNOFF 
CODE     (kg/ha)   (kg/ha)    (cm) 


16104 

5.155 

1.982 

25.509 

16105 

5.327 

2.032 

25.509 

16106 

5.501 

2.082 

25.509 

16107 

5.674 

2.132 

25.509 

16108 

5.849 

2.182 

25.509 

16109 

6.024 

2.222 

25.509 

16110 

6.200 

2.272 

25.509 

16111 

6.376 

2.322 

25.509 

16112 

6.731 

2.412 

25.509 

16113 

7.088 

2.512 

25.509 

16114 

7.447 

2.612 

25.509 

16115 

7.808 

2.702 

25.509 

16116 

8.168 

2.802 

25.509 

16117 

8.530 

2.892 

25.509 

16118 

8.892 

2.992 

25.509 

16204 

5.291 

2.101 

25.203 

16205 

5.473 

2.153 

25.203 

16206 

5.657 

2.205 

25.203 

16207 

5.841 

2.257 

25.203 

16208 

6.027 

2.309 

25.203 

16209 

6.213 

2.350 

25.203 

16210 

6.400 

2.402 

25.203 

16211 

6.587 

2.454 

25.203 

16212 

6.961 

2.548 

25.203 

16213 

7.336 

2.652 

25.203 

16214 

7.712 

2.756 

25.203 

16215 

8.089 

2.850 

25.203 

16304 

5.521 

2.249 

25.485 

16305 

5.713 

2.305 

25.485 

16306 

5.906 

2.360 

25.485 

16307 

6.100 

2.415 

25.485 

16308 

6.295 

2.471 

25.485 

16309 

6.491 

2.515 

25.485 

16310 

6.687 

2.570 

25.485 

16311 

6.884 

2.626 

25.485 

16312 

7.278 

2.725 

25.485 

16313 

7.673 

2.836 

25.485 

16314 

8.068 

2.947 

25.485 

16315 

8.465 

3.047 

25.485 

16404 

5.713 

2.305 

25.485 

17302 

5.069 

0.945 

32.421 

17303 

5.273 

0.999 

32.421 

17304 

5.469 

1.054 

32.421 

17305 

5.676 

1.122 

32.421 
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Table  D-2  continued. 


ACCESS   NITROGEN  PHOSPHORUS  RUNOFF 
CODE    (kg/ha)   (kg/ha)    (cm) 


17306 

5.882 

1.176 

32.421 

17307 

6.090 

1.230 

32.421 

17308 

6.297 

1.285 

32.421 

17309 

6.504 

1.339 

32.421 

17402 

5.017 

0.942 

32.385 

17403 

5.256 

1.002 

32.385 

17404 

5.483 

1.063 

32.385 

17405 

5.723 

1.139 

32.385 

17406 

5.964 

1.199 

32.385 

17407 

6.204 

1.260 

32.385 

17408 

6.445 

1.321 

32.385 

17409 

6.686 

1.381 

32.385 

18433 

181.520 

34.952 

56.079 

18434 

228.137 

43.989 

56.079 

18435 

254.040 

49.004 

56.079 

18436 

295.481 

57.035 

56.079 

18437 

316.196 

61.044 

56.079 

18438 

352.450 

68.070 

56.079 

18439 

456.049 

88.128 

56.079 

18440 

497.490 

96.160 

56.079 

18441 

549.284 

106.202 

56.079 

18442 

626.980 

121.246 

56.079 

18443 

720.216 

139.306 

56.079 

18444 

891.146 

172.410 

56.079 

18445 

979.206 

189.479 

56.079 

18446 

1067.251 

206.533 

56.079 

18447 

1248.550 

241.649 

56.079 

35105 

5.770 

2.170 

21.687 

35106 

6.432 

2.340 

21.687 

35107 

7.105 

2.510 

21.687 

35108 

7.774 

2.681 

21.687 

35109 

8.445 

2.851 

21.687 

35110 

9.117 

3.021 

21.687 

35111 

9.791 

3.191 

21.687 

35112 

11.146 

3.532 

21.687 

35113 

12.505 

3.863 

21.687 

35114 

13.866 

4.204 

21.687 

35115 

15.236 

4.544 

21.687 

35116 

16.602 

4.885 

21.687 

35117 

17.970 

5.225 

21.687 

35118 

19.338 

5.566 

21.687 

35119 

22.079 

6.247 

21.687 

35120 

24.827 

6.928 

21.687 

35121 

27.569 

7.600 

21.687 
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Table  D-2  continued, 


ACCESS 

NITROGEN 

PHOSPHORUS 

RUNOFF 

CODE 

(kg/ha) 

(kg/ha) 

(cm) 

35122 

30.311 

8.281 

21.687 

35123 

33.060 

8.962 

21.687 

35124 

35.803 

9.643 

21.687 

35125 

41.295 

11.005 

21.687 

35126 

46.781 

12.358 

21.687 

35127 

52.274 

13.720 

21.687 

35128 

57.767 

15.082 

21.687 

35129 

63.252 

16.435 

21.687 

35130 

71.489 

18.478 

21.687 

35131 

79.725 

20.512 

21.687 

35132 

87.961 

22.555 

21.687 

35133 

96.198 

24.589 

21.687 

35134 

120.908 

30.710 

21.687 

35135 

134.638 

34.106 

21.687 

35136 

156.605 

39.545 

21.687 

35210 

9.227 

3.106 

21.920 

35211 

9.924 

3.284 

21.920 

35212 

11.323 

3.641 

21.920 

35213 

12.726 

3.988 

21.920 

35214 

14.134 

4.345 

21.920 

35215 

15.554 

4.701 

21.920 

35216 

16.970 

5.058 

21.920 

35217 

18.387 

5.414 

21.920 

35218 

19.806 

5.771 

21.920 

35219 

22.650 

6.484 

21.920 

35220 

25.505 

7.197 

21.920 

35221 

28.357 

7.901 

21.920 

35222 

31.212 

8.614 

21.920 

35223 

34.075 

9.327 

21.920 

35224 

36.932 

10.040 

21.920 

35225 

42.652 

11.466 

21.920 

35226 

48.367 

12.883 

21.920 

35227 

54.090 

14.309 

21.920 

35228 

59.813 

15.735 

21.920 

35229 

65.528 

17.152 

21.920 

35230 

74.109 

19.291 

21.920 

35231 

82.690 

21.421 

21.920 

35232 

91.271 

23.560 

21.920 

35233 

99.852 

25.690 

21.920 

35234 

125.598 

32.098 

21.920 

35235 

139.903 

35.654 

21.920 

35236 

162.792 

41.350 

21.920 

35310 

9.558 

3.196 

21.987 

35311 

10.291 

3.384 

21.987 
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Table  D-2  continued. 


ACCESS   NITROGEN  PHOSPHORUS  RUNOFF 
CODE     (kg/ha)   (kg/ha)    (cm) 


35312 

11.762 

3.760 

21.987 

35313 

13.237 

4.126 

21.987 

35314 

14.715 

4.502 

21.987 

35315 

16.203 

4.879 

21.987 

35316 

17.685 

5.255 

21.987 

35317 

19.169 

5.632 

21.987 

35318 

20.653 

6.008 

21.987 

35319 

23.625 

6.760 

21.987 

35320 

26.605 

7.513 

21.987 

35321 

29.578 

8.256 

21.987 

35322 

32.551 

9.008 

21.987 

35323 

35.532 

9.761 

21.987 

35324 

38.506 

10.514 

21.987 

35325 

44.461 

12.019 

21.987 

35326 

50.409 

13.514 

21.987 

35327 

56.364 

15.020 

21.987 

35328 

62.320 

16.525 

21.987 

35329 

68.268 

18.020 

21.987 

35330 

77.198 

20.278 

21.987 

35331 

86.129 

22.526 

21.987 

35332 

95.059 

24.784 

21.987 

35333 

103.990 

27.032 

21.987 

35334 

130.782 

33.795 

21.987 

35335 

145.670 

37.549 

21.987 

35336 

169.488 

43.560 

21.987 

36108 

7.088 

2.512 

25.509 

36109 

7.447 

2.612 

25.509 

36110 

7.808 

2.702 

25.509 

36111 

8.168 

2.802 

25.509 

36112 

8.892 

2.992 

25.509 

36113 

9.618 

3.182 

25.509 

36114 

10.345 

3.372 

25.509 

36115 

11.074 

3.561 

25.509 

36116 

11.805 

3.761 

25.509 

36117 

12.537 

3.951 

25.509 

36118 

13.269 

4.141 

25.509 

36208 

7.336 

2.652 

25.203 

36209 

7.712 

2.756 

25.203 

36210 

8.089 

2.850 

25.203 

36211 

8.466 

2.953 

25.203 

36212 

9.223 

3.151 

25.203 

36213 

9.981 

3.349 

25.203 

36214 

10.741 

3.546 

25.203 

36215 

11.502 

3.744 

25.203 
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Table  D-2  continued, 


ACCESS   NITROGEN  PHOSPHORUS  RUNOFF 
CODE    (kg/ha)   (kg/ha)    (cm) 


36216 

12.264 

3.952 

25.203 

36217 

13.027 

4.149 

25.203 

36218 

13.791 

4.347 

25.203 

36308 

7.673 

2.836 

25.485 

36309 

8.068 

2.947 

25.485 

36310 

8.465 

3.047 

25.485 

36311 

8.862 

3.157 

25.485 

36312 

9.658 

3.368 

25.485 

36313 

10.456 

3.578 

25.485 

36314 

11.254 

3.789 

25.485 

36315 

12.054 

3.999 

25.485 

36316 

12.855 

4.220 

25.485 

36317 

13.657 

4.431 

25.485 

36318 

14.459 

4.641 

25.485 

37302 

6.044 

1.660 

32.421 

37303 

6.249 

1.715 

32.421 

37304 

6.445 

1.769 

32.421 

37305 

6.652 

1.837 

32.421 

37306 

6.858 

1.891 

32.421 

37307 

7.065 

1.946 

32.421 

37308 

7.273 

2.000 

32.421 

37309 

7.480 

2.054 

32.421 

37402 

6.492 

1.909 

32.385 

37403 

6.731 

1.970 

32.385 

37404 

6.959 

2.030 

32.385 

37405 

7.199 

2.106 

32.385 

37406 

7.439 

2.167 

32.385 

37407 

7.680 

2.227 

32.385 

37408 

7.920 

2.288 

32.385 

37409 

8.162 

2.349 

32.385 
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Table  D-3.  CREAMS-WT  hydrology  parameter  files  used  for  TCNS  cell 
simulations. 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:   131—,  141—,  151- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

80. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.50 

35 

.50 

125 

1.00 

225 

1.00 

335 

.80 

366 

.50 

0 

0 

0 

(d- 

tto  18 

times) 

-1 

0 

0 

81.3   80.7   75.1 


447.   405, 


370, 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:   132—,  142—,  152- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

82. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.50 

35 

.50 

125 

1.00 

225 

1.00 

335 

.80 

366 

.50 

0 

0 

0 

(di 

tto  18 

times) 

-1 

0 

0 

81.3   80.7   75.1 


447 


405, 


370, 
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Table  D-3  continued. 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:  161—,  181- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

80. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.44 

35 

.44 

125 

0.88 

225 

0.88 

335 

.70 

366 

.44 

0 

0 

0 

(di 

tto  18 

times) 

-1 

0 

0 

81.3   80.7   75.1 


447 


405. 


370. 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:   143—,  153- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

84. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.50 

35 

.50 

125 

1.00 

225 

1.00 

335 

.80 

366 

.50 

0 

0 

0 

(d- 

tto  18 

times) 

-1 

0 

0 

81.3   80.7   75.1 


447, 


405, 


370, 
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Table  D-3  continued. 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:  162- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

84. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.44 

35 

.44 

125 

0.88 

225 

0.88 

335 

.70 

366 

.44 

0 

0 

0 

(di 

tto  18 

times) 

-1 

0 

0 

81.3 


447 


,7   75.1 


405.   370, 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:  163-- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

86. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.44 

35 

.44 

125 

0.88 

225 

0.88 

335 

.70 

366 

.44 

0 

0 

0 

(d" 

tto  18 

times) 

-1 

0 

0 

81.3 


447. 


405, 


75.1 


370. 
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Table  D-3  continued. 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:   173- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

86. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.35 

35 

.35 

125 

0.70 

225 

0.70 

335 

.56 

366 

.35 

0 

0 

0 

(d- 

tto  18 

times) 

-1 

0 

0 

81.3   80.7   75.1 


447 


405, 


370. 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:   174— 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

90. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.35 

35 

.35 

125 

0.70 

225 

0.70 

335 

.56 

366 

.35 

0 

0 

0 

(d- 

tto  18 

times) 

-1 

0 

0 

81.3 


447.   405, 


75.1 


370. 


174 


Table  D-3  continued, 


CELL  HYDROLOGY  SIMULATION  FOR  ACCESS  CODES:  184- 


60000 

0 

1 

1 

12.36 

1.00 

.25 

.50 

3.3 

.43 

0.03 

.20 

90. 

.010 

2.0 

30.0 

.056 

.277 

0.73 

0.84 

.847 

1.00 

1.00 

60.2 

61.2 

66.4 

74.2 

76.1 

80.8 

82.7 

69.5 

65.1 

282. 

340. 

436. 

510. 

507. 

472. 

452. 

302. 

281. 

0.5 

1 

.10 

35 

.10 

125 

0.20 

225 

0.20 

335 

.12 

366 

.10 

0 

0 

0 

(d- 

tto  18 

times) 

-1 

0 

0 

81.3   80.7   75.1 
447.   405.   370. 
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Table  D-4.  CREAMS-WT  nutrient  parameter  files  (without  fertilizer 
parameters)  used  for  TCNS  cell  simulations. 


15  =  LANDUSE;  Access  code:  151 — 
0.0 

60000  0     0     0      1 

.40  0.27    1.0 

2 

.20  .200   15.0  .0003  .0001   0.12   0.12   99.0   99.0   99, 

99.0  1.0 
60365 

6  01    365 

450.  8000.     5.    45.    180    90.   100. 


15  =  LANDUSE;  Access  code:  152 — 
0.0 

60000  0     0     0      1 

.40  0.27    1.0 

2 

.20  .200   15.0  .0003  .0001   0.12   0.12   99.0   99.0   99. 

99.0  1.0 
60365 

6  01    365 

450.  8000.     5.    45.    180    90.   130. 


16  =  LANDUSE;  Access  code:  161 
0.0 

60000  0     0     0     1 

.40  0.27    1.0 

2 

.15  .150   15.0  .0003  .0001   0.10   0.10   99.0   99.0   99, 

99.0  1.0 
60365 

6  01    365 

450.  5000.     5.    45.    180    90.    85. 
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Table  D-4  continued, 


15  =  LANDUSE;  Access  code:  153 — 
0.0 

60000  0     0     0      1 

.40  0.27    1.0 

2 

.20  .200   15.0  .0003  .0001   0.12   0.12   99.0   99.0   99, 

99.0  1.0 
60365 

6  01    365 

450.  8000.     5.    45.    180    90.   100. 


16  =  LANDUSE;  Access  code:  162 — 
0.0 

60000  0     0     0      1 

.40  0.27    1.0 

2 

.15  .150   15.0  .0003  .0001   0.10   0.10   99.0   99.0   99, 

99.0  1.0 
60365 

6  01    365 

450.  5000.     5.    45.    180    90.    64. 


16  =  LANDUSE;  Access  code:  163 
0.0 

60000  0     0     0     1 

.40  0.27    1.0 

2 

.15  .150   15.0  .0003  .0001   0.10   0.10   99.0   99.0   99, 

99.0  1.0 
60365 

6  01    365 

450.  5000.     5.    45.    180    90.    64. 
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Table  D-4  continued. 


17  =  LANDUSE;  Access  code:  173 
0.0 


60000  0     0     0     1 

.40  0.27    1.0 

2 

.10  .100   15.0  .0003  .0001   0.08   0.08   99.0   99.0   99. 

99.0  1.0 
60365 

6  01    365 

450.  3500.     5.    45.    180    90.    42. 


17  =  LANDUSE;  Access  code:  174 
0.0 

60000  0     0     0     1 

.40  0.27    1.0 

2 

.10  .100   15.0  .0003  .0001   0.08   0.08   99.0   99.0   99, 

99.0  1.0 
60365 

6  01    365 

450.  3500.     5.    45.    180    90.    42. 


18  =  LANDUSE;  Access  code:  184 
0.0 

60000  0     0     0     1 

.40  0.27    1.0 

2 

.50  .200   60.0  .0003  .0001   0.10   0.10   99.0   99.0   99. 

99.0  1.0 
60365 

6  1    365 

100.  90.     1.    90.    180    90.     5. 
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Table  D-4  continued. 


Fertilization  for  row  crop  and  citrus  included 

13  =  LANDUSE;  Access  code:  131--,  141  — 
0.0 


60000 

0 

0 

0 

1 

.40 

0.27 

1.0 

2 

.25 

.250 

15.0 

.0003 

.0001 

0.12 

0.12 

99.0 

1.0 
60365 

3 

01 

365 

450. 

30000. 

1. 

55. 

180 

90. 

150 

50 

50.0 

33.0 

1.0 

150 

50.0 

33.0 

1.0 

250 

50.0 

33.0 
61365 

1.0 

3 

01 

365 

450. 

30000. 

1. 

55. 

180 

90. 

150 

50 

50.0 

33.0 

1.0 

150 

50.0 

33.0 

1.0 

250 

50.0 

33.0 
62365 

1.0 

(fer 

tilizer 
79365 

loads 

repeat 

for  each 

year, 

62-78 

3 

01 

365 

450. 

30000. 

1. 

55. 

180 

90. 

150 

50 

50.0 

33.0 

1.0 

150 

50.0 

33.0 

1.0 

250 

50.0 

33.0 
0 

1.0 

99.0   99.0   99, 
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Table  D-5.  BASIN  model  output  for  Taylor  Creek  —  Nubbin  Slough 
simulation. 


SUMMARY  OF  NUTRIENT  LOADS  AT  SUB-BASIN  OUTLETS 
ANNUAL  AVERAGE  OVER  20  YEAR  PERIOD 

Average  Annual  Rainfall  (cm)  :    118.85 
Average  Annual  Runoff    (cm)  :     28.32 

CELL  Area  =   5.0  ha 


SUBSHED  OUTLET  #  OF  FLOW  NITROGEN  PHOSPHORUS 

CODE  ROW  COL  CELLS  (cm)  (kg/ha)  (mg/1)  (kg/ha)  (mg/1) 

301  1010  210  273  195  27.9  20.69  7.41  5.62  2.01 

301  1020  220  257  256  28.3  25.41  8.98  6.99  2.47 

301  1030  212  232  866  29.5  6.90  2.34  2.08  0.71 

301  1040  220  256  1185  28.4  6.43  2.26  2.08  0.73 

301  1050  255  269  1697  27.4  6.79  2.48  1.99  0.73 

301  2010  261  268  1201  29.2  6.85  2.35  1.84  0.63 

302  0  292  297  1157  28.7  17.00  5.92  4.75  1.65 

303  0  298  302  1120  28.1  20.91  7.44  5.68  2.02 

304  0  310  327  363  27.8  21.07  7.58  5.33  1.92 

305  0  314  333  1240  27.9  12.07  4.32  3.21  1.15 

Basin  Outlet. . . 

298  302  9280  28.3  9.04  3.19  3.29  1.16 
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Table  D-5  continued. 


NUTRIENT  BUDGET  SUMMARY 

(Load  units:  kg/5  ) 

Summary  by  Sub-watershed:  Nitrogen 

---  LOADS  — -  LOAD  REDUCTION Bkgrnd  Outlet  Basin 

Code   Cell  Stream  LotRet  Impdmt  Ovrlnd  Stream  Load  Load  Delivry 

1  4223     457      0  0    102    568  25  4034  1379 

2  4323    4527      0  0    279   2124  58  6505  2100 

3  4972    3397      0  0   1003   1579  184  5971  3902 

4  6361    1611      0  0    336   1131  1117  7622  5282 

5  5494    6528     0  0    402   2845  2742  11517  8394 

6  4718    5166     0  0    670   3022  2041  8232  6257 

7  6979   22745     0  0    467  10922  1331  19666  17697 

8  8916   25935     0  0    305  11973  844  23416  23416 

9  3307    5065     0  0     55    995  327  7649  5198 
10    4824   14037      0  0    377   5436  1919  14967  10247 


Summary  by  Sub-watershed:  Phosphorus 

—  LOADS  LOAD  REDUCTION  Bkgrnd  Outlet  Basin 

Code  Cell  Stream  LotRet  Impdmt  Ovrlnd  Stream  Load  Load  Delivry 

1  1085  87  0  0  9  72  3  1095  1095 

2  1197  858  0  0  25  249  8  1789  1789 

3  1502  652  0  0  92  285  26  1803  1803 

4  2340  312  0  0  30  311  159  2470  2470 

5  2354  1251  0  0  37  582  391  3377  3377 

6  1423  1006  0  0  60  453  291  2207  2207 

7  2219  4334  0  0  39  1212  190  5491  5491 

8  2617  4960  0  0  25  1316  120  6356  6356 

9  1062  967  0  0  5  138  46  1933  1933 
10  1739  2690      0  0  33  685  274  3985  3985 
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Table  D-5  continued, 


Summary  by  Land  Use:  Nitrogen 

---  LOADS  ---  LOAD  REDUCTION  - -  Bkgrnd  Outlet 

Code  Cell  Stream  LotRet  Impndmt  Overlnd  Stream  Load  Load 

13  897  0  0  0      0  126  0  771 

14  903  0  0  0      0  155  0  748 

15  1454  2016  0  0      0  1371  0  2099 

16  15580  17495  0  0      0  9560  0  23515 

17  2353  1475  0  0      0  468  0  3359 

18  12984  0  0  0      0  3785  0  9198 

19  0  0000000 

20  0  0000000 

21  0  0  0  0      0  0  10592  10592 

22  0  526  0  0     121  53  0  352 

23  0  4330  0  0     501  1073  0  2755 

24  0  4431  0  0    1263  1320  0  1847 

25  0  6733  0  0    1763  1892  0  3076 

26  0  873  0  0     167  241  0  465 

27  0  584  0  0     184  88  0  311 

35  4864  22907  0  0      0  8405  0  19365 

36  2082  13608  0  0      0  4713  0  10977 

37  13000  14490  0  0      0  7344  0  20146 


Summary  by  Land  Use:   Phosphorus 

—  LOADS  LOAD  REDUCTION  Bkgrnd  Outlet 

Code  Cell  Stream  LotRet  Impndmt  Overlnd  Stream  Load  Load 

13  869  0  0      0  0  98  0  770 

14  894  0  0      0  0  124  0  770 

15  737  374  0      0  0  224  0  887 

16  6152  3353  0      0  0  1643  0  7862 

17  441  295  0      0  0  56  0  680 

18  2511  0  0      0  0  340  0  2170 

19  0  0000000 

20  0  0000000 

21  0  0  0               0  0  0  1513  1513 

22  0  101  0               0  10  5  0  84 

23  0  822  0               0  43  105  0  673 

24  0  849  0      0  113  161  0  575 

25  0  1286  0      0  160  226  0  899 

26  0  167  0      0  14  26  0  126 

27  0  112  0      0  16  10  0  85 

35  1420  4254  0      0  0  803  0  4871 

36  725  2608  0      0  0  472  0  2861 

37  3790  2898  0      0  0  1008  0  5679 
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Table  D-6.  BASIN  model  output  from  simulation  of  a  "maximum"  BMP 

scenario  applied  to  the  Taylor  Creek--Nubbin  Slough  basin, 


SUMMARY  OF  NUTRIENT  LOADS  AT  SUB-BASIN  OUTLETS 
ANNUAL  AVERAGE  OVER  20  YEAR  PERIOD 

Average  Annual  Rainfall  (cm)  :    118.85 
Average  Annual  Runoff    (cm)  :     28.24 

CELL  Area  =   5.0  ha 


SUBSHED  OUTLET  #  OF  FLOW  NITROGEN  PHOSPHORUS 

CODE  ROW  COL  CELLS  (cm)   (kg/ha)  (mg/1)   (kg/ha)  (mg/1) 

301  1010  210  273  195  27.3  10.64  3.90  3.46  1.27 

301  1020  220  257  256  28.0  9.46  3.38  3.35  1.20 

301  1030  212  232  866  29.5  5.51  1.87  1.66  0.56 

301  1040  220  256  1185  28.4  5.45  1.92  1.81  0.64 

301  1050  255  269  1697  27.4  4.51  1.65  1.40  0.51 

301  2010  261  268  1201  29.2  4.83  1.66  1.22  0.42 

302  0  292  297  1157  28.7  5.43  1.89  1.65  0.57 

303  0  298  302  1120  27.9  5.79  2.07  1.82  0.65 

304  0  310  327  363  27.7  6.87  2.48  2.35  0.85 

305  0  314  333  1240  27.8  4.77  1.71  1.38  0.50 

Basin  Outlet.. . 

298  302  9280  28.2  4.67  1.65  1.67  0.59 
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Table  D-6  continued. 


NUTRIENT  BUDGET  SUMMARY 

(Load  units:  kg/5  ) 

Summary  by  Sub-watershed:  Nitrogen 

-—  LOADS  ---  LOAD  REDUCTION Bkgrnd  Outlet  Basin 

Code   Cell  Stream  LotRet  Impdmt  Ovrlnd  Stream  Load  Load  Delivry 

1  4335     112   7338  166     45    249  25  2075  995 

2  4550     290   4724  370    121    740  58  2421  1283 

3  5173    1326   1531  221    437    852  184  4767  3734 

4  6613     458   1748  330    126    810  1117  6459  5049 

5  5714     630     0  283    155   1002  2742  7647  6935 

6  4917     915   2803  257    191    886  2041  5797  5264 

7  7553     422   4738  606     67   1101  1331  6281  6066 

8  9644     284   8062  734     35   1389  844  6484  6484 

9  3533     81   3481  256     14    259  327  2492  2066 
10    5161     484   1777  347     95    733  1919  5919  5412 


Summary  by  Sub-watershed:  Phosphorus 

—  LOADS  - LOAD  REDUCTION  Bkgrnd  Outlet  Basin 

Code  Cell  Stream  LotRet  Impdmt  Ovrlnd  Stream  Load  Load  Delivry 

1  1114  21   1419  42  4  45  3  673  673 

2  1255  55    913  78  11  130  8  858  858 

3  1557  253    296  75  42  204  26  1437  1437 

4  2409  87    338  131  12  273  159  2150  2150 

5  2415  120  0  140  16  403  391  2369  2369 

6  1478  175    541  78  20  236  291  1467  1467 

7  2371  80    916  181  7  307  190  1904  1904 

8  2809  54   1558  193  3  337  120  2038  2038 

9  1121  15    673  74  1  77  46  853  853 
10  1829  93    343  131  9  252  274  1713  1713 
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Table  D-6  continued. 


Summary  by  Land  Use:   Nitrogen 

---  LOADS  LOAD  REDUCTION  Bkgrnd  Outlet 

Code  Cell  Stream  LotRet  Impndmt  Overlnd  Stream  Load  Load 

13  897  0  0     102      0  98  0  696 

14  903  0  0     107      0  120  0  675 

15  1562  0  0     167      0  288  0  1107 

16  16344  0  0     805      0  2344  0  13194 

17  2367  0  0      17      0  73  0  2275 

18  12984  0   36205  0      0  999  0  2428 

19  0  0000000 

20  0  0000000 

21  0  0  0      0      0  0  10592  10592 

22  0  124  0      0     43  4  0  76 

23  0  549  0      0     90  35  0  424 

24  0  1566  0      0     422  294  0  848 

25  0  2341  0      0     642  328  0  1371 

26  0  211  0      0     30  23  0  157 

27  0  213  0      0     61  28  0  123 

35  5806  0  0    1155      0  1179  0  3471 

36  2567  0  0     387      0  412  0  1768 

37  13763  0  0     834      0  1794  0  11134 


Summary  by  Land  Use:  Phosphorus 

---  LOADS  ---  LOAD  REDUCTION  Bkgrnd  Outlet 

Code  Cell   Stream  LotRet  Impndmt  Overlnd  Stream  Load  Load 

13  869  0  0  96  0  86  0  685 

14  894  0  0  101  0  109  0  683 

15  766  0  0  70  0  119  0  576 

16  6367  0  0  339  0  959  0  5068 

17  445  0  0  5  0  20  0  418 

18  2511  0  7002  0  0  89  0  573 

19  0  0000000 

20  0  0000000 

21  0  0  0  0  0  0  1513  1513 

22  0  23  0  0  4  1  0  18 

23  0  105  0  0  9  8  0  87 

24  0  300  0  0  41  42  0  216 

25  0  448  0  0  64  54  0  329 

26  0  40  0  0  3  3  0  33 

27  0  40  0  0  5  3  0  31 

35  1654  0  0  192  0  194  0  1268 

36  852  0  0  88  0  92  0  671 

37  4002  0  0  230  0  480  0  3291 
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